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ABSTRACT 


inewadiation wheranpyes theres! are. sp manys? iustances: (yin 
which ipradiation of..the,,-lung is unavoidable. In 
particular, imelarce; fields irradiations, with, h—rayss « there 
is a potential for severe respiratory complications which 
mays bes fatal...) The occurence, and, severitys of, *= dung. damage 
ist telateds§ tojagthe avolume of \lune,iwtradiaited, therytotal 
dose, the dose-rate,. the. overall (durationrseot. ~radiation 
therapy, as well as the variable physiological response of 
individual patients. 

The present research deals with radiation—-induced 
pneumonitis in two ways: 

(1) A non-invasive diagnostic test based on X-ray computed 
Lomocrapnya CCl» ae quantifying lung damage early after 
irradiation, and 

(2)ieimproved: accuracy, )of, lung «dose, caleulations:) for 
planning radiation therapy. 

Theareactiony of the, lune. to. radiationsds monitored. dn 
groups of animals which have been irradiated to various 
single doses to the thorax. This test was able to detect 
the radiation damage, its time course, and dose response 
in a large number of mice. In order to compare this test 


tow otner “in-vivo sdiagnostic tests, such as) conventional 


X-radiography and nuclear medicine studies, however, 
larger “animals Gil.e. dogs) were subsequently used. To 
ficiher evaluates ther techniques ing a clinical, setting, the 
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densitometry test is applied to radiotherapy patients at 
different times before and after radiation therapy. 

Once the respiratory distress of radiation 
pneumonitis occurs, treatment® is “often *not™effective. 
However, the treatment of radiation pneumonitis is an area 
of active research. This provides sufficient incentive to 
develop an early quantitative assay of the efficacy (and 
perhaps mechanism) of new treatment methods. With a 
reliable diagnostic test, more timely intervention may 
become possible. 

The second way to deal with radiation pneumonitis is 
to reduce! its incidence’ by Limiting and tcontrolling the 
dose delivered to lung. For this purpose it is necessary 
toe have® an “-acturates*method = of secalcudarting® the dose 
actually: delivered tov lunge. *Thisyswork® allisoe deadisy with 
such radiation dose’ Vcalculations and dosimetry din large 
(half-body) fields. Various lung dose calculation methods 
simple as well as complex, are compared to measured doses 
in simple geometrical phantoms in order to study their 
range of applicability. These methods are also tested in 
more complex situations and in a humanoid phantom exposed 
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N T is the average CT number for each slice 
ahtnonie mid-lung. 


(A) Average CT number (Nc7) for the control 
SHoOUp aor Sil ce sas Va Bianieit Monitot tiime tatiter 
irradiation; (B) Average CT number (Nor) 
fom the *group of i mice irradiated to 7 Cy. 
Corresponding densities are shown on the 
right tordinate:. 


Average CT number (N ¢) for: (A) group of 
mice irradiated to 18 Gy; CB) “oroup of mice 
ivradi@ted to V3: Gy; “(C) group of mice 
irradiated to 14 Gy. 


CT scans of mouse lungs showing radiation- 
induced changes at 7 weeks, 16 weeks, 
18 weeks. 


The slung (density Lor Andivaduale mice 
invadiated torlorGy-e No; is the vaverage 
CT number for each slice through mid-lung. 


The percent of animals whose lung density 
increased by more than 40% is plotted as a 
function of time after irradiation for the 
FLOuUps .0i mice: de tadi ated "to el0! 12 5 enand 

14 Gy. The incidence is defined with reference 
to, ethie Simitital number of tanimals C002) in 
each dose group. 


The percent of animals whose lung density 
increased by more than 40% is plotted as a 
fimctiomfof idose ducine the acute (radiation 
pneumonitis) phase 15-21 weeks. The number of 
animals with severe effects increases with 
dose. 


The average CT number (Nor ) for. an 
individual mouse irradiated to 14 Gy. 


Chiscantonh aemouse Slunetiinradiated to 14 Gy 
CA) 23 Sweekis Gat ter tirradiation 

CBA) (2). JDeweekstlaftve G eirradiation 

The altered shape of the lung is apparent, 
particularly in the pleural region. 


Chtsdanstof Gtwoimicevtirradiated to 10 Gy 
S24 iweews@atter “irradiation GA) mouse 1 
CB) mouwseeaswStructural changes fare 
appace nite Only aii wGAy)s 


Comparison of in-vivo (CT) density and 
volume measurements to ex-situ mass and 
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volume measurements for irradiated 
mouse lung. 


Cis eameotiipa nat ihune CA) with a level 
Settanenor —250°H sand a window setting 

of «130 Ht 7(Bi) wat rthe.-same olevel ‘setting 

but iwith atlargengwindow visetting acf elO000) HH. 


Gheslicesof CA) anormal Gdoey lungs e( Biea 
dog lung after irradiation. An area of 
interest tinom which Na, is caleulated is 
traced. 


Tiesaverages Cl enunbem permsiliice SGN; icin 
dog tluneatonrevanious klong itudinal 
posttions <inSsthevthorax: 


CT slice of dog lung (A) near the apex 
(B) near the diaphragm. 


Average CT number (Nor) for each dog 
imnadiated tito s20aCy,< 


Average CT number (Nor ) forgeach dog 
irradiated to 10 Gy. 


Digital radiograph © scoutsscan ) of.dog 
thorax. The boundaries of the radiation 
field are marked by positioning clips 

on skin marks. 


Thesaverage Ch number.(N-7) for each slice 
in? lune fortiditrierent Longitudinal 
Positions in they thorax 

(A) 1.5 weeks before radiotherapy 

(B) eS “weeks Vafter the end of rnadiotherapy 
(C) 21 weeks after the end of radiotherapy. 


ThéesaveragevCT cnumber®( No, )sforteach slice 
in lung for different longitudinal 
positions in the thorax 

(A) 1 week before radiotherapy 

(B) 5 weeks after the end of radiotherapy 
(COs 20°weekswatter they end.of radiotherapy. 


The average CT number (N,,,) for each slice 
inglungatiortdifiterent _longitudinal 
positions’ in the? thonax 

(A) 1 week before radiotherapy 


(B) 5.5 weeks after the end of radiotherapy. 


Cilescansofehumandlunge(A) normal lung 
GB) tadiatiomwepneumoniitisfineslungs (The 
radiation damage is confined to within 
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the borders of the radiation field. 


Treatment schedule for hemi-body 
radiotherapy patients. 


Average CT number (Nor ) for sindavidual 
patients receiving hemi-body irradiation. 
Raghite and Het teiunee datalare! plotted 
separately. The general decrease is due 
to clearance of disease and absence of 
pneumonitis. 


Variation of dose with depth in a water 
phantom for a parallel-opposed pair of 
Padiations filed dis) G30x30) cm?) .) The 
uniformity improves with energy. 


Bending magnet and treatment head of 
a Siemens Mevatron VI medical linear 
accelerator. 


ochiematices diacramy of an ionization 
chamber with associated electrometer. 
The chamber can be used to measure 
dose rate (current mode) or dose 
(charge mode). 


Schematic of the instrumentation 
necessary to observe thermoluminescence. 


Definitions of tissue-air and 
tissue-phantom ratios. 
TAR=D/ DE, TPR=D/D , 


Geometry used for the Batho correction. 
OP2 and 03 are the electron densities 
relative to water. 


Parameters used in the calculation 
of scatter contribution in the dSAR. 


Definition of inhomogeneity correction 
factor ICF=D,/D, 


Variation of tissue-maximum ratios with 
depth in a phantom. The equivalence of 
these ratios for prestwood, polystyrene 
and water is confirmed. 


Variatdon® of tissue-maximum) ratios with 
depth in a phantom. The equivalence of 
these ratios at distances of 100 cm 

and 200 cm from the source! is» confirmed. 
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The dose increase with depth in shallow 
layers of cork. The depth of the maximum 
doseegtises id max Duis tat 6.5 ecm tfior 
radiation field sizes of 10x10 cm*and 
Poztoren?e 


The inhomogeneity correction for 6 MV 
X-nays  teyplotted  aseaiiunction tof 
parameter “Nalefor field sizev’(A) 5x5. cm; 
(BYSTOzcLOlemiy (OC) 35x35 ecmaslTha. cork 
thicknesses thliwas%.0+0.72 cm. The 
experimental set-up is shown in the 


inset. 


The inhomogeneity correction for 6 MV 
Merays is plotted asa function of 
parameter a’ for field size (A) 5x5 em? 
CB) 1 1LOx0> em: FAC) 95x95 Lem? seThewcork 
thickness, 1t ,fwase®s 2 ecm: 


Geometry for the Batho correction. Dose is 
calculated tatepoine iP, 

(A) when P lies in tissue beyond lung, 

(B) when P lies within lung. 


Geometry used to obtain the experimental 
data. Dose at P is measured as parameter 


a 1g) Qveesl Ql size I T@ ILS ims 


For 6 MV X-rays, the inhomogeneity 
GOrrectPoneractorvtisrplotteduas a function 
Ofs parameter “ay for £ileldy sizer CA) oxo cm 
(B)LOxl0 em §CC) 50x50 cm. 


For Cobalt-60 irradiation, the inhomogeneity 


eOnnpectiom factor is plotted as al function 


Of panpameter, -a (ror field’ size) (A) 5x5 cm 
(BON 1LOx Orem CC) 50x50 cm. 


Spherical coordinate system used for 
calculation of D, and Da . Photons are 
scattered, once at OQ) towards point P. 


GCeonetry used for calculation of ®=tirst— 
scatter dose to point) P from two separate 
layers I and II. 


The inhomogeneity correction for 6 MV X-rays 
at sapoant eo cme below “a 10 em thick cork for 


varying lateral dimensions of the cork. 


The inhomogeneity correction for 6 MV 


X=Laysceassa-tunctions, of parameters at 
forma mediastinum of lo tera) si2en 
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c=5x20 ecm? , and c=10x20 cm’. 


The composite “phantom” used in the 
calculation of the inhomogeneity correction 
Factor proposed) by (Lbulul and Bijarngard (69). 


The inhomogeneity correction factor for 

Ob MVe K=ray st ast ar function of parameter are 
The: field. size is fixed at 40x40 cm*. The 
lateral dimension of the cork, b, is varied 
CA) Sx > <cm- CB): LOx0 ccm? (C.). 20520 tem. 


Radiation profiles across the beam measured 
at several depths in a water phantom for 
GUMY Ri =tavec sate cecoucce-lo-surrace dis tamuce 
of 200 cm. 
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The? doser dtetriabut ion? in’ al CT isan of 
Rando slice 18 is shown. The inhomogeneity 
correction was performed using the 
generalized Batho method (109). 


Pnrewdoserdistributiion, in a-Cl iscan otf 
Rando slice 18 is shown. The inhomogeneity 
correction was performed using the 
equivalent TAR method (110) as implemented 
at the Ontario Cancer Institute. 
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isigRadiation® Pneumond tis 

in@radiatione therapy. there care’) manyoS instances’ in 
which itradtatrdioneio0t “the” Gdunpr is unavoidable For 
example, these include treatments of the lung itself, or 
of the esophagus, ory. ther®breastuut i nitthesel casessn? the 
tolerances ofthe Luneetovradiationi usually llamite= therdose 
which can be prescribed to the diseased site. However, 
Decauserponily al portion Sofh lung o ds oirradiated, severe 
respiratory complications can be avoided. 

There has been a gradual trend towards the use of 


larger?! ori oe magna’ radiation fields which encompass much 
larger volumes of both lungs. Such irra diatom, 
ord ginadily> dimeds ati treating: distributed! lymph* nodeses (as 
imeonodekine’ sitdisease) Wand reducing’ nheriincidencer: of 
metastases, has evolved into half-body and total body 
technigqueste, Halt-body irradiations) (HBL) has >> become ~ a 
standard procedure for treating widespread metastatic 
canceratand for?trelieke of) pat nwtor improverqua lity pol «life 
dari net either remainders “ofthe patientlse late 14100). Thies 
teehniqueldis+- also antvadjuvant form of radical therapy ~ for 
oat Ceuin Pearce 1 nomae as0e we lil asevEwing “st samcomal (128) 
Here; mint, thes CrocseiCancermedinstitutey HBIgcisi’ used “in 
combination with chemotherapy to treat oat cell carcinoma 
(GigZorcagilotal tbody arradiation techniques are).alsol used 
prior to bone marrow transplantation for leukemia patients 


1 


—— 7, a ; 

ss 4 = 7 7 
dizhaownes4 BP age Dh 
ai eeonarent, “iim 25% ived4 eiessi2 odsathes et 
227 efdabiovern @2 facet e440, 3  velsi thea ~aeeee 


. 


‘oe ,¥oeerk . xtel ett) ta tensp2e0%s Fbvise? goat? , SiqQanee 


ad. -S60e2  en5i2 of »SPag7d sit io ~.eupedyoee ~ att te 7 
weob ef3 @32eli. yiseans aatsalbes of garl_s02_ 36 sone ,elez a 
,secewor este baee bee xifs 07 Sets aonerg oi “as Watee 

SJevoe hbasaltieath) »ei “ans *6 iottiog, A eine senagced 


cinwe o¢ ax shedvavtiades Qaerasgqeg= 


to (seth wt? . ebtewW) ys bNes7 fephata 675 3e ‘aut * zedT- 

dene seeqanads iste, ehlsd* og iselhet "eigen" x “segtad 
,7roiselherss Aqué ad wat f3nd in (esmndey wogiat - 
gc) neben dgext Setedizisth @nezast4 38 hanks yileqiyt¢a ; 
ip seacehise! 849.9 anizabex baby Gbeneeeb «'@vagees ~at a 


Yhod-leg02 ban vhed-*2e08 coins y havinys cad ,#9es7e03S= 
e -seoned eed (180) mi Paalpaayt thad-iben se euptedses 
Siragee3ee feetgeablv: -anivests +04 srabssove brabaede: 
atti 26 oye) lawn) eotigat: <2 ndng To *giler vod has a Eee | 
eid? .(001L) stri e's retsey solide! rabnieaex waa ani ma 
ge4 ggeteds Pwotthes ead inne vuLsi ie Ble e) suplods —_ 
«{OS1) amvoree a igeaee $8 fine. foe sere bei - 
ai been of sla : tal (1agva e 


_ 


C645 9:)- 

intlarge titeliddinradiations; Varge “volumes. Sofe<cboth 
hunesigarempitrrnadiated €tol@®alshighetdose, and there isa 
potential for severe respiratory complication which may be 
fatal. The occurence and severity of lung damage is 
relatedeto "thetvolumesofelungcinradiated, the total dose; 
dosesprate ofivadeliveny; overall *adurations Soneradiation 
therapy, as well as the variable physiological response of 
tadividualbeapatients 8014535). "When “eomplications “develop, 
the response of lung to such irradiation has been called 


the acute pneumonitis syndrome” (28,30) and becomes 
evident clinically only several months after exposure _ to 
the radiation. In humans, the symptoms develop after 1 to 
7 months, with—a~peak incidence occuring near 3 £4months. 
bismoceurence ‘asnaaiunction of rdose Cdeliverned lim avsiingle 
fraction) to the entire lung has been documented by Fryer 
fiona. C30) gand® refined ieby (evan DykogetfalieG@h23);. ithe 
clinical symptoms include increasing cough with dyspnea 
and fever. Chest radiographs taken at this time reveal 
opacities. Once these signs have developed, the patient's 
condition deteriorates rapidly and death may ensue within 
Betvorssweeks< In some cases, treatment with antibiotics 
and prednisone diminishes the severity of the radiation 
response, and@eisiterotd twhthdrawalSimay “precipitate: thie 
syndrome. Chemotherapy, when used in combination with 
Gada tikon fn suas ohaneexaceanbha tine Mf acitio“n § (3.5). If Waprhe 


Acwitewsrddiations pueumonditis: ssubsideswor edoes not oceur; 
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fibrests in tunegamay stilidocecurvat tater timesmupito one 
yéar aftere inradidatione($86)s;8in whichecases this becomes 
thewiimnitine@tactor intirnadiationttherapy . 

Mieroscopic fLindingsmineyhumans catacautopsy.« show’ oa 
eombination of the following features: atypical septal 
cell purolifteration, ivascutartchanges ; cand eat later times 
widespread hyaline membrane formation (ll a3 59). Some 
degree of fibrosis may be observed in every lung exposed 
toamtherapeuticy? dosesveort) radiations fatetimes of 69 sto 12 
months after radiotherapy, bit’ Bthetmanet dence] #oingacutée 
radiation pneumonitis is less common. There seems to be 
little correlation between the acute syndrome and_ the 
Mater fLibroticeechanges (14), e"succestingsthatedifterent 
cells may be responsible for early and late damage. A 
good review of damage to different cell types in humans as 
well as animals has been given by Gross (35,38). nate 
lung, the cells with the highest mitotic rates in which 
radiation-induced chromosomal damage is expected to be 
most pronounced are: the bronchial epithelial cells, the 
capillary. endothelial cells, and type=-2 pneumocytes. 
Although these histologic changes in human lung have been 
observed repeatedly at autopsy, no systematic study has 
been (ethically) possible, prepneumonitis changes in human 
lung have certainly not been available. Thus, most of the 
knowledge obtained on radiation reponse of lung has been 
acquired through animal experiments and it is not known 


how well events occuring in animals correspond to events 
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occuring in humans. Indeed yhtherewis indication that the 
timeghcoursesgs,ofaradiationgdamageriniluns noteonly ditfers 
between humans and animals, but that it. alesorndiftters 
within species and even among different strains of mice 
C249) 

Until now the management of radiation pneumonitis has 
consisted of prevention by reducing the dose to lung or 
reducing the volume of lung irradiated. As can be seen 
fromsthesdata of Van. .Dyk.~(128), the incidencewot;radiation 
pneumonitis rises sharply as a function of single dose _ to 
total lune, reaching 1002 for 11 Gy. However, “evaluating 
dose,to lung-aceuratelyjisya difficult. dosimetry problem 
since there are large differences in lung geometry and 
density from =pattent to “patient, “and “there “are “also 
variations of lung density even within a patient 


Gerson. 


l-2 Animal Findings 

In order to better understand the processes’ which 
give rise to radiation pneumonitis, animal models have 
been developed. Histologic “changes pipet mouse lung 
following irradiation are well documented, and are listed 
in Table 1.1. The time course of these changes has_ been 
divided into three periods: acute, intermediate, and late 
(or chronic). The intermediate phase overlaps the acute 
andmlatie phases. JyAccordingstoaelravisry Gli6), thesextenteof 


the histological lesions observed in each phase is 


? rs = . 
5 ; rs 2 
all a | ‘ 


@se?ith cals 24 ~ Fes eee’ /elenine — it Ss 
wotaks sotaris sabeplats grows asye , bee catseqe ln hapie : 
teat 


ead ataanenven ce tinier Der epamegares on? wom fbaae 


sé poet. cz 000! wa gnlewbet ee: syAeepyes to 6a eres. 


ane > al ban = Sf)  begeh beast goes. 20 ou¥is® 2 gatsvber 
ausgathes! tel vanasim® Blm (Rit MAS aby jn étebogay weed 
hi» Geos. 5ighue % ‘Iannat-» eb Zbqiaie snefe Siiironese, _ 
Saljseutavn ,ce°our: eo jf ss) BOUL ancnseae anal faao? : 
SwideGe) arfedbeud:- 21031 128) 6 StViesthee2P soul na wen 
bar 2a ebAae ail! © 2i=ana=sb193h agasl e236 apeas ‘snake 


Bale) Be @2037 34 tnsig@hg ¢9 3081405 weet ¢s teen 


7 
7 

sneiieq «2 b¢n9tv 7899.9 ¥Ptehed tard 1 ene l2Hesey 
otest CP) 7 


agathut® Dumiod fet 
Gd 
Hatdw ser ee yo? eis tal bani 163760. oF gabio ni 
aved afeboc teclirs avieaienting nolvealbe: a4 wei ivi 


at estes oil raciaaiert 900 
| sakin’ 


ee 


desie=dependent; ‘although )°the '<time  ecowrse! sis snot «The 
threshold dose for acute and intermediate effects in CBA 
male mice was determined to be 11 Gy, while the threshold 
dose for late effects was determined to be 13 Gy. hor = aa 
dO0SEM 40,0415  Gytnonialtityioecurred sin tthe aAcutearphase for 
aun andmaits +673 ,106)'.. Ghoweveripcthe radiation. “sensitivity 
of  *nowsie *eluneeimay becditierent uimtdiftiierent tstrains sof 
mucew= elhe sD@5 0/41 60 Randal DiGs0/1G80))i tt or Elhine SW airradtatton 


i Sidiiterent strains Softenmicevhave been teported veo tiie 


Table 1.1: Radiation-induced changes in mouse lung 
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Phase Time Peak Activity Lesions 
(weeks ) (weeks) 
Acute PO = 3.0 G22 Edema 


Macrophage intiltration 
Pabrin® im aa7r space 
Alveolar walls thickened 


Lumina reduced to a slit 


Intermediate 22-54 36 Increased septal cellularity 
Large celissin air tspaces 


and debris 


Late Pio permanent Collagen deposition 


damage Cells diminished in number 
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betweentrieoara yi) OnlaGy aC LOGR1S2).gebeath trates tof tdaninal 
eroupsesis pthe bassayéanddLD C50 /60)vand UD (C50/89) fare tithe 
lethal doses required to achieve a mortality rate of 504 
assayed at 160 and PelLoGy daysiy erespectively.© atter 
imnad tation. £ihorebalbicanicetused iat our lneethtutes (02), 
the LD(50/160) was determined to be 13 Gy, in agreement 
with the above data. Within the acute and intermediate 
phases, there is a period during which the maximum effect 
is-observed. For the acute phase, this peak response 
occurs between 16 and 22 weeks. The observations during 
thvsiepeak “<pertodPricorrespond schosely  yto” those lof ethe 
“tadiation ipneumontiris tsyndrome™ inimanh@0e52)s SoThe speak 
acpi ty for the intermediate phase oecurs at 
approximately 36 weeks. No peak develops for the late 
stageebecauseriibrosismis «ai chronic efifectiresulting “from 
irreversible and cumulative lung damage. 

Detection of some lung damage at times earlier than 
the acute phase was achieved by Sharplin and Franko (102). 
They found an increase in the mass of lung at 6 weeks 
atterwairradiatton, toarkoweredoses tas Lowlas 2 Gy. bheir 
experiments suggest that measurement of lung mass in the 
“prepneumonitis” period may be a more sensitive and 
earlier indicator of lung damage. The mass increase was 
dose-dependént 1eup to sa adosie sof 10 sGy. | fAmerimerease of 182 
in lung mass was detected following irae ddkatirone, wistih, = ya 
dosie@ of HaveGy, Gand €anvincreasie tof 28% was idetected with va 
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Ultra-early changes within 2 weeks after irradiation 
are possible by electron microscopy as discussed for dogs 
by Moosavi et al.(/77). The capillary endothelium appears 
to MpePphehersinitial (site of damage. As in the mouse 
studies,the time course of events is not dose-dependent, 


but the severity of damage is. 


1.38Diagnosticetests 

Various tests can be used to detect radiation damage 
tor Viung; and these are listed in Table 1.2. They are of 
several categories; morphometric, physiological, 
biochentieal “andteradiological. Thescihianges occurring in 
lune,eltstedtin Table 1.1), are all expected to produce an 
increase in density of lung. Therefore measurement of 
changes in lung mass and/or density is expected to be a 
good indicator of lung damage. For this reason, the 
application. Of = toray | computed tomography to early 
detection of pneumonitis was developed in this work (see 
Chaptere2)< 

NoOsernoheethet tests* listed) ing’ tables 1.2. show a 
dose-related increase in severity of the effect, and there 
is’ {aeethreshold™ doserebelow “which) eno, changées®vcan be 
detected. Thiswethreshold dosevyiormdetectionvisenot the 
same for all these tests because a different aspect of 
pneumonitis (i.e. specificity) is being observed and the 
tests have different sensitivity (i.e. inherent 


precision). For example, cellular changes are detected 
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Table WaZiat Testsimseed toy detect® radiation damage to) dung 


Changes in lung Instrumentation Reference 


A. Morphologic 

1. visible changes in cells light microscope, 5 Doo eee Lei 
at histologic and electron microscope 
ultrastructural level 


Bree Pihyisiiomoed c 


1. decrease in lung volume densitometry, lunge mass, 457,102,104 
and compliance spirometry 
2. decreased ventilation 133Xe gas inhalation 86 


Sevescwlar, chaneessisich?, ais 


perfusion reduction nuclear medicine 29% DOs 5.80 
increased permeability alveolar lavage 
Awe enictriedse: ime ibreathins plethysmograph ES) 


frequency at time of 
pneumonitis 


5. changes in arterial arterialy pH, pol. pGco Sor 105 
blood gases 


C. Biochemical 
1. surfactant release radioisotope labelling matey e Moras hiP Ake ite) 
and alveolar lavage 


2. phospholipid increase, radioisotope labelling B6:5,07.,30 
excess of protein in 


alveolar fluid 


3. collagen deposition hydrolyzation 63 
(late) immunofluorescence 76 


D. Physical 


1.2 guantifdcatdons of edema lung weight, densitometry 45,102 
2. changes in density X-ray computed tomography 2 J 2 
Compton tomography NHS eee UE) Ke 


E. Radiological 


i opacities of planar X-rays 65, 
structural changes conventional tomography 
X-ray computed tomography ues) 


Compton tomography 9 
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early and at low doses by electron microscopy, whereas 
radiologicak *changesrsappear wmuchinlatervand for greater 
doses. 

Most of the tests detect the lung damage during the 
acute pneumonitis phase. At this time, treatment is often 
too late and the course of the syndrome cannot be halted 
or reversed. Therefore it might be useful to have 


“prepneumonitis"™ tests that detect lung damage in the 
“latent” period of several weeks before the onset of 
pneumonitis. Such sensitive tests are the biochemical and 
electron microscopy tests which can show changes several 
days after irradiation. Also, lung mass increases’ have 
been observed in mouse as early as 6 weeks after 
irradiation and at low doses (102). 

The testsiewithhivearlkiest detection, (surfactant 
release, electron microscopy, and lung mass increase) are 
unfortunately invasive and not applicable to the 
Management of patients developing radiation pneumonitis. 
The earliest in=vivoytestiwhich can betapplied totpatients 
to monitor lung changes has been described by Prato (86). 
Regtonalsbloodvftowyeing lung «is (assessed byvsr,ingecting 
radioactive Xenon-133 dissolved in saline and monitoring 
activity over the lung using a gamma camera. A_ reduction 
im@sactivityehindieatesea reductiontin ntimberseof particles 
seaching PZunghiwhichtiniturn sienifies tavreduction in blood 
ProweGgcotitthematunge Aestudyvontradiotherapyepatients has 
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weeks after the start of radiotherapy. The limitation of 
theomtestiis “that (in patients with lunge tumor, if cannot 
be distinguished whether the changes in blood flow are due 
to radiation damage or progression or regression of the 
tumor itself. To overcome this, a radiograph may be taken 
toglocaltze thepetumorerand, excludes .thisgeregiongaqduring 
analysis of blood flow. 

Another nonsinvasive (method) which «can, slocate, the 
Gumonr leas i nwells,s,asadetect lung density «changes iin=viveo is 
X-ray computed tomography (CT). This method has been used 
by@severalyvauthors in, the, study ‘of human. lung-and found “to 
be superior to planar chest X-rays. Encouraged by the 
early increases in mass observed by Sharplin and Franko 
€1,02>', thepipresentelwonks wass sinitiated= ~ingrordex to 
determine whether the increase in mass observed early 
aftersitradiationgmightsalsopbe detected asvan increase-in 
denshtyeiby kusings suche+Cihw scans. This § besteawtiliebe 
described finidetail intGhapter,2. 

In abl tests vethe™ <aiime (fsieeto “quantify p,earckyeyt lung 
damage non-invasively, and to establish whether. such 
damage is predictive of later treatment complications 
C1 Je) 2 It is assumed that an early signal of severe lung 
damage exists during the “latent” period. Then procedures 
andagdrugssedesienedanton modifyerthiss response. «might be 
activelyesoughtasoythattthesclinicalasymptoms oferadiation 
pneumonitis may be relieved or avoided altogether. The 


non-invasive methods would then monitor the effectiveness 
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of any strategy used to reverse or suppress radiation 


pneumonitis. 


1.4 Treatment of Radiation Pneumonitis 

Since an inflammatory neaction is part» of» radiation 
pneumonitis, corticosteroids are expected to suppress this 
Reaceton,, because these drugs are Aantiainat lanmatony, 
suppressing connective tissue reactions and decreasing 
cellular vtinititrativons. sindeed;. prednisone “administration 
prior Ato scirmadiationwand untilvG0 tdays aver was found «to 
reduce smortality “due to radiation “pneumonitis (CASE 
Cortisone has been found to reduce the severity of the 
radiation response in lung assessed by measuring changes 
iio cOralesechest) compliance: inmmrats: (7.8). Cpe sterodds are 
withdrawn in patients already on steroid therapy, the 
symptoms OL radiation pneumonitis are aggravated. 
Antibiotics have no effect on the radiation response but 
anemabenehicial nti) minimiizine *sopportunistic *fintections 
which are often concurrent with pneumonitis. 
Radioprotective or radiosensitizing agents which render 
ine Stiimons: Nore csensitivemto, radiationyvthan™” mormal »* lung 
tissue can be considered to reduce incidence of radiation 
pneumonitis. Fractionated radiotherapy instead of single 
doses causes more damage to tumor relative to normal lung 
telisisate td uieirto fedifierent: tmepain eicapaci ties) » (3:21, 84))i< A 
preliminary ‘study using a Chinese herb medicine named 


HAeoOvoritommodifysradiation’ response off lungiwase reported 


| 7 
eeise2eaar 30 eq, of Ambauend! sips Wms Lon) + aan a 7 

siagtze pu sbi éceseonlives ;bit beomuang, a 

penotaeedy )t 2244 i sarOdQwh AaB" « er copes. sealed 


ocltaes? heel, ee yaRres aetenotieg ie, 


Seschaees Ohta vebyiies 


S54 507 10i n)wre grips ine ss , ue a hd 


svat. Of Lope Sek euljeithzal oF volsq 


Ss 


bles 606".94! J 


Qe dan +7 Sah be? 


aheheicte sav 1 a59% 

cae) sic inuwy Dy uargloethne O28 out giitesaee soothes 
ers Yo Yv13 jf Soyhsd, os (beist ape) eer avast seeo 
cusadc. sities) 88 Seadeeee Bans. Gb vaneqtea oncsesbay 
me ee ee 4°) «dev Bicecapdlqaes , sears lerad eh 


1° , 46 7th ied 9 OoOSe ) vig tts ogns tka na eer beh he. 

: : + 7 
.oe3) Goth aed AS at¢areorushg no tag ies i oresq¢eED 
>is! enw lates nupse fbew Bile ae a teal wu oven wnteerd aed: pea) 


Seeitawin) ‘010.1 fuszodam npdeda lite ei cata tewed: 96 | 
Sibpswosaing . 940) Seeaneatieg — URES “MerRS # _ 
Sctces Arye pera pagent: aér evans 

et ere eee 
Bélietor? iv soanbiens - bei” sian a 


7 ni i 

) inom 

Danial te ees ’ ad iq ae air aa 
anpt 3 ee pe a : ia papain a 


& nd Toes. « sheer at a 
— ao Se 


dy goctl teetos . oad notre 


>. Sh ae 


a. 


by Shen and Liu (68,103). Radiation changes were observed 
in rat lungs by microscopy and electron microscopy. Lung 
damage was observed to be much less severe at all times 
amtenainradtatton leinntthe)@treateds: ratis:< The proposed 
mechanism iiof Ggaction “ofipthis,s medicinesaibs: prevention tof 
blood vessel spasms, Ano ante —intlammatoryes sact1 on and 
inhi bieion of collagen production. None of these 
treatments consider wall aspectssofs« nadtation “pneumonitis 
and therefore radiation pneumonitis and its treatment is 
stillvan area of .cactive, research. ihis® (has “vprovided 
sufficient incentive to develop early detection techniques 
which will permit more timely intervention, and assay the 


efficacy of treatment methods under investigation. 


1.5 Guide to Thesis 

fnethis thesis Gradtation pneumonitisuwis dealtewithoin 
two ways. In Chapter 2 a non-invasive quantitative test 
based on X-ray computed tomography for detecting lung 
damage early after irradiation is investigated. The 
reaction of the lung? to radiation is. monitored in  egroups 
of animals which have been irradiated to various single 
doses to the thorax.It is determined how well this test 
can detect radiation damage, its time course, and dose 
response. ) FOr this part of the” work, large numbers of 
animals were required so that mice were used. inorder. to 
compare this Des CaO OO Crm LV Vio bes CSc such as 


conventional X-radiography and nuclear medicine studies, 
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however, larger animals (i.e. dogs) were required. To 
further evaluate the technique in a clinical setting, the 
densitometry test has been applied ale) radiotherapy 
patients at different times after therapy. How well CT 
monitors lung density changes in these cases is reported 
inyesiection 2.4. 

in Chapter 55, radiationt dosimetry in’ larger halit-pody 
fields is dealt with and various methods to calculate 
accurately the dose to lung per se are intercompared. The 
various lung dose calculation methods are compared to 
measured doses in simple geometrical phantoms in order to 
study their range of applicability. These methods have 
also been tested in a humanoid phantom. A detailed dose 
distribution “ins ay transverse slice of the phantom is 
calculated and compared to measured dose values, and to 
the “manual” planning methods presently used for hemi-body 
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2. DETECTION OF LUNG DAMAGE BY COMPUTED DENSITOMETRY 


2-1 X-ray Computed Tomography (CT) 


Zalsl Principles 

A-fray computed tomosraphy (CT) is > (ae body imaging 
technigue that’ permits viewings “of “a  efross—section or 
“slice” without interference from adjacent layers. X= ray 
transmission measurements made at several angles or 
projections through a plane of interest, may be used with 
an appropriate algorithm to synthesize the contents of the 
Mayer of interest. “This technique is sieniticantly better 
im terms of tissue contrast resolution in comparison with 
radiography and conventional focal-plane tomography. 
However, it does not compete with film imaging techniques 
im’ terms of “spatial resolution Cat high contrast). Thus 
computed tomography is particularly suitable for imaging 
soft tissues (e.g. Lune) whale ti ime radtLoctraphy and 
tomography are more suitable for diagnosing bone 
fractures. 

Eney mathematics of construction) of “any “images of an 
Ob,eECcCey Lroml projection» “data -obtained at a number of 
angular views was first developed in 1917 by Radon (87). 
Other techniques were developed in 1956 by Bracewell 
WOLKA0g Mine radLoastronomy Glo). In 19/754. Gordon et. «al. 
cited more than 1500 references related to image 
CONS ErUCct Lon, (345) 


The first medical application of image synthesis from 
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radiological’ orojections#iwas publiished’*® by’ Mldendorf in 
129.644 <8 1s Other important experiments were carried out 
ine L963" by Kuhl vand colleagues utilizing a tomographic 
image seagmeraror | radionuctide detection “~“(61), and by 
Gormack using radiographic techniques (18,19).  Hounsfield 
pioneered the first clinically useful computed tomography 
Systema and it was’ installed in 1971 in the Atkinson 


Morley Hospital near London, England (2,47). 


2.1.2 1tansmission Data Acquisition 

The simplest CT system consists of an X-ray tube 
collimated to produce a pencil beam incident on a 
detector. The plane to be scanned is placed between X-ray 
tube and det@€ctor as shown) in Figure 2.1. The first step 
is to move the X-ray tube and detector assembly linearly 
across the object for a fixed angle. Readings are 
obtained at a large number of “points = during. «the Linear 
motion, and the flvence profile of the transmitted X-rays 
is detected to form a projection. When a linear scan is 
complete, the tube and detector rotate together through a 
small angle and the linear scan is repeated. In this same 
Manner, many sets of projections are measured. Li one 
defines a coordinate system X-Y through the object, then 
each epoint, “along "the tay is specified by its “coordinates 
Sy veand pa Lunegt ion aC x.y) is defined “at every point. For 
X-ray transmission measurements, this function is actually 
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integral of Ge seas) Abonce aeapatnd ds 5 21s “called <a 


ray-projection p: G16) 


p= fucoy ds 


For monoenergetic photons, the transmitted beam fluence 6 
t 


is given by: 


6 = a exp (-f uceydas) 


where: ) is the incident beam fluence. 
fe) 


X-RAY TUBE 


DETECTOR 


Big. 2.1 Simplest Cl system consisting of "X-ray tube 
and detector moving linearly across the 
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Therefore the ray-projection is measured as: 


p = -ln Dele) 


The ray-projections p from many directions can be stored 
in a computer and the problem is then to invert a large 
number of such measured ray-projections, p;, to solve for 
WGC, y. Actually the rays are of finite widths and the 
u's obtained are the average u in a volume element (voxel) 
Of finite volume. =slhesreconstruction process then, will 


create a map of y's throughout the irradiated plane. 


2.1.3 An Image Reconstruction Technique 

Suppose there is an object with a dense pin at its 
centre. Views from several angles will give the 
projection proficveswshown in  Eisures 2-2. An ficst “attempt 
at reconstruction would ) be to  ~baek=project "these 


profiles (i.e. evenly distribute the profile value along 


aera yi. summing back-projections (from all directions 
gives a star-like image as shown in Figure 2.3. helisn, 
however, is@eame Dliurred | image. Vof the. “pins the image 


Contrast 13s impaired and) the “true [natures of the real 


object “ls obscured this  tmage can be edeblurred “or 
deconvolved (by= —(pre-ftiitering) each profile pELOT to 
back~projection and summation. in our example, hileure 


owe ne promi Le WOuldmebe = ihltered so. that Lt had | both 


negative and positive lobes as shown in Figure 2.4. When 
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Pige. 2.2 Brogectione protiles#ofeantobject with a 


dense pin. 


Fig. 2.3 Image of a dense pin formed by simple back- 


projections. Note the, star =lLike artefact. 
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Pigee2.4 Protile arter modifitcation with, a filter 
function. 
the filtered back projection values are then summed, the 
negative values cancel positive values and the spokes of 
Paigure (2.5 disappear. 7 lt the ytilter function: is correct, 
the cancellation jwiiInkl* be ‘exact, and the reconstructed 
image will be a corrects image of the object (31). This 
procedure is now used on all commercially-available CT 
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2.4.47 GEe Cl/T. S800 Scanner 

The scanner used in our experiments is a "“™ third 
generation” scanner which uses’ an incident fan beam of 
La—rayspeauult. pe waetectorcmand: rotationad motion “only to 
acquire data. Elimination: “of the linear motion makes 


scanning possible in less than 10 seconds per slice. This 
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type @of Siscanmer Tits pudiusitrated’ ein’ Pigure. = 2.5. The 
detectors are arranged on an are which rotates with the 
k-Tay. Vituber P about itiherpatient 4) The patient™ isivon a-‘tabiwle 
which can be moved longitudinally through the opening (or 
"Lonuth \atunm thetscantry tom producer a sertess of slice 
images. 

The X-ray tube delivers a pulsed beam of X-rays 
collimated to a wide fan beam. The generator-tube 
combination permits variation in the current up to 600 mA 
pero * pulse, “but “the voltage is fixed “at 220° kVp. The 


number of pulses per rotation may be selected as either 


a een 


X-RAY TUBE 


OBJECT _7 COLLIMATOR 
COLLECTING ELECTRODE 


HIGH VOLTAGE » 2000 V COMPUTER 


AMPLIFIER 


Fig. 2.5 Third generation CT scanner using fan beam 
OLeXerays,, array of detectors*and rotational 
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285.) 04 99/60, "and the pulse duration its variable up to 3.3 
msec. 

The original (Hounsfield) CT scanner used one sodium 
fodide detector coupled to a photomultiplier tube. Other 
scintillator crystals such as cesium-iodide, 
bismuth-germinate, and cadmium tungstate have since been 
used to achieve better light conversion and response 
times. The main disadvantage of using crystals and 
phovomultipliers, however, is) sthat ~ ae ts ditticult ”§ to 
“pack"™ many of these closely enough (eg- 2 mm) to achieve 
high spatial resolution. More specifically, 
photomultipliers cannot be made smaller than 1 cm in 
diameter and are expensive (53,54). For this’ “reason, 
photodiodes are used instead of photomultipliers on some 
scanners. 

Anvalternativewappuoach 1ss to wsie arian. array of don 
chambers filled with pressurized Xenon. Each detector 
element is about 25 mm long with an aperture of 1.05x1.05 
mm2 . Each “on “ehamber “tconsists of a planar coltecting 
electrode between two tantalum or tungsten plates, £idled 
with Xenon at 27 atmosphere or more. Such a chamber can 
absorb more» thang 7/07 of the incident radiation energy 
G53',0 40)". 

The General Electric CI/1T 8800 “system “used in our 


experiments uses the components listed in Table 2.1. 
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Table .2 0k: 


scanner 


Tube type 
Detector type 


Detector elements 


Rotation 

Scan time (sec) 

Number of pulses/rotation 
current/pulse 

pulse duration 

voltage 

image matrix 


pixel size 


Components of the General Electric CT/T 8800 


Rotating anode 

Xenon ion chambers 
D2pr Of WHICH Ceeende detectors 
provide reference readings 

360° 

doe O Le oO 

208 0r 376 

up to 600 mA 

lelesOG cece OTs omens ec 


120 kVp (fixed) 


3203320 


2-1.5 Tissue Densitometry 


Rather thanespeciiivime, the 4s, a 


the CT number Nor 


u = 


related quantity, 


has been defined as: 


(1) 


the CT number in Hounsfield (H) units 


linear attenuation coefficient (emo) of 


imaged tissue 
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— 1] 
U = linear attenuation coefficient (cm ~*~) of water. 
W 


On this scale then, each Hounsfield unit is equivalent to 
0.1% of the attenuation coefficient of water. Thus, the 
CT number scale is fine enough to distinguish fractional 
percent changes in the attenuation coefficient of tissue. 


Air with negligible attenuation has oh = -1000, and dense 


Y 1 


t+ - 
witnin 


"Sott™ tissues 1: 


ay 
m 


bone approaches +1000. Most 
-100 to +100 Hounsfields, except for air-filled lung which 
has a typical CT number value of -/750. 

The two-dimensional matrix of CT numbers may be 
displayed as a grey-level or colour image on a video 


" 


screen where each picture element (or “pixel”) is assigned 
a shade of grey or colour, depending on its CT number 
value. The range of CT numbers displayed on the grey 
scale is determined by setting a window, while the leve 
specifies the CT number value assigned to mid-grey tone. 

u as well as CT number depend on the type of tissue 
and on the incident photon energy. However, an advantage 
in the use of the CT number instead of u is that it is 


less dependent on the effective energy of the incident 


beam. This is because the ratio p/ varies much less 


al 


with change in photon energy than does y . In the 
diagnostic energy range used in computed tomography (60 to 
80 keV), there are three important mechanisms of photon 


interaction with tissue: Rayleigh scattering, 


photoelectric absorption, and Compton scattering (54,75). 
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The relative importance of each of these processes for a 
typical’ enersy of 65t keVe (ec ee X-ray tube operated at 120 
KVp J@facer shownteint Tabler 252. They’ totalieattenuation 
coefficient due to all three types of competing 


interactions is 


Table 2.2: Relative importance of Rayleigh scattering, 
photoelectric absorption, and Compton 


scattering’ tos total) attenuationvat 6545kevV. (95,75 


fyperoteainteraction, “contributd onto depends 


Potala eeenmulaltdomeate oo kel 


muscle bone 
GLERU,) 
Rayleigh 4% 6% 0 72 
e 
Photoelectric WE Dey 0. 74 
Compton 897% 67% 0 70 
e 
where: 
0 = electron density (e/cm?) 


0 = physical density (g/cm?) 


N = Avogadro's number 
Z = atomic number 

A = mass number 

Ph = (NZ/A)p 


E = photon energy 
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uta at [o, (EZ) + o° (E,Z ) ee) (2) 
where of and Z are the electronic Cross section "(cm "fed 
and pipe centre atonie mumber “for “each type “of Unteraction. 
then, susine equation, (ll) Sand 2) one” “can” "obtain ‘the 
electron density of tissue, relative to water, p. ase G99)s 


be =e On / 0 = CL OO) ae OOO IN 3) 


e e ew cr? 


where Ro ise the ravLo ote total eleceroni a sicrocss=section 
LOtewewater wlOustne mtotaie electromic "cross = SectuL1on LOL 
tissue. in orders ycomconvert —Cly numbers to relative 
electron densities, Ry Must be known for different tissue 
Compositions. Kor air, Lung, and muscle RG Shey oholoigens 7 Ce iS 
DiUsy. el Ome Somt “Cissives: with atomic tcompogmtion Similar 
to) that of water, the Cl numbers can be? “converted to 


relative: electron, densities according to: 


' — 
Dae 100074 +0 00% Nor (4) 


The above equation has been tested for various 
tissveesubstitute “materials, (82), and for tissues in-vivo 
(9) For water-“like;) tissues, it was found experimentally 
to “be accurate to better than 32 (ll). For tissues which 
differ from water in chemical composition (e.g. bone ) 
however, the influence of atomic number is seen. For such 
tissues, ae Must be pre-established with appropriate 


tissue= ‘substitute materials or measured in-vivo using 
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dual ‘energy CT *scanning %( 9). 
The *sravinetric density “of "tissue “1s “obtained * irom 


the electron density using the following relation: 


eee 
(5) 


(NZ/A) | + sue 


If lung tissue is considered to be a mixture of air and 
water-like tissue, where the fraction of air by weight is 
negligible, then: 

(NZ/A) 


water See 


(NZ/A) + sue 


ana“ tromeequatt-omss(4): "and: (5), Lt follows  sahac 


p = 1.000 + 0.001 N.. (6) 


Because CT numbers can be converted to tissue 
dengireies = ii “Pung, Sone, eas) for, Sethe  Preirse te “time = the 
possibility of "obtaining quantitative values “for lung 
density in-vivo. However, the lung is an inhomogeneous 
organ, with local density varying between that of air to 
that, “ofS soft) tlsisue*and blood: the composition includes 
Dlood*vessels; "air spaces “of the alveoli and” bronchi ; and 
lung parenchyma. Therefore one is faced with the problem 
or choosing ‘which *Vune "region to “sample “in order to obtain 
a representative average density. In@addttLoneto spatial 


density gradients throughout the lung, density changes 
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OcCuresgwith timesduring) thes breathingewand cardtacs cycles; 
such variations have been measured by CT for LUSDl rat 1 On, 
expiration, and continuous’ breathing by Van Dyk et al. 
(129) <.elneour.experiments;yseveral methods,» lofusobtaining 
an average density for lung were investigated, including 
sampling of .smali»areas.or linear»\ profiles, sof the, “lung 
image. However, outlining the complete lung area, and 
averaging density therein, provided the most reproducible 
(although not the most sensitive) method. This method can 
also be extended to yield a measure of the proportion of 
overall. lunge areas whiehJehas sangalteredadensity. «This 
Subject lis discussedyturther-in® section, 2.1. 6. 

Lucene tnOdewOfesatantitatives eC. densitometry is 
applicable to diseases affecting the entire lung such as 
interstitial disease which includes pneumonitis and 
fibrosis, and diseases affecting the general fluid or air 
eoutent fof thealune) 16136). The method would be less 
suited to measure densities of very localized lesions 
because.it would be difficult to reproduce the ‘sean’ site 
and to avoid partial volume effects. 

In summary then, CT densitometry of lung can be used 
to diagnose lung abnormalities such as _ interstitial 
disease, geradiativon” /injury, @¢pneumonitls » sand ssfibrosi's ); 
edema, and emphysema. CUWeca na gals 0 Pibe Qusied {to cdrecord 
changes in«phystological function of lung since diseases 
affiecr pulmonary perfusion and/or” ventilation . CT 


densitometry has already been used in animal and human 
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studies as a non-invasive sensitive method for detecting 
small changes in lung density (see Table 2.3). 

Recently CT has also become a useful tool for 
planning radiotherapy. As well as providing detailed 
three-dimensional internal anatomical data fox the 
individual patient, it gives the electron densities of the 
inner structures which are necessary for accurate dose 


ealcudla tions: Gull 2.9). 


Zot. 6 Precision 

Forw «each, scan “session, pins of known material 
surrounded by water were placed in the scan field. The CT 
numbers of each pin as well as of the surrounding water 
were measured, this, providing a check on long-term 
reproducibility and scanner stability. Results are shown 
inte Beloit Cpe 2.5 OF Each point represents the average CT 
number for a particular pin for about 30 slices taken on 
the date indicated. The number of pixels per sampled 
rection waseu20.) The "slice dose eperescan iwas, 3.975 gt 9 0 31 
cGye The standard deviation in the average CT number for 
VaArLows matenlals .iorn 30 slices was found to Demet. 
Which. “COnmesmonds, Sto same uncertainty, of 0.6% forsa lunge 
densitye Ot .04e/.cm oe lhis sig phe sinherent precision or 
“noise” of the scanner when used for densitometry of 
inanimate specimens. 

Dive Tease Sens OS Sammi Hep he C lgn On associated with 


reproducibly definins the sample, region of interest. For 
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Tablerz.4. Densitometry of -Lung 


Species Lung status Remarks Reference 

human emphysema and normal Rosenblum (92) 
human normal Rosenblum (93) 
human SnoOr mann Vane Dyke 7 -GL2.9%) 
human abnormal lung Wegener (136) 
human normal and abnormal xenon inhalation Herbert (44) 


pulmonary ventilation 
human abnormal Compton densitometry Dohring (sia) 


computed tomography 


human normal changing air volume Robinson C91) 
human radiation damage Nabawi C79) 
mouse radiation damage Pela ia ape Ole) 


Vane Dyes ely.) 
Labpit radiation damage compare Ci to X-ray Hermann (45) 


and scintigraphy 


dog edema oleicVacid ina, y Hedlund CaOs) 
dog and changes in density due to changes in Hedlund (42) 
baboon blood, extravascular fluid,” or eas feontent 

dog changes in density due to Hedlund (a) 


hemorrhage and rescuscitation 
primate regional ventilation Xenon inhalation Gur C390) 


and sheep 
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an objective trace of the lung outline, it was necessary 
to standardize the viewing “window" and “level” settings 
(60). Using simple test objects, the best correlation was 
found between the traced area and the true area for a 
level setting half-way between the CT value of the 
Structure of imtereste and, -that. of ~ethne wesurrounding 
material. The window setting did not affect the area as 
significantly. Several tests were then performed in order 
to determine the most reproducible method of obtaining the 


average CT number (N..) of the sample region of lung 


GT 
Ina Vivo C2, 129 13.6) ...5 Lhe followings meth oder: illustrated 
invFigure 2../,) were .considerned: 

Ca) adjust the level as™ above and traces close to the 
visible boundary—of lung (2 users), 

Cb)? places a> rectangular’ cecion (of interest, (ROM) box 
inside the lung and obtain the average CT number therein, 
Copp place, a ROL inside Wunc, obtain the averace CI Pnumber 
therein Nes a Oo ) where o is the standard deviation. Then 
find all pixels having CT number values within Nop + 2 One 
and recompute the average. 

(d) Print, outsthe Cl numbers for the) Lung. region and 
calculate Nor along a linear profile through Lung: 

These four methods were studied on a normal and an 
abnormal mouse lung. Results are shown in Table 2.4. The 
O reported there is the standard deviation of the mean for 


numerous ir Las of each method and is thereby an 


indication of ther reproducibility of, each tracing method 
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2.7 Various sampling methods used to determine 


the average density in lung. 


NOILISOd 


32 


oy 


A[TrTaeayrqae ussoyo 


a[tyzyoad jo wuorqyrsod 


Ssun~T ~[Tewzouqe uf 


SSAn}jONAAS Opn [oxXs Aew 


¥SOAaq,UT FO uoLTsaA Jo 


uotytsod uo yuapusdsap 004 


AATTLqQronpoidsda yseg 


2 


Bun], eAtTyus ayzduwes 


sjzueumwog 


C15 Co OL 


GL OURS 


T7209 9c Ser — 


OG OG a7 O1Gis 


0 Ly 


sun, [Lewiouqy 


Ste 9 © 


93276 


Lament 


Ge Tae 


St CGS = 


SOURS ee 


Clee leGi 


Lo 


SuUN[T [eW1ION 


$3Ss0}3 3urrdwes snoriea Aq pautTwieszyap se Aequnu jy aeseraAy 


CP) 


C2) 


(q) 


Ce) 


pour en 


<2 Co ee 


= i, _ a a 7 
eee? qoeilques ausiver vd bettevisnb-2e vadmaun T2923 
2 4 : 


= ras 
- = 
” 

~~ 


= — 


= 7 = 


per se. Eeom- this Tableyerthistapparent thattmethods (a) 
is the most reproducible and objective for normal as well 
as abnormal lung. 

Toethes@dogs, ~the reproducibility tofhthe@scan sites as 
well@ias@etheP tracing. procedure) was *veritiedt by repeat 
scanning on the same dog on the same day. There were 
eight len silicesa -Ehrowgne thet thoraxg Lirome fapexg to 
diaphragm, the sternal notch being designated as origin. 


The average CT number N was obtained for right and left 


GT 
itunes ineeach slices@ifor bothe scansHacivyaluesoriors rcrisnt 
and left lung as well as area traced were compared for the 
eight nominal positions. The maximum difference in N of 


CT 


a given slice was 27 H, but most differences were within 
10 H. The average CT number Mo, formats @alices hore, the 
dopEewase-JOSoe olen eS eine thew nicht lumens amde—7 fi). 72245). U1 
im the lett lunes "On the repeat “seam the average CT 
numbers were within 10 H of these values. The difference 
in average CT number for the entire dog lung therefore was 
108 Hounsfields;, “which represénts 5% in®@lune density. If 
one slice through a highly vascular area of dog lung is 
traced steneetimeseeaccordins™ tol method (Cagne fabiee2 4, ha 
standard deviation of 4.3 H in the mean CT number for. the 
slice’ was obtained, indicating the reproducibility of the 
tracing method alone. 

imethepvhunan, vamrepeat trace ofethery sameteslice "ten 


times according to method (a) gave a standard deviation of 


U.4ehHs in’ thetmean CT, number=for= the slices hisses lower 
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thani@g@for Bthe edo; Bdticletosthe tiactrthatethe «CT stice Loft 
human lung is more homogeneous and the boundary is better 
defined. Therefore we conclude that the reproducibility 
of method (a) is very good and this is the method adopted 
throughout ethisi ework. however; athe’ sensitivityecot this 
methodrinedetectings localized  =changesuis  fnoteavery  e2o0d ; 
simcesgtheses vanes diluted@eby lathes Laveragingseprocedure. 
Vitsualaciginterpretationgeans besmorevappropriatestor Tthis 
case because pattern recognition is exercised by an 
experienced radiologist. CT densitometry is more suitable 


tosdetect: subtler difiuserchaneves throughout the lune. 


Zelient ekeccuracy 
To obtain densities from CT numbers the _ following 


steps discussed above are used: 


R (NZ/A) 
Oo water 
a TS eR 
Hy, Nor Degas (NZ/A) .. p 
tissue 


Thegerrogs sinvolvedgin obtaining ssCI) numbers Proms diay i 
will be due to beam hardening effects, system noise, and 
Motlon=-iuadueeds, artefacts '<«(64;,24,753). These must be 
checked for each system and may vary depending on size of 
object; imotion;j and’scanning technique Ceg. kVp). The 
conversion to relative electron density is performed using 
equation G3) a) Crayvimetriiestidénsttties, sare) obtained from 


tewative ‘electronadensitiessusingtequation:(3.) efwseetion 


(NZ/A) 
2120 oe Valuesgo ft oR. ; Een ret eoletet in tebativer telect hou 


(NZ/A) = osue 
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density, andeeeravimetric density for water, air, muscle, 
and some’ tissue-substitute plastics at an effective CT 


enersy «of 65. keV vVare given in Table 2.5. Often R_ and 
fo) 
(NZ/A) 
water 


(NZ/A), . 
tissue 


are obtained simply from equation (6) of section 2.1.5. 


are assumed to be unity and densities 


The errors introduced by-= this. simpler procedure | for 
various materials are shown in Table 2.6. The densities 
calculated in this way are all underpredicted, relative to 


the true density. However, for soft tissue (e.g. muscle 


and lung) tthe sainderprediction is very smali (0.72). ee 
greater 2CCULTACYs a Sam Ceqim red, this, Pofiset™ could be 
conrected ~ for: ine tics) iawork,, fractional changes in 


density, Ao / 0 , were measured and this correction was not 
applied. 

In order to test the above theory, materials of known 
densities were scanned (see Table 2./). Results are shown 
plotted in Figure 2.8. The cork and polystyrene values 
were obtained from several slabs of variable thickness 
(e.g. Oe em) and Waitn small inherent density 
differences. These slabs were piled and scanned. The 
cork slabs were variable in consistency and density, as 
well as having uneven edges. This made comparison of the 
measured density (sample volume 200 cm?) to the CT 
calculated density (sample volume 5 cm?) difficult. The 
composition of cork and wood is unknown, therefore AG and 
(NZ/A) cannot be determined but setting them equal 


tissue 


to l as in equation (6) gave acceptable agreement. iss 
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is reasonable since wood and cork are probably water-like 
biological materials. There was also an uncertainty in 
the physical density measurements (+/7%) and variations in 
density from slab to slab. The polystyrene used was_- the 
opaque polystyrene containing a TiO, filler whose 
COncentration may svarysup to 7 44 by welsght le(137)). The 
exact amount is unknown. ime lable. 2a -demeatvevaluwese for 
polystyrene are calculated assuming filler concentration 
of 1% and 42%. 

In order to study the influence of sample size on 


accuracy of CT number, a test was performed in which three 


Ta Dive 26 0.2 Conversion of CT numbers tio denstties 


overall error, 
if equation (6) used 
—_—_—_—_—_—_—_—___———— i 
Hey Ee ome DO! Ce 0 instead of 


equatrvons (3) and! ©) 


ACE YE +17 -11% -10% 
Muscle +0.32% —12% -0.72 
Lung +0.3% he =O 0% 


Polystyrene 


clear -6% 1/5 = 9% 
Poly et lis TiO, —47 aie 187, 
Poly.+42Z TiO, sey a eaA -47, -3% 
Lucite Shy =O -6% 
Nylon -47, =. % —~5 2.54 


Polyethylene =6%, api moh 5 
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Sualliecorkuecylznders ofediameters 6 mm, 1.25 cm, and “1.9 
cm were inserted in plexiglas sleeves of outer diameter 
2.54 cm. These were then inserted in a water phantom and 
scanned. Resultse indicatedtvariations of less than 12. H 
forrsampled regions tanecings from OF 045scm? toeQ0ee75 cm. 
Asdagfinatetestge4tcylinders of matéeérialsdof “similar 
volume as mouse lungs were substituted into the mouse 
thorax. The mice were placed in the phantom and scanned 
in the usual manner. Results are included in Figure 2.8. 
Thisee indicates, that the eald bration linearity of 
CT-derived density is maintained whether the samples are 
Within» the mouse: or within. calibration, phantoms. The 
errors), involved, in  obtadninge \lLung’ densities - from, CT 


numbers are summarized in Table 2.8. 


Table 2.8: Summary of errors involved in obtaining 


lung densities from CT numbers 


Source lie ie Oe 


Scanner instability +062 
Reproducibility of the 

tracings procedure for: 

mouse lung +7050037 
dog lung ee LO Ona 
human lung +e 0.0074 


Accuracy of the conversion 


of CT numbers to densities ca Ne A 
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2.2 Experiments with Mice 

This study was initiated in order to investigate 
whether the increase in mass observed early after 
irradiation, (102) micht. also be detected in-vivo as an 
increase in density by CT scans.. The mass increase of 18727 
observed at 6 weeks post-itrradiation to 7 Gy= by Shtarplin 
et al. (102) would be easily detected by CT, assuming 
lung volume does not change. The work described in this 
section (2.2) has been published in the International 
Journal of Radiation Oncology Biology and (Physics (C25). 

Female mice: dfethe: Dolby Cee Chumstrainme were obtained 
Erom the® SmallgeaAnimale, .Propramie of. genes University of 
Alberta. Mice weighing between 19 and 23g were irradiated 
at 2 weeks: “ot, —iace:. invonder. to minimize the risk, of 
infection, Omtetracycline  at® a vfcoteentration of “425 
mg/liter was added to the drinking water. Tetracycline 
eontrols most bacterial infections in Laboratory mice; but 
does not control mycoplasm - a common virus infection in 
pneumonitis. Autopsies) of sample animals, however, 
reveaded” no wactecial)»ior mycoplasm) Hinmiections. Thus’, 
density increases could be attributed solely to the effect 


Oe elev GLieieel(lSL ele Ik @ial c 


DteeleMe ChOaS 
ine miceswere irradiated with an X-ray unit (Picker 
Vatgiand) SuSime tie parameters given im Ficeure 2.9. The 


nicer were, restrained at’ thes necks in: a plastic jig 
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3mm lead 


Pio weseo Micewun the irradiation jig.) Lie thorax 
was irradiated, while the abdomen and 
head were shielded by 3 mm of lead. 
M=ray. technique: Tube voltage: —Z60. kVp; 
bubermcurrvent: 919 mA; Halt—value layer: 


[teen Cn teld  siees 20r20) ems). 55D! 50 cms 
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described elsewhere CLO 2s) e. The unanesthetized animals 
were placed securely in individual "“half-tunnels” Figure 
Pinos ALT parts | of s the mice “except the) (thorax were 


shielded by 3 mm of lead to reduce the body dose to 4.1% 
of the lunge dose. The mean dose rate was 0.73 Gy per 
Minuterwit i sasdose:unitonrmity of 0.04 for ehe 1. 20x20 cm’ 
Paces Sie. dies This was measured using an air ionization 
chamber in conjunction “with. an » electrometer (Capintec 
Model 192). 

The phantom which is recommended by the American 
association: (of. Phycdictisots™ ain @ Medicines tor, Cl scanner 
performance evaluation (55) was modified to hold the mice 
during the imaging session (see Figure 2.10). An acrylic 
Resin (iucite) ender platew wlth? nine setubesem™ 2.6 -'cm in 
diameter, amd Beachy .s)5 cm lomnes was» consittructed. In 
addition, achyilic resinercylinders of “smalier diameters 
were constructed; the mice slipped easily into these 
inner cylinders, which were then inserted into the outer 
"sleeve" tubes. The mice were reasonably well immobilized 
without anesthetic and without having any part of their 
body tied or compressed. 

Five slices, each 1.5mm in thickness, and spaced 2mm 
apart spanned the complete mouse lung longitudinally from 
apex to diaphragm. A slice through the phantom with all 
Hicesw 4 aplace is shown ain BPigure 2.11. The appearance of 
ChemGlgslicesmate Various: posduLons in. the lungeis shown in 


Figure. (2.12% The data presented here were obtained from 
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Pic.eceteu Appearance of the ‘Cl “slice at different 
longitudinal positions through mouse lung 
CA) TVearechemapex. (3) and (CC) mid—lung. 


(D) near the diaphragm. 
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thewmid=—lung slices, "the Location of swhich was judged “in 
terms» of largest cross-sectional area of lung im the CT 
image. The mice were scanned with the CT scanner 
described in section 2.1.4. Starting at 6 weeks after 
irradiation the mice were scanned once every two weeks 
until 14 weeks, after which they were scanned once per 
week. For mouse lung, the viewing “level” was set at a 
value Ot 2 UR HOUnsttLeldcm and thes window Jat 9100 
Howms i teldg Seomracmhi tates turacine= Ofe thes lung “outline 
(60). Under these conditions, the lungs appeared black 
and the surrounding tissue appeared white. An area of 
interest was traced as close as possible to the boundary 
Opec iemLlung. The only possible disadvantage of this 
tracines spechnique, is that radiation, -efiects near, the 
pleura may be excluded, and thus changes in lung density 
were probably underestimated. 

In mid-lung the boundary between the right and left 
lung was not always clearly visible. However, since both 
lungs were irradiated and showed similar effects, the 
derinmveLron, of (the ara ont-rette boundary wassnot critical. 
Examples of an unirradiated normal, as well as an 
imradiated, -mouse™ Lung’ are, shown. in Figure 2.13 for 
different window settings. Note that the density 
increases observed during® radiation “pneumonitis are 
localized as dense patches. 

The radiation dose received during one scanning 


procedime: was) o0.0S9+0.001 Gy (112). Over the whole 
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scanning regimenyot 25orweeks o(16 sprocedures), the mice 
WOoUldImthererorewereceiverdantsadditionaludose oiseeo2 Cy. 
This dose is certainly significant in magnitude but since 
the lungs have capacity to repair sublethal damage, a dose 
Obe OS OZ nGy line loOvfractions ovyerrZ5eweeks sisti@expecteds to 
produce nosnadditionals effect C132) 2enhHowevery lthermice 
usedoin, our studies:did exhibit much less’ tolerance to: the 
radiation than was expected. The LDG50/160)° Band 
LD¢50/180) forfiune irradiation in. different strains of 
mice have been reported to lie between 13.5 to 16.1 Gy 
CNOGR132)— LEor Balb/e micewased iniother Texperiments “at 
Oure rinstitutes (lO2), PthesLD(50/160).wass determined to be 
13 Gy, in agreement with the above data. These mice were 
irradiated (but snot oscanned )»Jandtkept) undececondPetons 
similar to those used in the experiments reported here. 
lt.) ts) speculated ithat the additional dose of 0.62 Gy and 
the stress of weekly handling exacerbated the lung damage 


or interfered with recovery. 


twliec Results 

A total of 142 mice were used for these experiments 
with the dose received reported in Table 2.9. 

The average CT number for an individual animal, Nak 
was measured in the right and left lung. The average CT 
number for an animal group irradiated to the same _ dose, 


was also calculated and is designated Nog » torn Right cand 


left lungs separately. 
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The Cojsnimbenssitor control lungs of the individual 
animals Saray were found to lie between -475 and -550, 
corresponding to a normal range of lung densities of 0.45 
COn0n 2 ecm « eA pl Ob. OF Nop for individual animals’ i's 
shown in Figure 2.14. The average CTI numbier for this 
control “group UN) was -S25 with a standard, deviation of 
30 Hounsfields tor the group. No significant, deviation 
from these values was found at any time for mice 
ipradiated to doses of (Sifand a. Cy. AS pl Ofte aot? wet nicws: GL 


number for, -abl we control, fanimals (N as well as those 


cpt)» 
iyradiated to, 7) Gy are ‘shown Gin Figure 92. U5. For the 
Soniroly “groups! ‘and. “thes groups trradiated) to Sand: / Cy, 


the standard deviation from the mean was maintained at 30 


Houwns toenidisi. Thts demonstrates fhiaite thiew Cal nam! bens tor 


Table 2.9: Number of mice irradiated at 12 weeks of age 


to the various radiation doses 


Dose (Gy) No. of Animals Surviving sat 
25 weeks 
0 7 23 2 
5 18 18 
yf 28 28 
10 28 5 
12 tage: 0 
ts) 15 0 
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animal lungs were reproducible to within 6% Cis er 
307. S000 This Vata ton is Stcributed omaandy to 
biological differences within the animal groups since the 
feproducipility of “the CT scanner data alone is50.6% for 


lune-sSsubstictute. materials, as noted’ in section 2.1.0. 


Lntiot. lL Aciwtes Phase 

Mice receiving doses of more than 10 Gy showed a 
gradual increase in density beginning at 15 weeks 
Following dtrradiation, with targe varlations in) todividual 
animals, as shown in Figure 2.16. 

The time course of density changes for mice receiving 
13 samdeoa4 Gy are also shown. (These curves are similar, 
with more severe changes observed for increasing doses 
(see Figures 2.19,2.20). The appearance of the scans as 
the mice go through the pneumonitis phase is shown in 
Pisure™= 2.0/7 forms individual mice. At Yoweeks no effect 
is visible or measured as a density increase. At 16 weeks 
no effect is visible but a significant CT number increase 
is noted. At 18 weeks CT number changes are both measured 
and) Sclearly visible. For the mouse irradiated to 7 Gy no 
change is visible or measurable at any time. The large 
individual variation 2s, Lllustrated in Figure, 2518 where 


N SS plotted fore. sindi vidual animals irradiated: [tor 910 


CT 
Gy. 


Student “s.) t-test was performed to ers the 


Statistical sipniticance of the density increase: at 16 
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OCH ren 10 GRAY eright lung 08 


‘left lung 


ATI OSS ce 
Cee CPP Oce ue eeo, a lomo. 22) 25 
TIME AFTER IRRADIATION (weeks) 


Average CT number (Non) COL wsCe I eroup of mice 

bi rcadiated StomlGyGy, .GB) growprot mice drradiated 
Domes Gy, (C) srounpotmemice irraddated, to 14. Cy. 
Corresponding densities are shown on the right 
ordinate. Typical error bars represent the 
standard deviation from the mean of the group. 
Error bars increase during the pneumonitis phase 
(16-22 weeks) because of greater variations in 


response of individual surviving animals. 
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weeks “after darradiation to 10 Gy. The level os 
Sicniticance "was ) found, €o0) be larger than 995997. The 
changes observed at later times and at the larger doses 
were as significant. 

Deaths occurred between 19 and 24 weeks. Mice whose 
lung density in both lungs increased by more than 40% from 
apnormate walue Sof) 0. 5s 2, cme tor moehO0sl7 s picmmewoald 
generally die within the ensuing 2 weeks of such change. 
The above data are summarized in a different manner in 
Eieures 2.309 sandy2.20.,° lm Picures 209 * the time course of 
the changes in density is shown. The percent incidence of 
animals with a threshold increase of greater than 404 in 
fone = density ts plotted as ta Eunmctifon of time. Figure 
Z.19° shows that major changes! occur in the time perdod of 
16 to 21 weeks with more severe effects occuring slightly 
earlier at the higher doses, in good agreement with the 
data of Travis ‘etueal. CL i6e wel 8*)i. 

Figure 2.20 presents the same data plotted instead as 


a, function Sot radiation, dosemiorm the racutes “phacsewot 15 
to 21 weeks. This plot suggests that the severity of the 
effect increases with time as well as with dose, again in 
agreement with the findings of Travis et al. lV Oey). 
However, the 40% increase in density chosen is an 
arpucratryeenceshold. Thus, ©£teures 2.19 and” 2.20 emiteht be 
different if other threshold values are chosen. 


Two of the mice were autopsied and the 


hdstopathological findings in. lung corresponded to the 
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2.19 The percent of animals whose lung density 
increased by more than 402 is plotted as a 
[UnctL1on Of timer artersitrradilation ton the 
PLuUUDSsODeiCe. trraddactced to 1:0. 82 eal. 
and 14 Gy. The incidence is defined with 
reference to the initial number of animals 


(100%) in each dose group. 
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2.20 The percent of animals whose lung density 
increased by more than 40% is plotted as a 
funetion of dose during the acute (radiation 
pneumonitis) phase 15-21 weeks. The number 


of animals with severe effects increases 


with dose. 


classical description of radiation pneumonitis. 


Histopathology report as obtained from Digs L.D. 
Armstrong, Animal Disease Section, Laboratory Services, 
University of Alberta. 


“Mouse No.l - The lung tissue examined was much _ more 
cetbulari®@thansenormal Juwithrtlargse (numbers tottalveolar 
macrophages present in the alveoli. These macrophages 
were swollen and vacuolated with some golden brown 
staining pigment. Interstitial areas or alveolar walls 
themselves were thickened due to cellular infiltration 
which was primarily macrophage and mononuclear in 
nature. Some interstitial areas showed evidence of 
fibrin accumulation and this is what was judged to be 
in some of the macrophages, giving them a ground glass 
appearance. Occasional double nucleated macrophages 


were encountered. Very slight mononuclear 
accumulations were noted adjacent to some of the 
airways. Airway epithelium itself was moderately 


uniform with moderate anisocytosis of nuclei and with 
SOMea er jumbLline setup) tobecellis into multiptes@layers.. “The 
lumenal border of these epithelium cells was jagged 
with ~individualexvcetis eprotruding &Tandehavings rounded 
borders (Clara cells). A few macrophages were evident 
in the lument-of airways. 


Mouse No.2 - Was similar although more pronounced in 
these changes, large numbers of macrophages ranging in 
Gaze s el rom ynormaletoseisantic Twerne jahound Gwithin Cthe 
alveoli. These macrophages contained dark golden brown 
pigment or were very large and pale staining with 
apparent protenaceous edema fluid. Some edema was 
Haund Fim this panticulan lung’ talons Awithe mere atocal 
distribution —~of nypercellularity to the alveolar walls 
with the predominant cell type being macrophages’ and 


mononuclears. Similar changes were found in the 
epithelium of the airways, being slightly hyperplastic, 
uniform in outline and numerous epithelioid or 
maégrophagevetype acells iin» the Lumen. Mononuclear 


accumulations were more prominent adjacent to airways 
in this mouse.” 


The lungs of several mice were also examined post mortem 
by microscope. The lungs were severely edematous with few 
functidnalbpalvyeocli.s. uihhe=appeatancée wwas dark «red; weand - ino 


ous, Anormtally Jindicativesof infection, was observed. Thus 
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the changes are entirely attributed to radiation effects. 
Only one animal irradiated to 14 Gy survived the 
early acute) phase Tybut) (died at 309weeks. Data for this 
anecdotal case are shown in Figure 2.21. The CT numbers 
for the right and left lung are plotted separately. The 
density increases during the acute phase were confined to 
the apex and diaphragm CT slices for this mouse. There 


was no density increase in the three mid-lung slices. 


Therefore the volume, “of (lune Sattecteds bya caditation 
pneumonitis was much smaller for this mouse, which 
explains its longer survival. However, the density 


increased progressively in mid-lunge with “the Left lung 
showing more pronounced effects. At the same time the 
appearance of the lung took on a different shape as is 
shown in Fireure: 2.22; with the pleural boundary Dednoe? very 
irregular. nis case suggests re-organization of slung 
structure @tow’ cope with damage. However, despite this 


remodelling, the mouse died during the late acute phase. 


2.2.2.2 Intermediate Phase 

At 30 weeks after irradiation, there were no 
survivors among mice irradiated to doses of 12, 13, and 14 
Gyo However, at) a sublethal dose of 10 Gy, 5 of 28 (182) 
animals survived the acute phase. Between 22 and 39 weeks 
after irradiation, the average CT number for this group 
(N 


cr resolved to near-normal values, lying between -420 


tow -460ehounsftedds, (i.ces=mdensities in’ the range of 90.54 
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mL Clescanmor anmouse: lung), trradtacedm@ lo aia) Cy 
(CA) Zoe weers atiter irradiat.on, 
GB) Zoo) weeks arcter irradiation. 
The altered shape of the lung is apparent, 


partwewlabiy in the pleural region. 
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toro). 5 0 eaey cm). However, with so few animals surviving 
the acute phase, data on the intermediate and late phases 


are very limited. 


262e2.5 Latemenase 
Two mice irradiated to 10 Gy survived at 52 weeks 
after Lr radiation, ater woicnums tamer” aa -final scan was 


performed. The following values were obtained: 


Animal R-Lung Lune 

Nera o Ce cme) Nye O (eyica ) 
1 =19 2 0 9 Ol 08 = 3 13+ Seg mee Oe 
2 = 6 On) 2m plo =S 4270/3056 517. 


in comparison Por control density values apetchis time, 
Onto OD)  Oses47 9 =517 2.48 OnGi83 

Thus the density is elevated with an appearance as_~ shown 
in Figure #2725.) One Otethnemmucesonowed slight structural 
changes, again indicating lung remodelling. For both mice 
the: lett Gliunew wasewmore jatiected, than tie right. The 
reasons for Ciiemelight= lett asymmetry saves not. clear, 
although  tLtglismeknowno that structurally Seach side is 
different in the number and size of component lobes. 
These mice were subsequently sacrificed and preserved for 
analysis by immonufluorescence_ to quantify collagen 


deposition during late fibrosis (76). 
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Fig. 


Dwell scans OLetwoemtce Lmradiated to 


10 Gy 


Deo weermse alter ipradiation, (A) mouse. le 


(B) mouse 2. 


DEEUCtural changes are 


apparenty only: in (CA). 
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2.2.3 Comparison with Mass Measurements 


2e2.36l Raditation-Induced Changes 

No increases in density were found at 6 weeks after 
itradtation, the time at which an increase in lung mass 
had been reported (102). Assuming a constant lung volume, 
an increase in mass was expected to produce a 
corresponding increase in density. However, a comparison 
of mass increases with CI density increases for 9 mice 
Meradiated to ll. 5 Gy at Vly .5 (weeks) after irradiation 
inagtcated an VeryamepOoOLanCOlTre vation. haisivdd spa a ittye 1s 
probably due to differences in volume and density of 


excised lunes in "comparison with in-vivo” lung. 


2.2.3.2 Adrenaline-Induced Changes 

A further comparison was made by inducing changes in 
lung density “in a2, more controlled way, without usine 
radiation. An injection of adrenaline produces edema in 
the lungs within minutes, thereby increasing lung density 
(94). Mice which were injected, sacrificed and had their 
lungs weighed showed an increase in lung mass relative to 
controd, mice, and this effect lasted about 15 minutes. 
ALter this time, the edema began to clear slowly. In 
Order CO. -Optadan —ditferent mass increases, various 
concentrations of adrenaline were injected until the 
desired range of lung mass increases was achieved. The 
concentration of adrenaline, (Epinephrine Hydrochloride 


Solution, Parke-Davis lmg/lml) was varied between 0.07 
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mg/ml and 0.2 mg/ml. Thirty-six animals were used in this 
experiment. The dose, normalized to the animal body mass, 
varied between 0.69 mg/g to 1.58 mg/g. The animals were 
scanned before injection to obtain control density values. 
Then the animals received the adrenaline injection into 
the thigh muscle, and were scanned immediately afterwards. 
Iwoe tom three! CT@eslicesa#throughy lung were obtained! andy ct 
data were calculated in the usual manner. The animals 
were subsequently taken from the phantom and sacrificed 
immediately so that the total time from time of injection 
Lowe time? Sofe vdeath* wast less’ than) 20 minutes, i.e.) within 
the time that the adrenaline effect lasts. The excised 
lungs were weighed, and the mass increase determined 
according to the assay developed by Sharplin and Franko 
(102). This mass increase was then compared to the 
inevavo anereasel ing Civedensity for individual animals. 
The inéreasetin®ClT density* is relative to thei mean for the 
COMtrols® (0.47508 e/ cme)! and “the inerease’ in?) Gess® “Ws 
relative tor kas “baseViner lungs’ mass, calculated stor ia 
esitandardys mouse @wedphniner e211 “gmi (999. Comparing the 
percent increase in post-mortem mass to the percent 
imerease intin=vivo Cit’density and) analysinge® they data 
using*® Iimearr regression, a correlation cocfifitetent of 057 
was determined. Assuming a constant lung volume, a 
StraLent Adanelwith antintercept) of, 0.0nand a slope of 1.0 
was expected, but not observed. A 25% increase in mass 


OnLy nereasedae ther i dens itty — by? 625. ¢which “ean only be 
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detected gwith emodest (672%)° “confidence: “by CT. The 
measurement OF mass, although invasive, is a more 
sensitive indicator of lung damage since the entire lung 
is sampled. The CT tracing procedure avoids the pleural 


boundary of the lung, thereby excluding possible pleural 


changes. Thew density eincrease ‘which’ “oftens ,occurseas 
localized patches (see Figure Zee) is thus 
underestimated. Also it is speculated that the volume of 


lungaechanges with ethet-adrenaline “(and «possibly twith 
imradiation) ; and that a comparison of changes in mass to 
changes in density is not possible without simultaneous 
Measurement of volume. Furthermore, there are reasons to 
expect tihatisuchs changes will certainlyecbes dittetenmt: Vin 


excised as opposed to in-vivo lungs. 


2.2.3.3 Volume Measurements 

In inanimate objects of regular_ shape, ints is 
possibile to obtain accurate measurements of volume by CT 
scanning a series of contiguous slices (see section 2.1.6 
and¢ (1605) In=thei Live mouse, ithis isynuchemorepdithel cubt 
due to voluntary and respiratory movement of the mouse, 
and the, vankablewishapetotmthe lungs fromeslicegtogsiiice: 
Toieis tequabblyse@edititcult.Aato. obtain reliable volume 
measurements ex-situ due to problems with excision - 
possipléesleakagsetoGin bDhboodicor alr,” dnd ithe s@diririculty 
measuring the small volume using Archimedean methodology. 


Nonetheless an attempt was made to correlate density and 
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volume measurements by Cr with mass and volume 
measurements ex-situ. Measurements were obtained on 16 
micél itradiated » to —S Gy at.6.5 weeks after irradtation. 
The datavare 9 plorred— “in Eteure 2-24. The method of 
obtaining lung volume ex-situ was as follows. The excised 
lungs were tied at the bronchi to avoid the escape of air 
and were immersed ina vial with saline. This vial was 
weighed, lungs were removed, the vial was refilled with 
saline again and weighed. The difference between these 
two values gave the air volume according to Archimedes' 
principle. The lungs were then removed, bronchi cut away, 
and weighed. This weight divided by 1.05 (tissue density) 
gave the tissue volume. The ex-situ lung volume then was 
the sum of tissue volume and air volume. If both methods 
measure volume and density accurately the values would lie 
on the equal volume and density lines. However, volumes 
as measured by CT were lower and CT measurements 
consistently gave higher values for density. 

im order to inprove accuracy, larecer animals ssuch as 
rats (Sprague-Dawley) were used for mass, volume, and 
density measurements by CT scan and ex-situ measurements. 
The rat lung (2-12ml1) is more than ten times greater in 
volume than mouse lung ( 0.4 ml) as seen in Figure 2.25. 
Agreement of in-vivo and ex-situ measurement of lung, was 
Detter for ratsetnan, [or mice put Cl density Sand volume 
values were still generally higher. These differences are 


too large and variable to be explained by systematic 


70 


(ee: Sy ese sats sented ak 
Me intaew wih ake apegee et faszelq 
unilet en aev hha inde eae 
sy 24 ehn men oho @iave O8 tiienetd eft ou bats avew egaut - 
sar) «6-epedae. apie Opleoe (ot hwo rsannt exav) baa 
‘jona aay  dpekyoeds ,bevemet gaee epguat cbsdgtew 
eyed ssagreteah aft -Setgdew bee atege wnrtee 


Seolian-e walt 


“ae 
a) i ee ie @etbontoe . soweey T1e es pray esvibv 0¥2 | 
dy wit Lisqgeac ,soedeuc’ weds ease “gees afl «siqtsarsg - 
Le accei)) 20 i yd bebtval sigtev wide Lbedtgrew tne 
cut i. vl pajenagacst? = <@aclgs seurit? sft, > eveg a 
api?) « : v seule’ 22a, Yop “senlsy ewreal? to aoe edd 7 
at. om wouced atc) Cid ceesed Crteneh bas suetho <a 
oo) ae pie) Med ganti vatereh one eautov Saupe > 

Cieaniccanen SO bak Sabet. Sage 120 xe aptaak om a 


yt henehoted peplie tadgad weap, ha 
ee 2544 -Loeddia seQter Gs TI87 ees vargak as /a2ebt6) ; 
Loe ,comdon . year) 3e% heey ove . . 
Glerec wecen » giti-te bor ages SO: xe 
GP) ddiseay obese Sed DeereTe et 
2566, Vaetl a me fo ® 


pay .ace! ‘bY 


ia 


*dun[T esnow pojJeLTpeadiy AJOJZ sqyuowaanseow aunToA pue ssew nqTSs-xa 


O27 SjUusWsiNnseswW oWNTOA pue ARTSUSP (JD) OATA-uT Jo uostiaedwoyg 7°7 


(,.w26) ALISNId 1D (<4) JWNIOA LD 
SO TAY em Te 


TUES AG se eo ie 


ALISNIQ 90 IWMOA 90 


(42 6) ALISNIG G3LVINDIWD NLIS-X3 


SNM] ISNOW GILVIGVAY 


2 


(42) JWNIOA NLIS-X3 


¥ eyed 


EDMONTOKR 


- 
+ 
| 


HE 


biel 


i? 


(B) at the same level setting but with 


e 

m 

on) 

S 

(=) 

es 

Qy 

(S) 

jen fore) 

(=| 

(@) ro 

w v 

— Ww 

oD] 

Q4 WY) 
o) 

3 

00 e) 

c ns) 

4 ‘| 

om) a 

v 3 
fe) 

(e3) 4 

co) 

3 o0 

e) +a 

A) ise} 

c ml 
oct 

3 fou} 


Eo mel Gis candor eae orate Lounom cA eewith) a2  Lewel setting of —250 Handa 


FT] 


é 


spane” @SCey's 


j poe? ‘ 


diiferencesstiin @datageacdquisitioneser + analysissebyae both 
methods. Lung volume and density are highly dependent on 
aitaconutente(ventilation) rand bloodscontent (perfusion) in 
therlunge whicheéing turn’ ois related) ito< bloadk@pressure. 
These are much different in the live breathing animal than 


inetheseexciseddulunga (42); thiswetsPethe Vmostpettkelsy 


explanation for discrepancies observed between the in-vivo 


CT and the ex-situ measurements. Since the air content 
waseleontroliled Vin our) texperiments Wby sealing wofiethe 
BracheasQedifierencesiring lungs periustonyiwere further 
investigated. 

ingorder ®toetest whether va change? in’ \bloods presenure 
affects lung density, the following experiment was 
performed. Efenht= mice were anesthetized with an 
intraperitoneal fingection  otfiIrnibromoethanol.) Bhey@adsio 
peced VedRangeinyectiongmwofe.O0<27 mie o8—"a0( €anticoacsulant 
(Heparin Sodium, Hepalean) into the thigh muscle. They 
were then scanned in the usual manner. When mid-lung 
slices were reached, ther ,animals sweres,abledettrom the 
jugular vein and re-scanned immmediately. Mice: «biked by 
Varihous@nanounts, ebutasalhawere still alive atvehesgendhot 
thegexperiment.enlhreercontrol emicet were ‘salsom scanned ; 
these mice received the same injections of anesthetic and 
anticoagulant as the experimental animals, but were not 
bled. Thesertugnicershadugslicghthy eehigher \saverages lung 
densitidsa(O0ssl2te/cmy)yethan theecontrol value sestablished 
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A definite decrease in lung density from baseline was 
produced when blood pressure was reduced by bleeding the 
animal. These amounts varied from a 9% decrease to a 6/74 
decrease depending on blood loss. This effect has already 
beenSreportedsin dogse(41) 55) Therefore, thesalosstof peblood 
dnd@ resultant vdecreasecine blood pressure inviange produced 
whenithe lung is excised’ from the» body; explains ®therlower 
density measured for the excised lung as compared to 
in-vivo CT densities. At present, The CT measurement 
Sives+Stheseamost Saccuratesvaluesftorsine=vivo lune dansity. 
Any method measuring density and volume in excised lung is 
therefore subject to a systematic error because of blood 
Foss durvingewextision.» This conclusion, however, does not 
invalidate any experiments performed measuring post-mortem 
changes because relative (%) changes are usually measured. 
Such changes are independent of the absolute value of the 
parameter being measured. 

An attempt was made to correlate the blood loss and 
the decrease in CT density. Mice were anesthetized and 
bled + frome@asveini nearsthe) eye, whererthe »blood scould) “be 
takenVv@iup: Tinieatesyringee and weighed. This £irsegattempt 
failed as the anesthetized mice showed patchy density 
increases throughout the lung, and it was not possible to 
measure an accurate decrease in density. THiispawas nl 
Contrasts touethe  prevyviousy experdmentsopwhereqe@mice had 


Feceivedeasneparineingection in addition to the anesthetic 


and a slight homogeneous density increase had _ been 
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observed. 


2.2.4 Conclusions 

This partlof thetstudythas indicated thattan imerease 
imetdensityevot mouse? Lungr “causediby irradiationgcan be 
detected easily by CT scanning» The time course of the 
density changes, the time of peak response, and the dose 
dependence correspond well with the histological findings 
CLr6) = therefore YthescT scan is *a good indicatereor the 
€hangessoccuring during @themacutevresponse?’ Of -thestune™ to 
irradiation. imdeedgfa “change* offmore than’ 402 adiuelung 
densityeis “predictivesof > mortality. Moreover, ® Ch @data 
were obtained weekly from the same animals, so that they 
servedvasetheir—own controls during progression “of” the 
effects. This non-invasive technique has an advantage 
over other approaches in which data are obtained from 
different animals by excising the lung post-mortem. As 
wellycT densitometry gives the ‘most accurate value) for 
Lievivowelking Gdenstty sulhis*non-invasive method tcotidebe 
applied to monitor lung tissue damage in patients who have 
received radiotherapy. Pe tfiacel Site thas @galrcadyeepecn 
demonstrated to be advantageous over conventional 
tadvography Binaponis cesard (79,89). In otr Lustitute we 
havewedestoened 2such a, clinical ‘study £Gr patients, 
receiving upper halfi-body radiotherapy. Results are 


presented in section 2.4. 
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2.3 Experiments with Dogs 


22321 introduction 

Theil estudiess, withesemice reported ein’ the “previous 
sectiong* proved suseftlein establishingythe, validteytot CT 
densitometry as a method of assaying lung damage in-vivo. 
However, several difficulties were encountered due to the 
small size of the mouse lungs, including the definition of 
thes¢actuail »lung»svolume. To overcome these problems, 
experiments with a larger animal was desirable in order 
EO: (a) sample a larger volume of lung, (b) study changes 
mear the pleural boundaryroftllungyand Gc)ycoupare togother 
in=vivo diagnostich tests, such as®standardsxXstadiography 
andes nmuctearsaescintieraphy™= hom the intercomparison 
studies, the animal was sacrificed as soon as the dog was 
in the pneumonitis phase and all tests became abnormal. 
Histology and electron microscopy were also used to 
analyse the lung sections in order to correlate the extent 
and) jseverity oftlune damage with’the’ diagnostieliindinegs ; 
ihetoniy Limitations tortheidoe studyawas © thessatestriction 
Lnetthemntimberfor animats thatecouldybe manageds  “Lhus® the 
tinmeLeourse: Pots i pncumond taspjawase studieds sinititadividual 


rather than groups of dogs. 


2 oMeCmleradtatton 
Ten dogs were irradiated with a single dose of 20 Gy 


Come tnesemid—-wuneke innietheninferiorazoneror ther righte lung 
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only. Two dogs were irradiated to a lower dose of 10 Gy. 
The Cobalt-60 tptadiation was delivered with a 
patallel=opposed “pair “ofs 5x7cem*? fields tol Sachteve’ “a 
uniformity**of Vapproximately 102 in lune dose .~8ihe dose 
Was, chosenthisgh®enough sop as’ toe *prodtce way sieniticant 
likelihood of damage in a restricted volume of lung, but 
irradiation of heart or mediastinum as well as esophagus 
was minimized. Tniseisestentiiicantly diftterenGguarons one 


irradiation protocol used for the mouse work. 


Ze oso *Diaonos eictlesits 

The animals were mongrel dogs of both sexes with an 
averages body® wedieht of 921 °ke." “During all procedures the 
dogs were anesthetized with pentobarbital (30 mg/kg) lying 
Supine and breathing spontaneously. Besides CT scanning, 
conventional X-radiographs of the *dog»®elungs wernesealso 
obtained. Pn caddtteton, clearance of radioactive sodium 
pertechnetate (teoy) aerosol from Lung; cand Lunes pernsusion 
using 99mTc-labelled macroaggregated albumin (MAA ) 
scanning were monitored using standard nuclear medicine 
instrumentation. Details of these nuclear medicine 
studies are reported in a separate thesis submitted by I. 
Ahmed (Department of Medicine). 

The dogs were CT scanned every two weeks starting at 
Sueweeks Lpost=iLrradiation. Pre-irradlation@ascans® were 
obtained and served as controls. A digttaleradilLographetor 


iscowt, tscans OLunthemEthOrax was first )»petiormed efor 
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positional referencing with respect to the area of 
irradiation. ADP Lox ima Delve Om slices. each lcmethick 
and spaced 1.5cm apart were acquired in order to cover the 
complete doe lJlunewhronmmapex to diaphragm...) in each slice 
the complete’ lune areapwassoutlined® Stor eleie ands right 
Pune, and the average density calculated therein. The 
best correlation between traced area and true area for dog 
lung was obtained at a level setting of -370 Hounsfields. 
A window setting of 100 Hounsfields was used, although 
this setting @istnotvasmcritical. | The areator interest was 
traced as close as possible to the pleura of the lung. An 
example of an irradiated and unirradiated dog lung and the 
tracings are shown in Figure 2.26. 

In the dog-there is an intrinsic density gradient 
from: apex to “diaphragm= as well. ase) frome anterior to 
posterior presumably due to blood pooling caused by 
gravity. This gradient is apparent in Figure 2.27 where 
Nor is plottedmacmae rc UNnctEi one Ofte diStanucces2,0m apex to 
diaphragm, relative stor theiisternal notcheat, 0 cm. A slice 
near the apex and one near the diaphragm is shown in 
Figure 2.258 tome demonstrates, the density cradient in a 
pictorial manner. The denser areas are more vascular. 

The right and left lungs do not have the same _ shape 
and longitudinal position in the thoracic cage. Therefore 
thes bert)» lung slices at "the same Longitudinal position 
could not be used as controls in these experiments even 


though they were unirradiated. In the interfor zone’, 
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there” isan additional problem®due to the diaphragm which 
interferes@with the*®leftsside**lung* ‘slices, Wicures2e257 
The radiation “dose “received during each" ®scanning 
procedure depends upon the technique used. The X-ray tube 
VOLtage *was \fixed“atel20 kVp} "but "the® tctrrent sand @i pulse 
WidthPewas "varied “depending ® upon “the “size ofthe dog. 
Representative maximum doses measured at the dog skin are 
shown in Table 2.10. These doses were measured using a 


specialized aire ionization chamber CPE2): 


Table 2.10: Radiation dose to dog lung and technique used 


during the scanning procedure 


Dog l Doge: 2 
Current 160 ma 250 ma 
No. of pulses/scan 576 Ske 
Pulse width 3.3 msec 3.3 msec 
scan time 307 mA-s 480 mA-s 
Dose Ze lieecGy JOO CGy, 


2.3.4 Results 

As noted earlier, ten dogs were irradiated to 20 Gy, 
while | two dogs were irradiated to 10 Gy. Pre-irradiation 
scans were performed on the dogs to establish control 
Valweswort density. -lypically the inferior 4 slices: of the 


right lung were through the irradiated zone. CT number 
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values were averaged for these slices only. Lung was 
outlined in the same manner as with the mice, and an 
average density for the slices was calculated. The 
control values were obtained for the same volume of lung, 
prior jfto irradiation. Control values were=alsomoutained 
in a different way by sampling the unirradiated left lung. 
mnus,. gato each) scan) ftime,, 7a) valwe {AN = Nop Cleft) = 


CT 


Nap (Ftshe) was determined. The data analysed using either 
the gipre-trradtation™ ricnt lung -scans.som theaconcum: ent 
left lung scans as controls gave similar results. Only 
the results obtained using the pre-irradiated right lungs 
as controls are reported here. 

The average CT number for the 12 dogs «was ))-755229 
Which) corresponds to a,density of 0.25+0% 0s sc/cme emetic 
value is considerably less than that of miice,* and» “much 
closer to values reported for humans (129). The CT number 
and density values for the individual dogs are shown in 
Figures “2.29 WA B siand 72.3508 tor different -timesmad ter 
irradiations horn the dogs irradiated to 720°) Gy sone) doc 
showed a significant increase in density as early as three 
weekseatcen. tpradiation 9(C-568)).. Another dog (C-694) 
showed a significant increase in density at 4 weeks. Two 
dogs died of unknown causes before tests became abnormal. 
Of the remaining dogs, all but one (C-566) showed 
significant changes by 10 weeks after bucadiatvon. 
Infection was ruled out since dogs were given weekly 


innectioner ots Peni ct line 
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IRRADIATION OF INDIVIDUAL DOG LUNG 
/B DOSE 20 Gy 
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Fig. 2.29 A,B Average CT number a) for each dog 


onadtate dato 20 Gy. 
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Of sthe does irradiated to the lower dose (10° Gy) > one 
showed a significant density increase at 4 weeks after 
irradiation, which increased further at 6 weeks but 
returned towards normal at 8 weeks. After this time this 
dog showed normal results on all tests, and it seemed that 
the, Carly effect had resolved. This supports, thevconcept 
of lung repair. The other dog became abnormal at 8 weeks. 
Both dogs were sacrificed after ten weeks. 

The onset of radiation pneumonitis seemed to occur 
earlier in the dog than in the mouse. Because the animals 
were sacrificed when tests indicated abnormal values, long 
term follow-up was unavailable. Therefore it is not known 
whether the initial pneumonitis phase resolved, andy pli 
severe early changes are predictive of severe 
Somplicatitons sin “thes later portion of the acute 
pneumonitis phase. 

the earliest onset of abnormality as observed von Yall 
bests, (ls  snown ain Table 2.11... A visual, inspectionsor eene 
dog CT scans was done by an experienced lung radiologist 
and visual onset of abnormality was recorded. The values 
for the radioaerosol and MAA perfusion scan as well as 
those for X-ray were obtained from I. Ahmed (private 
communication). Im this study, the objective was ton isec 
which. test.was the earliest indicator of radiation damage 
Ceol ip. It is seen that a visual inspection of CT images 
is advantageous because of overall pattern recognition of 


damage. However, visual CT is subjective and varies with 
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experiences @andseimare* viewing Sconditions Ceses level). 
Moneovertiiepice not® userulperoroequantitative® Steltow—-up 


studies. CT densitometry is more objective since it uses 


babsoluteosch, values: éburrsttrict@criteria Lor stne, steacineg 
procedure, level and window settings and boundary 
definition must be established. This method of averaging 


CT numbers, however, dilutes the density increases, which 


Can oceurminevery Localized’ patchés or "near? *the® ap leura< 


TableeZ~lischaritest onset “of abnormality (weeks) in 
dog lung as observed on various tests 


(M.Se. 8thesis L. “Ahmed, “U~of A) 


Dog Radioaerosol MAA Xray CT 
(clearance) (perfusion) visual density 

C=566 2 9 8 6 = 
€=5726 Z 6 i 6 6 
C=577 i i 10 5 8 
C-694 2 9 9 4 4 
C2620 Z 14 13 10 10 
C-564 6 tl ko 10 ING) 
C2568 Z 4 4 3 3 
D119 3 6 6 A 4 
D-134 4 6 v7 7 = 
Mean 2.67 8.00 Sie zz ayemey) 6.15 
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The Cl seani, ts expected to be superior to conventional 
X-radiography because the superimposition of structures is 
avoided. GI "is “sable “to “detect smaller difterences: in 
density, both of —smalii. low-contrast. "rtocal lestons mand) Of 
more global increases in the density of the homogeneous 
structure. Therefore obtaining an average density “for a 
siice of lung should detect such changes. 

The earliest indication of radiation damage was shown 
ii a, Change: Sine (clearance of radiogerosol.|§ Vasual and 
quantitative CT showed changes 3 to 4 weeks later. The 
MAA perfusion scan and conventional X-radiography were the 
last to indicate changes in lung at roughly 2-4 weeks 
after ihe Cl resuitse. In this study the MAA perfusion 
scan was not an early test for radiation damage, as. had 


been expected from’ the work of Prato’ ¢36)- 


23 oe lamt tat ons 

As shown previously, the density of dog lung is more 
Variable than in the mouse or human, with a eradient 
towards higher values from apex to diaphragn, the lower 
lung being highly vascular as shown in Figure 2.28. This 
makes the tracing procedure more difficult since we need 
to exclude large blood vessels. 

Radiation changes were observed mostly near _ the 
pleura, and may be excluded by the CT averaging procedure. 
Therefore, the CT number averaging procedure used here may 


be a conservative estimate of true damage. Furthermore it 
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isvdiriilculescomiudge  pogition of the irradiated= zone trom 
Scan to scan Bina dynamic orean, such as the lung. Skin 
markings mark the area on the dog, but the lung seems to 
wander in and ouvtwolernts marker Pierre...) Lrradiation 
of the entire dog lung would avoid some of these problems. 

Another unknown factor is the Short-term and 
long-term effect of anesthetizing the dogs weekly. It has 
been pointed out by Hedlund et al. G41) Bthat.s atelectasis 
May occur in lung of anesthetized supine dogs breathing 
Spontaneously. A long time was required to perform all 
tests so that the dogs would be scanned following these 
tests after being under the influence of anesthetics for 8 
hours. There are probably some effects due to the 
anesthetic, observed occasionally as dense local patches 
im’ the Lunges suggestive of lung collapse. Anesthetic 
administered) justapriotreto Cr scanning shadpepeLoduced  Ysuch 


an, effect invmicemCsect donee. 2.3) 
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Two uot tar adhocmapnmeG scout scan.) of “dos 
thorax. The boundaries of the radiation 
Eleld@are marked by positioning clips 
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2.4 Clinical Study in Humans 


2.4.1 Retrospective Study 

In order to determine whether a densitometry study in 
humans would be weefud§ 70 ithorax’CT scans@of patvemes at 
this Institute were analysed retrospectively. Large blood 
vessels and obvious. carcinomas in lung were excluded 
during the tracing procedure. The average CT number for 
lung for 52 patients (of both sexes, aged 50 to 80 years) 
was) f0und to bey O03. 5434.0 H which comvesponds: toegam lunge 
density of 0.197 g/cm>. These patients were scanned while 
hotdine their breach atjrull) inspiration; this reproduced 
the lung volume for serial scans and also reduced motion 
artefacts. For each patient there were usually 10 tto 17 
sfices, each 1.0° ecm thick through lunges. These deneditives 
are 12% higher than those published by Van Dyk ClZ9)) £or 
this*® ace = group at full inspiration. Sihiis may  bemdue to 
the fact that Van Dyk excluded patients with carcinoma in 
Lune, whereas in our study these patients were included, 
although visible tumor was circumvented during the tracing 
procedure. 

Three cases have been selected to illustrate how 
qduanGieative ) Ch “compares=with other tests such asc yvisual 
imtenpretation | of. the, CL scan by a radiologis &. 
conventional chest X-rays, and bronchoscopy. All these 
caseqc aad etumor. Ine lone. Relevant Sellinical jdata, eis 
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before and after radiotherapy were available for each of 
these cases. 

Case lL: the tiret, Ci scan. was done 1.5) weeksaipt.0o. 
to radiotherapy. The second and third scan were obtained 
8 weeks and 21 weeks after a completed course of 
radiotherapy. CT number data, plotted as a funmct7on) of 
Yonga tudinal position in the thorax, are. shown inwecieare 
232 9A, BaC. theme oriein is) at the sternal. 10.¢ nae line 
Hongi tudinal border of the radiation field is iandiecatedsin 
Pi SULre 2.52. Oa Lie lateral boundaries of the taqtation 
field nave also been taken into account. Thus only) that 
part of lung which lies within these borders was sampled 
in each slice to obtain the CT number data. 

inwthespre-trradiation sean (Picure: 22.325 "A Sidemadny 
valwes “ins the Veit lunge and lower right lung arelnormad.. 
Density values in the upper right lung (which corresponds 
to thes position fot the tumor) are, higher < ObvLousiisolid 
tumor patches in this area were excluded during the 
tracing of the lung. Eight weeks after radiotherapy 
(CPaAcuren2no2 BB), theres is ancinerease ine densi tye eine tue 
left Lung throughout the area.of the radiation £ilet@. Vin 
tiiewecichte tune. at) the position of the tumor, thenveme cea 
decrease in density in some slices and an increase in 
others. Vatacin, Luts saned pica dLoriLcoul st to interpret as 
there is radiation damage superimposed upon tumor 
regression. However, there is a definite radiation 


reactionnine the opposite normal left lune. These findings 
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Lung Density (g cm-?) 
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are confirmed by visual interpretation of the CT scan. No 
recurrent disease was detected by bronchoscopy. 

Twenty-one weeks after radiation therapy (Figure 2.32 
CG). the radiation change in the left lung has resolved 
somewhat although density values have not returned _ to 
normal. There is a further decrease in density in the 
raehe0 Lune, indteatings further clearance of disease. 
However, normal density values are not observed, 
indicating residual disease and/or superimposed radiation 
damage. These findings are consistent with the 
Padiologist. s)interpretataon of the CT scan. 

Case 2-5 The first ‘scan was done one “week “prior to 
radiotherapy. The second and third scans were obtained 5 
and 20 weeks after the end of radiotherapy. The CT number 
data ares plotted: aseal function of Longitudinal postetonein 
Chew BNO a xe ities te ipure. 626359. “A BCs Five weeks after 
radiotherapy, Crynumbers, are sitehtly higchersin them. ont 
lung. The tumor in the left lung seems to have cleared, 
Duibsthewluneetissues1seot shagher density (0.25 92) cm )etnan 
normal, indicating possible radiation damage. No tumor is 
detected “by ~visual “inspection of the CT scan’ or on) the 
chest X-ray. Twenty weeks after radiotherapy the density 
throughout, the > rieohtelung remains the same. However, in 
the left lung there is an increase in density at the 
Vast don or tne primary — “tumor and within the radiation 
Eveld, escelvatltvyer to, the previous scan. This scould be due 


toaradiation chauge and/or recurrence of tumor; on visual 
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Lung Density (g cm-3) 
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inspection, this “haziness" is suggestive of radiation 
changes, but residual disease cannot be ruled out. Thus 
the CT number analysis is consistent with other tests, 
pParticitarly wthewm Cle cimape inspection. Comparisons to 
previous scans, however are essential for interpretation. 
The, visuad- inspection "of the. CI ‘image ®is ‘subgecrive, 
whereas quantitative CT is more objective. it) provides 
information on the extent and severity of lung damage. 
Case: There is one clear case of radiation 
pneumonitis. The diagnosis was made by the radiologist 
trom themeCl = scam Gand= confirmed by «a pathologist at 
autopsy. This patient received radiotherapy Ore 
epidermoid carcinoma: of the Lung to a dose of 30 Cyminee 2) 
i eacte Ose -OV er a pertod of 35 days.5 Iwo Barges fields 
(U62.50%17 5 em] ) were wsed inttially in a parableleopposed 
arrangement. At the end of this treatment, a “boosting” 
dosenor LO "Gy ine Oemtractions “over 7 = days." iueince = two 
smaller Siteide s@l6xt4" “cm )ewas delivered. Iwo CT seans 
were available for this patient, one before irradiation 
and another five anda half weeks after treatment. For 
these two scans, CT numbers are plotted as a function of 
longitudinal lume position (in Figures2.34 A,B. Before 
ivnad@tattron, theGl number values in regions with no mtumor 
Tew che normalearance, buteatter irradiation Cl numbers 
within the .rradiated zone increase dramatically. This 
ertect eis conranedmtuo the radiation field as: indicated in 


Figure 2.34.82, However, an) unexpected result is also 
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observed; lower density develops im the wnivradiated 
lung. This is “probably due to a compensatory “effect 
causing the undamaged lung to overventilate. The patient 
deceased eight and a half weeks after completion of 
irradiation atpewnpches times "all “clinical. ‘sympiomice of 
radiation pneumonitis were present. The appearance of 
radiation pneumonitis in humans is shown in the CT scan of 
Figure 2.33 ee lhiseparticular,, case is” an @vextrememmcace 
where radiation pneumonitis was clearly visible and CT 
number analysis would hardly be necessary. Nevertheless, 
the severe changes were apparent three weeks before the 


patient's death. 


2.4.2 Prospective Study 

We expect CT number analysis to be more useful in 
cases where the entire lung is irradiated and diffuse 
changes may occur affecting the entire lung. Therefore a 
study Mewasmeeinittiateds inesorder ~to # monitor Anne densi ty 
values, 2taswarious= Sintervals ay patients receiving 
hemi-bodyemeirradiation at this, Institute.) Thespativents 
used in this study were diagnosed with small-(or "“oat") 
cell, lungs “cancer which is one of the most aggressive and 
bethal stumours “(l22). 0 .the treatment strategy may iAanetude 
chemotherapy, systemic and local radiotherapy, or 


Compianations “thereof Cl 21122). 


24a e radiata ome CUHBL) P.rotoco. 
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Fig. 


Jes Gl scan Of (human lung CA) normal lung 


(B) radiation pneumonitis in lung. The 
radiation damage is confined to within 
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In-order to obtain the large field necessary, for 
uppet Nat itobody sirradiation ~the patient Lies on gaecoucn 
Near the floor at a distance (SAD) cf Z00cm “from the 
source. Arms are at the side with forearms crossed on the 
chest. The field size is chosen large enough so as_ to 
enclose the | entire upper body superior to the umbilicus. 
Dose is delivered with parallel-opposed 6 MV X-ray fields 
along the anterio-posterior direction. The 
anterio-posterior separation distance is measured in the 
mid=sageital plane of “the  ‘spatient cat 10 semeinteryvais 
starting from the umbilicus. A mean separation is then 
calculated, and is used in dose calcudations®. A shield is 
used for the entire lung, thus allowing generally higher 
doses of radiation to be delivered to the disease sites. 
Ties tryin al: fteld? sd zieu ds) 40/5240 or picmel,» » atidesm hel wotae 
diosies 10. 12:0 Gy? Hsecdeliivered iny 8 fractions; secthel Luncisiose 
Gumcon rec ted; liseer section: Ss. 2ei0fs thisrthesis).,, ismkeptagmat 
obeGy/6) Gracthonse at ithe: position: of mid-plane. Thewlung 
shields are beam divergent blocks of “cerrobend” 26mm in 
thickness Gk HVL=15 mms s comsittucted) «to © matichy ithe 


indiwddualapatdients by- shape. 


Lewis 2 lt Denes tometry 


Scans were obtained according to the schedule shown 


Lane Pde uirrer 215: 6r. 
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Fig. 2.36. Treatment schedule for hemi-body 


radiotherapy patients 


The actual time of scans deviated from these nominal times 
because of difficulties in having patients comply owthic 
schedule; this group of patients is very ill, and some 
deceased during the course of study. 

During the scanning procedure, five transverse 
slices, each 10 mm in thickness, were obtained through the 
thorax. The CT number data are recorded and averaged as 
described cipal the animal experiments (section 2.2). The 
sternal notch is the reference landmark and anatomical 
Orie atone. mM)’. The CT data are presented in Figure 
Jee Ce control, density value dis that. “given for) sthe 
pnormal “patients of this age group as published by Van 


Dig Ke 2198) 6 The. “dataw are. plotted. for the> individual 
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patients, with Hightweluas (KR) “and lett Vlunce cl). 
separately. ilwoltoft@theypatients, “(S)) had rsungeryeas twalid 
and one notes that the remaining good lung is of lower 
density, presumably expanding in order to compensate for 
the missing lung. 

Forrvall patients there isan “initial” decreases in 
density which» is attributed to clearance of thevprimany 
disease, as seen repeatedly in the retrospective study. 
Thereafter, all values return to normal and remain there 
for the duration of the scans. This is an expected result 
Since the dose delivered to lune (6 Gy Cuncorrected) ins 
Evactious),. ls not expected to produce radiation 
pneumonitis. At this dose level, there appears to be 
dittie synergistic effects of radiation and chemotherapy. 
Since this study was completed, the dose to lung has been 
increased to 10" “Gy 9 in. 8 fractions “by  “decreasingsgthe 
thickness of the lung shields for all subsequent patients, 


With the intent of improving control of the disease. 
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2.5 Summary - CT Densitometry 


The CT densitometry studies are summarized in Table 


Table 2.13: Summary of CT densitometry studies 
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Species Number Sample Volume Doses Observed Change: 
Mouse 142 cevalal Die say 15 - 24 weeks 
Dog Ja) lower right 10>. 20) Gy 4 - 14 weeks 
Human 3 irradiated zone 60° Gy/30 fractions 5) weeks 
Human 6 all 6 Gy/8 tractions none 


From the mouse study one notes that CT densitometry 
is sensitive enough to detect the histologic changes 
observed in lung during the acute pneumonitis phase in 
animals. The time-course and severity of changes 
corresponded to times of observed microscopic changes. 
Even though the mouse lung is very small, only 0.4 ml in 
volume, significant CT number changes were observed before 
any visible changes were apparent. The CT number changes 
quantify radiation damage; a change in density of more 
than 40% was predictive of death in the ensuing two weeks. 
Not ‘enough animals survived this phase, SO > that.) «che 
fibrosis (late) ‘Stage "or lune damage could not “be 
observed. Due to the large number of animals in the 


study, good statistics were obtained. 
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From the dog study it appears that radiation 
pneumonitis occurs at earlier times in this species than 
inieeive: mouse Due tosthe fact that the dogs “received fa 
Rkeh Local @cdose™! torthe Lower “zone of the riche dung, Cl 
visual changes were observed as early as CTI number 
changes. Radioaerosol clearance was an even earlier 
indication. But since there was no long term follow-up of 
the dogs, it is not known how well these early tests 
quantify later radiation damage. Since many cell types 
are involved in lung damage, it is not clear whether the 
Celigthat “indicates angeariy change in “function willie be 
Biers (celle) type finvolived in) the critical maltuncevonimner of 
therilkung idurinugethe clinical manifestations of “radiation 
pneumonitis. Lome” “term follow-up and irradiation, to, the 
emepretkbunge would berdesirable in order to investigate 
Which @ 8b the tests most accurately predicts a severe 
effect in the late portion of the acute pneumonitis phase. 
The most useful test for radiation damage to lung may not 
Demetiveteariiiest™test “ut ethe quantitative test which ewr el 
give information as to what level of damage will predict 
severe complications. In addition, we wish to obtain a 
dose-response curve so that doses’ to lung can be kept 
within safe limits. 

In the human (retrospective) study, eae number 
analysis) yielded quantitative information about radiation 
damagertasiawell eas clearance or recurrence of tumor, and 


these findings agreed well with independent tests. In the 
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human hemi-body (prospective) study, CT numbers indicated 
initial clearance of disease and remained within normal 
limits ™thereaftter. This is as expected, since “dose to 
lung was too low to produce severe radiation damage. In 
subsequent patients, dose to lung has been increased to 10 


Gy ShEtactionssinworder to improve the “therapeutucmratioe. 


of disieasescontrol, to complication. 


FUTURE WORK 

Since much work is being done in mice, and radiation 
effects in mouse lung are well documented, it would be 
desirable toe follow “av “sroup- of mice:  ipradiared to 
sublethal doses into the intermediate and late fibrosis 
stage and continue the scanning. Such avVstudy is now, )in 
progress: (Dri. GCG. Miller, private communication) < 

For dogs,“whole lung arradiation would  taciilicate 
interpretation, “and Vong “term “follow-up miehe sindiicate 
which test quantifies relevant radiation damage, and 
whether the test showing a severe early response predicts 
the later response. A group of controls on whom all tests 
are done at the same time as the irradiated animals would 
be necessary in order to show what effect the weekly 
anesthetic or the tests per se have on the dogs. 

faethe humans, doses are kept low in order “to ~avodd 
fadiactvwons, damape to lung. ibe study is ongoing with new 
HBL patients receiving higher doses. But many patients 


are. quite’ ill and very few complete the study. Possibly 
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patiemts a ceceitving iMmuchws higher doses locally; using 
muita field irradiation, <could be monitored: 

From this work on CT densitometry of lung, the ideal 
experimental “animal "tio uweie imilght “be the (natyefiomabne 
following reasons: Lys The lung of the Cait is 
sufficiently large and has a nice homogeneous appearance 
Onwetne Gly scan asashiown iui cure e2). 25°: ine Rater Teiou ld 
probably be well immobilized» cfor «CT ‘scanning wi Choate 
anesthetic. 3. Facilities for keeping rats texts t tin many 
laboratories and many animals could be used in the study 
Worlobtaim-datatwith vcood statistics: © Asi well, tefteces #05 
Eracdtionatifon, chemotherapy, and radiation response 
modifying drugs need to be investigated to obtain reliable 


dose-response curves. 
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3. RADIATION DOSE IN LUNG 


Icom introduction 

Inémany, radiotherapy: centres;«radiation fieldse which 
encompass one half vorrtthes entire body aresused) to) treat 
systemichor,disseminated= cancer. Total) bodyasdrradiarion 
(TEID)teissused prior!’ to” bones marrow! transplantatdony inthe 
treatment of acute leukemia. Radiation doses to the total 
body in the range of 8 to 10 Gy destroy leukemic cells and 
also suppress the immune response in order to avoid graft 
rejection. Body doses of less magnitude (1 to 3 Gy) are 
used to treat radiosensitive malignant diseases such as 
lymphosarcoma, Hodgkin's disease, and Ewing's sarcoma (6, 
b2e>o).@ehaltiahody arradiatton, a(HBIl)eeieetused? “Comm temcat 
widely disseminated and advanced disease (120,121), or as 
adjuvant!’ therapy for tumors with poor’ prognosis-« In these 
cases, theemain result is rapid relief of pain; §sometimes 
withtnw24 hours,sinepatientsiwith severe 'paingwhich cannot 
be* (controlled even with»high doses of analgesics. )This 
patlidativeybeneftitioftenteontinucsaforevtheseremadnder: Sor 
the patient's life. in addition; the tumoreand’ metastatic 
burdényis decreased, thus resulting» in» improved ‘squality 
andsequantityaegol isurvivalaetortethe “patients SOF Ghar, 
elintiealilresultsearesprelininary for*)both) TBLleprior i ttso 
bone marrow transplant and@viors SHB igepUuUteathey are 
encouraging, (275100). 

Different treatments techniques: arestuseds in the 
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various radiotherapy centres. These differences occur in 
the energy of radiation used, the patient positioning, the 
dose prescription, the dose rate, as well as the use of 
beam modifiers such as tissue compensators and shielding. 
Therefore direct comparisons of results are difficult. 
The techniques used in North America and England 
(reference 59), and in European centres (reference 6) have 
been reviewed elsewhere. 

Patients receiving TBI or HBI treatment are often 
sick and uncomfortable; it is desirable to give the 
treatment as quickly as_ possible. In bone Marrow 
grafting, the patient should also receive the marrow 
transplant as soon as possible. Thus, a single dose at 
high dose-rate (eg. 200 cGy per minute) is more 
desirable. However, from radiobiology studies this is not 
optimum from the viewpoint of normal tissue recovery. The 
maximum therapeutic ratio (i.e. “Rill” -'ef i *tarcgetd@tcett- 
relative to normal tissue damage) is the goal. 
Radiobiology studies indicate that the therapeutic ratio 
may be greater for irradiation at low dose-rates or for 
total dose delivered in multiple treatment sessions 
(83,84). All of the major organs subject to radiation 


toxicity during TBI ‘and HBI- have. ‘capacity to’ «repair 


sublethal damage between dose fractions; this is often 
not the case for the target leukemic or cancer cells. In 
upper HBI, the lung is the dose-limiting organ. 


Clinically, the discrepancy between practicality and 
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radiobiological findings is only being reviewed now for 
doses fractionatwon inv large field radiotherapy ©CYZ0p- 

Bor eEBl tow be “ettectives “a> “sufficiently “nvenh and 
uniform dose should be delivered to the entire body. 
However, significant dose inhomogeneity will occur because 
of variable tissue thickness and composition throughout 
the body. Therefore achieving a homogeneous dose 
distriibdthon  eoversthe entire body is a difticuwlt spiyctcal 
problem. Sometimes it may be desirable to obtain higher 
Orw lower” doses to various parts of the body depending von 
Bue locationsol critical or previously—treated onreans. 

Generally, conventional treatment machines are _ used 
for these treatments. In order to produce largesfticids: 
extended source-to-patient distances are used, while dose 
uniformity is achieved with parallel opposed beams. 
Acceptable dose rates and dose-uniformity are obtained 
Whol "nich  Senersy X-rays produced in linear accelerators 
Ghinads) 7 Cy). In one eiustitutes a <specialized™. Cobale—oo 
machine “Shas = beenyedesigned specifically tor TBI and HBL 
€b30)).. 

Ror Iargersfields and distances, the (dosimetry fdata 
obtained for smaller fields may not be appropriate and 
musuebetverified (124) Discrepancies have been reported 
because of differences ineepackscatvtrering.| selectnon 
ecntaminatton,Gdose ©Lratios (erg. “tissue-air ratios); and 
corrections for tissue inhomogeneities. The following 


parameters are affected: 
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(a) Inverse square law 

The inverse square law may not be appropriate because of 
substantial photon scattering from the surroundings (e.g. 
Bloor) tand)  dose-rates must be verified vat the large 
distance. 

Gboelissue-ain ratios and Tissue-=phantom ratios 
Eheseindependence of TAR (tissue=-air ratio) and TMR 
(tissue-maximum ratio) on distance to the source should be 
Vewmliied. | sTAR sgand | IMR s "are usdally measured sin, va 
semi-infinite phantom: in Large field irradiation, the 
radration (field tymay -be larger than the body.) “and Sia 
reduction in dose due to scattered radiation occurs 
Pevatal vewito Biherigiwiliiatscemi—~intinite phantom: According 
toy) Faw and Glenn (26) the dose depends on the field size 
or the phantom size, whichever is smaller. 

(c) Electron contamination 

Electron contamination of the surface dose due to 
eobEimetonms cand 9‘shielding block trays will be Less of ta 
problem since the patient surface is further away from 
them. However, subpiaces doses, may “still be@sitenimcant 
because of scattering Brom surrounding materials, 
fvcluddie backscattering» and scattering in the intervening 
volume of air. The position of maximum dose Sac should 
be verified. LiETsutitetent backscattering material ts mot 
provided behind the patient, the exit dose may be low and 
that part of the patient may be underdosed (32). 


(d) Inhomogeneity corrections 
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Sineée InTHBE or TBI the complete lung will be irradiated, 
lung dose corrected for tissue heterogeneities must be 
determined accurately for individual patients. Die eis 
crucial because the radiation pneumonitis syndrome is 
Significant when the total lung volume is irradiated to a 
Hach | Pdose@sin s @aasing le Cotmea (few )merracttonGs). sells nas 
proven, totbesktheimajoreGand@ fatal) toxicity for uppers eibd 
CLOO) < A slight increase in dose produces a dramatic 
inerease inwincidenceseot weradtation pneumonitis, “but va 
sufficient dose must be delivered to achieve _ the 
Ghenapeutie. ehiect. There are several computational 
methods available to correct for inhomogeneous tissues and 
thesia will beediscussed “and > “validated in» sections, 3225. 
Dose perturbations due to bone inhomogeneities at 
megavoltage energies are usually small and neglected in 
HBI. 

AL©Our@institute> “HBL “is sperformed sat. “ay sources. to 
Patient. Surteace distancey CSSD) of 200ncm andawi theagsi old 
size dependent on the size of patient. Dosimetry 
measurements reported in this thesis were performed for 
fteldgsizes pap sto!’ 60x60% cm2e and») at small iand= vlarge 
distances *®e*frome’the® »*source. A 6 MV linear accelerator 
GSiemens,Mevatron’ V1) ise used for the “HBI «treatments: 
Therefonet most’ tof thesedosimetry was performed on this 
treatment unit. Some measurements were also performed on 


a Cobablt-60euniteas¥a reference “standard” radiation. 


TS 


Pin | 4 A or 


bared bases os ews sala eo 7 
oe tae vei kane gereiad: auester 402 baysenes 
ep weet , esad2eq airs. 
sgosherve  altiauaesng avigettse: sif3. sauhoed 


at 
jeer gt 02 sealow geet Leres) 943 43. wade foavs 
7 ;ciscitcaw? (et « ae) ataate & ‘focde J 
oy 75% Yelatzed OFeset ben) 70hee a3. ad)w2 nevotd. 
sewbésy, peob al sese3cn! sdpiie “s) af ¥, 
anaveng o#¢lsatbat. 16 ssasbiont. et seteront 
ag Seseriiod -o4 aue? euob smatoriten | 
r laroveG 924 63527 -ayetts ‘Aieoeguaens | 
passveandeal gui 4587300 2 sideiteva vein 
FPegabsadaye” bar baeevsaih od iriwi : 
‘oi gage ‘as enh anp)Téas¥sseQ 
fp ‘ ico Shane gbfeuee ese eetgushs wes 


iv 
ree ¢ 

7 
ré hauws60968 a? 248 sesurisest eee 
oo a2 GbS/te (G28) @2neseth seaieethnes 
TL se snnssag jc eels (‘S22 €o sosbwegeb a 

i 7 7 
pd, Meuseieeg Stee ae pues eta i. aaeeesey v . ee fa 
ayes bee Sthee. ye has eo ‘Bande $a ge 
= | oy 
te laxelepoReeenes VEY ‘a eeren | S os% 
; io : = tt eons es be 
esngn isn? 124 (ste ae aa Ripe eM | 
erat os sae tous. of 


aa ne 


7 — 


3.1.1 Linear Accelerators 

The (major ‘source sor radiation for radiotherapy — since 
1951 shas been the photon emissions (1.17 and 1.33 Mev) of 
radioactive Cobalt-60. However, in the last decade the 
higher energy X-radiation produced in linear accelerators 
has predominated. The penetration characteristics of 
higher energy beams are better suited to the treatment of 
deep-seated tumors, while sparing surrounding normal 
tissue. Ideally, such a tumor should be treated with a 
radiatlongaose that e1gezero sate the  “surtace, Bis es tow a 
Maximum Gat “the: tumor location; and falls off capidly 
beyond o-. hor ans “Opposed | parallel “pair om tadiation 
beams as used in TBL and HBL, the variation ot dose with 
deptneiunma 82 Che patient, “for: sevemalyeote Our eared t ment 
machines is shown in Figure 3.1. The more uniform dose 
distribution is obtained with the higher energy machines. 
The location of maximum dose is deeper for increasing 
energy. his is caused by the “~build=up of electrons set 
in motion by the incident photons. The dose at the 
surface and at shallow depths is affected by low energy 
electrons produced by the beam-shaping collimating system 
and any beam modifiers (eg. compensators) placed in the 
beam (5,04, 72,100). thesdosesat deeper depths is ¢reatex 
for higher energies “primarily” because of the decreased 
attenuation of the primary beam. 

The high energy X-ray beams used in radiotherapy are 


produced when an electron beam is accelerated to high 
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Fic. oul Vantatnonyotedocseswilthe depth inva water 
phantom) for ay parallel—-opposed pair of 
radiation tields (30x30) cm*). The 


uniformity improves with energy. 


Ley. 


velocity and made to strike a target where it produces 
bremsstrahlung photons. The linear accelerator consists 
ofvan electron guny which injects. electrons “into “the 
accelerator Space, “and “a “suitable target Ghor “= ray. 
production. The source of electrons can be a hot filament 
or cathode in an evacuated tube. These electrons are 
accelerated to a final energy of the order of 4 MeV to 35 
MeV) ing medical. ~ linacs =. Details of the accelerating 
principles have been reviewed by Karzmark and Morton (57). 

On emerging from the accelerator, the pencil beam of 
electrons is bent by an achromatic 270 degree bending 
Magnet, Figure 3.2. The beam is thereby focussed in 3 
dimensions) onto a very Gsmaliy “focal pointy... Ltecan be 
used asrsuwen) after scattering on a foil or raster scanning 
tem ispread the beamy Hntom ab univorm electron iiveld.. “A 
retractable X-rayatarcet can also “be inserted nto” the 
beam for photon therapy. These photons produced by high 
enereayvelectrons are primarily: “emitted sain “the "forward 
direction Srelative, to  sthe incident felectrvoneg beam and an 
inhomogeneous field intensity is obtained. Chase pmo Dem 
is) “overcome by placing a) stlattendine weil ter dan the x-ray 
beam, as shown in Figure 3.2. If properly designed such a 
filter will produce a homogeneous photon fluence. 

The “output"™ dose rate is monitored continuously by a 
system or Louization chambers. This system verifies field 
symmetry and also dosimetry to shut off the machine for a 


pre-determined dose. A light. and mirror combination are 
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used to indicate visually the area of the treatment field. 
Ae qualityieasstrance program isi necessary to veriiy that 
these systems are operating reliably for clinical use 
Cl oe Ole 

A primary collimator limits the maximum field size 
Wotensm@eane Dew produced) for X=ray  therany.= The actual 
treatment field size is further determined by secondary 
collimators, “consisting “of tour thick metal blocks -often 
made of tungsten. There “are -support “trays ewhich are 
designed to hold beam modifying devices such as wedges, 
tissue compensators, and shields. 

The Siemens Mevatron VI linac used in our experiments 
is powered by a 2 MW magnetron, the accelerating structure 
is 0.82 meters long, of standing wave type. The isocentre 
ioe -2aerelLoo> cumands the: maxinum tieldwsizgerat  bhise distance 
fsac0=40 "em G57). Ine this) @worke slLarcer ffi e@lds@awere 
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EngapCopalt—-60fundt used Aniradiotheraoy” the» ‘source 
Diaeldtgation Vs a radcToacrive 1 sOtope., pALcypical cobalt 
SOULGCOECORLAINsS approximately 278 TBq in “a volume of 5 
cml.) BAand@eproduces= eau dose Srate sof? approximately 200 
cGy/mans at SO cm from the source. This source is. placed 
atthe centre of a Vead—-titied steel container. During 
treatment the source is moved opposite an aperture so that 


the radiation beam can emerge. Various safeguards have 
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been developed so that the source will always return to 
the®"shteflded +oftf*s "position, sin the event of a machine 
failure. 

Collimators are used to define the field size and a 
ELAy PecoPriioice beans mod41 1 oT sei seatsotavaiwapies in. the 
experiments reported here, a Theratron 80 and a Theratron 


780 (Atomic Energy of Canada) were used. 


sabes Dosimetry 

The objective in radiation dosimetry is to measure 
accurately @eradt ations dose, at <a) point, im. 4a) phantom oc 12 
ain @tonteabiveatcion. The dosimeter should. ter “small, aa 
order to minimize the perturbation in the medium caused by 
its presencents): The instrument of dose measurement used 
in oun Cexperiments, was” primarily the alr-donization 
chamber. However, for measurements in a humanoid (RANDO) 
phantom, this dosimeter is impractical. Thermoluminescent 
dosimetry (TLD) capsules were used instead because they 
can be inserted into small cavities. Both of these 


dosimeters will be described briefly in this section. 


Ae bes eakalonizatton*Chamber 


Radiation absorbed dose is defined as 


De = dE», /dm (1) 
where Da = absorbed dose in the irradiated medium 
dE /dm = the mean energy imparted to a small mass, 
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In radiation dosimetry using ionization chambers, the dose 
to the medium is determined by measuring the ionization in 
aheaspeibbeds cavity * introduced into the medium. The 
principles, involved are ithe Brage-Gray theory Tand its 
refinements (4). The Bragg-Gray equation, appropriate for 


“small” chambers is: 


Doe=.Q/m (W/e) eS (2) 
where: 
Do = absorbed dose in the irradiated medium 
Qie— tient zatiougcharsce collected 
m = mass of gas. 
W = average energy expended by electrons in the gas 


to produce an ion pair 


Sersacharnge- iol; theselectron 


003 


= ratio of the mass stopping power of the medium 
to that of the gas averaged over the energy spectrum 


For the telectrons crossing the cavity. 


In order to determine the absolute absorbed dose _ the 
= Jil 
diantbities /’m, WW, teiy and aC must be known, and the charge 
liberated, Q must be measured. The mass is obtained from 
a knowledge of the volume of the cavity as well as the 
density of the gas during the measurement. Generally, the 
accuracy of ion chamber measurements for megavoltage 
photon Madtation ts 2.372 (66), while the precision is 


ex@el lent 1 ¢€<0452Z ).. 


22 


F : oe _ ; 7 7 
_ ; - - 


7 ; a a 


gedbrnilt diedidity WEES eedaal: geten aakeck, 
vi pedgectiel sft gettelene edenkwreveb ot” } 
wit obscene oe) opal) pesabapeet qotver Se, 
2! ter yadda. gewiepgasl | ere “Ste esyIevat 982 
» os, teeaged -kubtaspe geeeggie’ wit, a (Rete : 
sok e¥etmeds “Ils 


so 


“ye (eG) | ole ee 

evitea basatunsti ode ct aeob Bedzoads a 
betoalias agtada notseatnad: © Ry 

28 topadenaee 

J+ o) veovoSais. YR OUbnegrve Teter ae 
viag nol of ar AD 

aoviseie 203 2. 7) 

Los ups te. tle goaggode ene, aria Beers a) 
166) Sg%ege edt BS¥e\FOSTETS GEE aA3. Io sens 03" 

es eueg G87 grteaat> epossaate; Sn; . 

y °° a 

wc oth badeegda ,wvebeuds tj osnt 
eeQudc evs eon gtwent of tenn bat ae oneal: 


— 


a 
¥ 
- aa 
“:.* : 
a 


a 


- 


eoo4 bewtset@e @2 ween. s#T «bh 
of? at tisw wa retwes she V6 7 
aaa «i Peers 
oe acl . mht ; i. 


a 


7 nals rene 
a Roe 


Two (or more) electrodes of the ionization chamber 


. 


define the “sensitive” volume of the chamber and collect 
the charge produced by the radiation. Le ‘these radiation 
eCous ists "or \sphotonus, these first interact witn) thewwa il 
and gas of the chamber to produce electrons. These 
electrons then imteracet through numerous Coulomb 


eollisions within, the schamber “cas -and ‘strips ‘“ellectrons 
from the gas molecules. The result is the production of 
positive sons ands cnock=on- relecuronsm wht cia mc ata tran 
themselves to other molecules which become negative ions. 
The electric field across the chamber electrodes collects 
these ion Dads, thus producing “ay “curpene. The 
electrometer attached to the chamber then measures. this 
current Cdosemeracep “or “thel ‘total scharce produced sine. 
certatineperiod) Cdoseo, as, in) Figure 3.35 

in practical use, LG sy elas Lest to Cala prate 
Lonization chambers) in “a known radiation field such as 
Gobalt-—60 radiation [rather than “to determine “all =the 
faceers of equation (2). After —such calipration the 
absorbed dosevat' a point in a medium irradiated with 
radtatton of quality A’ may be expressed as (49)): 

(D Dy = RN Nop Cy (3) 
where 
hrs the Lontaatton weadine of the chamber—-electrometer 
Daal, 


N “Galipraction factor for thes reference radiation 
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CusuallysCobalt—60)), 
Nop: temperature-pressure correction for the unsealed 


chamber. 
The, factor Cy includes: all other factors necessary tio 
eonvert the reading for the reference radiation into 
absorbed dose in the medium for the radiation energy being 
used. ihe joverall uncertainty ine thissdetermination oF 
absolute absorbed dose is 2.3% (49,54,66), mainly because 
of this Dae tO. The ionization chamber has ereat 
precision or reproducibility (<0O.52) and can bes usedl vover 
a wide range of photon energies “C0 KV to 50) MV) and 
fluences Gice. cGy/min to ton ceGy/min) and various’ types 
Ola badtatlous. 2ncludanese Leceronsr an duneuitt0 14) eho Caawildecon 
anda hire rent Cy factor is needed. 

In our dose measurements a Capintec PR-O6C air 
Lonizatvion chamber was ased. Lhe characteristics or, tnis 
enamber Fare, listed in) Table 3.0: This chamber was 
connected to a calibrated electrometer (Capintec Model 192 
ADs Measurements in tissue-equivalent phantoms were 
Ovtsinede wate 82e 00.5 em | buddd=up “caps, Sinmeplace sn in 
line -eqguivalent: cork, ‘thes build-up cap) wass removed (in 
SOudermto observe the effects of electronic disequilibrium. 
Under these circumstances, the diameter of the air cavity 
WAS SeMM ys wlth al tadioloectical wall thickness of only 0).-5 


mm. 
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fabiesS.!<*Characteristics of the PR=06C air ionization 


chamber 
Sensitivity: 02200 ne/cecy Ctypical) 
Precision: 0.001 cGy, cGy/min 
Chamber Material: air equivalent plastic 
Active Volume: WI ose aviik 
Diameter: 7.0 mm 
Length: 22 mm 
Wall thickness: Oe 269m sem mey em: 


Se. lose thermoluminescent Dosimetry = (1LD) 

PLD Sis va solidestate: dosimeter available wins the storm 
of a loose powder, solid chip or rods or 
tet lon-impreenated rod “or flat dise.. A thermoluminescent 
phosphor has regular Crystal structure’, but when 
impurities are included imperfections arise in the 
lattice. Energy traps arise from these imperfections, and 
when the phosphor is exposed to ionizing radiation, many 
of the freed electrons (or holes) become trapped at these 
locations. Subsequent heating of the crystal will raise 
the energy of the electrons sufficiently to be released 
from these traps and return to stable energy states with 
the emission of light. If the light fluence is measured 
and plotted tas a9 function pol “temperature, the graph 


obtained is called a “glow curve”. There may be one or 


more maxima on the “glow curve” as traps of various energy 
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depths are emptied. The relative amplitudes of the peaks 
arenaneindicatvon, Of sthe relative populations Yor (trapped 
electrons. The total Jieht emitted durine pant or att vor 
the “glow curve", or the height of one or more of the 
peaks is related to the absorbed dose in the phosphor (4). 

TLD dosimeters are not absolute radiation “detectors, 
and must therefore be calibrated by exposing the dosimeter 
Eo known amountsesct radiation. Lhes, quantity. -disece.y 
measured by TLD is the energy absorbed by the phosphor 
i.e. the dose to the phosphor. This response may be 
related by calibration versus an ion chamber to the tissue 
dose in Gy. 

The instrumentation needed to observe TL is shown in 
Figure. 347. A heating element is used to heat the 
phosphor electrically while a photomultiplier tube 
Measures the emitted optical light. The photomultiplier 
current integrated over time is related to the absolute 
dose. TL sensitivity is defined as the amount of light 
released by the phosphor per unit of radiation dose, the 
lower limit depending on the type of phosphor and the TL 
reader. Most commercial TLD readers can measure doses _ as 
bow vas) 10) meGy. | The fupper limit of the usefwl ranges is 
eenerally lajmited by thes phosphor alone; LiF is linear up 
to 107 Gyox lll). Mor normal dosimetry applications, ae 
phospnor must have good light conversion efticiencygand be 
able to retain electrons trapped for reasonable periods of 


time at room temperature. The properties of LiF (Harshaw 
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TLD=700l) "eS the phosphor used in our dosimetry are given in 
Table 32. 2%. 

The (precision and reproducibility of the TL system is 
not) as “good as for the fon chamber. With best of care the 
precistom-sis 34. Thecaccuracy is of similar magnitude. 
However, the small size and useful large dose range make 
the TL dosimeters very suitable for many elLinical 
applications. TLD crystals can be made as small as 1 om 
or less. The powder fused into teflon discs or rods can 
be * used -/ron direct, insertion into body ‘cavittes andwtnius 
monitor the “patient dose in-vivo. TiD® jis wsie Gal in 
regions where the dose varies rapidly in space, as near 
iperstlmlal Or Mintracavitary sources  orvingpche)) build=—ao, 
region or penumbrae of megavoltage radiation beams. The 
vargee®iusetul’ jrange make it, applicable. for radiation 


protection as personnel monitors. 
Table *s.2:) 8Properties of Lif 


Dens ty : 220402 hen 
Effective Atomic Number: Dam a8e 2 

fe) 
Temperature of main glow peak: 190 C 


) 
Useful Range: LOSMieGy Oe Sec y 
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3.1.4 Tissue-Substitute Materials 

In order to obtain dose measurements in materials 
which can be related to doses in the patient, 
tissue=suabstitutes with “comparable. radiation “absorption 
and scattering properties are used. lf materials e are to 
absorb and scatter photons and electrons in the same 
manner as tissue, all of the. following, quantities must be 
matched: mass attenuation coefficient (U/0), mass energy 
apsorption, Pcoehiicient CL es electron mass stopping 
power Cs 0D electron mass angular Scattering power 
(02/0) and sass densttys..0 "G67 Jz 

Since tissues are composed largely of water which is 
also, (readily ‘available, this liquid is used routinely for 
measurements of absorbed dose. Lie O1Lr experiments, 
however, complex geometries and phantoms made of 
muscle-equivalent as well as lung equivalent materials 
were required. It would be inconvenient to immerse 
lung-substitute materials in water. Therefore we 
considered solid tissue and lung substitute phantoms 
rather than liquids. These can be cut imto various “shapes 
and sizes. The materials used as tissue-substitutes were 
polystyrene, prestwood, (for muscle), and cork (for lune). 
Tneswpniysi.cal parameters for these materials are Aistedwin 
Tablerm3..3). $/p is the average mass stopping power for the 
spectrum sor electroms setin motion by photons of energy 
eee 62 / pl describes multiple scattering of electrons 


travelling wa osmallepath in the scattering: material and 
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being scattered with a mean square scattering angle @ 
(es Us The densities (Oo ote cork were. U.s7a enone 
g/em?, and of prestwood 1.02 g/cm3. The relative electron 
densities Cp!) as measured by X-ray computed tomography 
were 0.52 0.28, and) 1.02.. The composition Vol “presewood 
and corkyis* variable and fauncertain. In order to validate 
the use of prestwood, tissue-maximum ratios (TMR's) were 
measured in prestwood and compared to those measured in 
wacersCsee section 323.) Tissue-maximum ratios. are a 
ratio of a dose in the phantom at depth d to the maximum 
dose in >the phantom at gdepth d> (54) £00) a) potur eae ea 
pexeds Mdistance. | trom (tives Source. Tissue-maximum ratios 
were measured for field sizes from 5x5 cm? Eo. SO0xs0R ecm 
for both prestwood and polystyrene. Values were found to 
beswithin be of those tor water. We therefore conclude 
that polystyrene and prestwood were water equivalent in 
their tTadiation ” absorption, attenuation “and | scattering 
properties at the radiation energy used in our experiments 


Ca-€.— SOUMV X-rays). 


Sal.) Tissue-Air Ratios, Tisswe-Phantom Ratios 

Doses along the central axis in water-equivalent 
media are wsually tabulated for different depths and fietid 
Sizeseas tissue-air-ratios (TAR) *or tissucs-phantom-ratios 
(TPR) “Lhese “are defined “as shown in Figure 3.5. 


TAR 


Dose (D)/Dose (D.) 
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EigseaqssoeDefinitions of tissue-air and 
tissue—-phantom ratios. 


TAR=D/D TPR=D/D_. 
a p 
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If the reference dose Mes is at the point of maximum dose 
inane 2 ipnanieom, thre J8TPRR tisiccablbed: a tissue=—maxinum—ratio 
(TMR). Hisitordca liste SuAR,’ s Piare used for Cobal t=60 
radiation and lower energy beams, but TAR's are less 
easily wdefined for higher Keray ) eneretess and “therecore 
TMR's are used for the higher energy beams (46). Once the 
EMRiYs for TAR 's stare *knowmy -the:sdosie at any Ipotat talone (gene 
central axis in the homogeneous phantom can be obtained by 
MULL cep Lyd euethe aelMR kCor: (SPAR): (atie thats depen = bre eee 


reference "“calibration' dose, dD. maximum dose in the 
phantoms (Cor dD» the tdose "to gaiesmaikl@ mass.) “of tissue?) in 
ac intr )s. 

wo idsisi tals or possible: tt o ase pavate (the) contribwti omietot 
primary and scattered raddia Giom by defining a 
Scatctier-airaratLorms¢SAR)) om ra csicat ter—phantom—naitito e aGSPRUr 
These itean® berrtcaleulated’ from the *eisisuesairlora tiascigor 


tissue-phantom-ratios in the following manner: 


SAR(d,W,) = TAR(d,W,) - TAR(d,0) 


SPR(d,W, ) DER Cdi Wiehe SOBER Cd 30) 

where: SAR(d,W,), SPR(d,W,); TAR(d,W, )>5 and es are 
the ee eres scatter-phantonm-, Gis's wea tr 
tiissue—phiantom—ratios atitdepth, id, sand for «field usiize We: 
TAn Caso cand abe R¢d JO) face the tissuesaita; tissue=phanton-— 


Tattoos, at. depthyed jmandeioreaepencil*beam@of “zero” field 


Si te 5 4.) 
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3.2 Summary of Calculation Methods 

The human body is not homogeneous in external shape 
onveineeinternalin densittys “Thesexternal eshapea cof ithe ibhody. 
eanebetaltened effectively@tbystusingtteither fabholisgicor 
external beam “compensators”. However, internal 
ditterences in density affect ‘the dose distripution ins a 
complex way which cannot easily be modified by a simple 
external “compensator”. In this thesis, we are concerned 
with the accurate calculation and measurement of dose 
within and beyond lung. We are specificallyvinterested ein 
dose calculations for a lower density object introduced 
into a water-equivalent mediun. Lis low density 
inhomogeneity will affect the primary and scattered photon 
RadiLatvon, as well as) the transport of “electrons, sete) in 
motion. WLthetimproved. accuracy in. dose calcudattions, the 
Severity of perturbation of dose distribution “by, ~tidscue 


inhomogeneities can be judged, and = Sethe possi bidiviy. of 


constructing appropriate “compensators” becomes viable. 

In order to calculate doses in inhomogeneous media, 
thes usual procedure, is to os first: calculate thepdocesin 
water, which is relatively easy to Viauleda tie 
experimentally, and then multiply it by an inhomogeneity 
Corpectrton Tactor, CLGeE )i 

Dose in heterogeneous phantom 


ICF = ee 
Dose in homogeneous (water) phantom 


The available inhomogeneity Connection methods are 


summarized sinvTable 3.4, in order of increasing complexity 
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tic i Csteewo metiOds Only Ccorrecu Lor tne seirect “Of Eue 
inhomogeneity on primary photon fluence. They are based 
on Calculation of a “water-equivalent™” or radiological 


thickness. Suppose the primary beam passes along a path 


Table 3.4: Inhomogeneity correction methods 


Metnod can take account or. 
Algorithm Longitudinal Lateral extent Electron 
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position OL Structure transport 


OLS Erle tu ne 


Rian OtaeacA ly NO NO NO 
Cor@Efrective’ SsD) 
(21) 


Batho YES NO NO 
C18 1095 


Lulu and Bjarngard 
(Batho) YES YES NO 
C67) 


Scaled TMR YE's YES NO 
GEhiswwork) 


DEeE trem emitsdranl! 


SARs YES YOR, S NO 
Gus. 20p 

Equivalent TAR YES YES NO 
CLO) EEO 

delta-Volume Vino YES (NO)* 
Cia) 

(Goynve) ilw ye sloin NOES YES YES 
C7) 

Monte Carlo ABTS YOR anne 
CLA 34) 


*under development. 
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of heterogeneous tissue with a spatially varying relative 


electron density, pic). Then the water-equivalent depth, 


d'= foscan. 


thes taissue-air Sratiow method sivess Jthe correction 
factor jas: 
TAR(d",W,) 
pi ee, (4) 

TAR(d,W,) 

where d' is the water-equivalent thickness, d is the 
geometric sthickness of material and Wa is the beam size at 
depo. oNiis COmLeCtiLOnmhact ObMaCCOUneS wf O DOL bie dam 
Size “and depths ot). the. calculation spoint, (butenot stne 
DrOominitty, of ther tahomogenel ty. to the caleulativon spout. oF 
the ‘lateral ‘extent of the inhomogeneity. The effective 
Sov “Csource-skinwdistance), method “is “equivalent. “fo. Vtne 
ratio. of PIAR*si*except that it uses percent “depth doses and 
aw einverse square correction instead of LAR s) (212 

A refinement which accounts for the proximity of the 
inhomogeneity, was developed by Batho CHS It SCL el 
Semi—empirical method, and its» implicit assumptions and 
derivation are not well documented or understood. It uses 
TAtioemandem Delative: selectron densities: and thereby 


indirectly. accounts for) the “changes “invithe, sceatterin: 
power” of the inhomogeneity. Originally this method was 
developed for calculation points beyond an inhomogeneity. 


It has since been extended to consider points within an 


inhomogeneity (109). Lhe sceneral. form, of the? Batho 
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[ tar(a,,w,) | ae 


ICF= Po (ey) 


1 a 
TAR W | 
[tar(a, .W,) 
where Po» P35 di; d,, and Wa are as given in Figure Shee ome 
inte method will be discussed in detail in section 3. o. 

A method which accounts for the shape of the 
inhomogeneity uses the volume integration of differential 
Scatter-adr ratios (doa) (13520). Ss thescorrect ton nict os 
is given as: 

TAR(d',0) + >_AS 
ICF = SS (6) 

TAR(d,W,) 

d 
d' is the water equivalent depth. LARKGd 4.0) yued oeee ine 
Zer0-area  (uissve—ainm ratio, and the integral LAs sums 
Scaucer CoOMperlputcitons from all” volumey velements in the 
irradiated volume. The contribution from each individual 


element is calculated as follows, Figure 3.7: 


“H,(a'-a) 


As =e [asar | ps { e (7) 


aera 
where: aSAR, 1s) the wditterential «scatter=adr Strat vomeror 
volume element, AV, assuming the surrounding medium is 
Witter etneceremalining@tetims correcthethiseevalue.} Lor the 
presence offs S5Urroundine structures vofdnonszwater density. 
a' and b"' are the water-equivalent depth along the primary 
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Fig. 3.7 warameters used insthe calculation 


Of Scatter cont niyDmWweLvon in) the dSark. 


asisumime —slingle scatter only. a and b are the “teal (or 
geometric) path lenghts. P 2 is the relative electron 
density in the scattering volume element, AV. u and 
Wu C8) are Sthe Linear attenuation coecifictents 0 wace.r sor 
ther primary “and the scattered photon (13). This method 
lends itself to npixeloby—-pixel Sleorrecti ons mein ii. 
COMpiketen eveetronedensity ) distributions. sobtained 1) VaGw 


scanning. This method was “ahead of its time” and was the 
first  tCo™separate the ) primary “and scattered radiation 
components in a heterogeneous phantom. However, all of 
Ene multiply astattered tadiation 16 scons deredmracu. aan diets 
behaves like single-scattered radiation; this. mays account 
fo cep OO ad ose results obtained Lot large PAsericd 
irradiations of heterogeneous media (10/7). 

The simple ratio of TAR's~ method (equation 4) has 
beem ~moditied by Sontac and Cunningham @l10) in® order to 
avsomacecount, for ~bhe. size, shape, and) Spos Ltilon so fest he 
inhomogeneity. This method is called the equivalent TAR 
method. An effective TAR(d' Wi, ) is determined where d' 
is the water-equivalent path discussed earlier, and Wi is 


aneecet Lective: 


bis = Dac RW BA hem cormtect it one & act Onmmes 


TAR(d",W1) 


LCF (8) 


TAR(d,W,) 

d 
deeismeasytOvObDtaitgnebut Wa PSsenore "der lcwltm tow 0 boa 
Since the entire environment of the calculation point must 


be considered. An equivalent field width is determined as 
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Eo) LILOWS © 


he 9S 5 (rau 
Wi nw (p)) (9) 


Pp. is the weighted average of all the relative electron 


densities in the irradiated three-dimensional volume. 


e 


o' = fe: CAN 2) vias ce yn zt) eV, (10) 


ime WCx,y,2) Lear set of weightings factors which, depend won 
distance of each scattering volume element from the point 
Gr (aleulLatid on. This method also lends itself to 
ape by pa cele COME ect LOlsT Geum re 

A method which relies on “ray tracing” between the 
calculation (poimt=™)-and “the scattering) sites has) ) been 
developed by |/Wouew G41), sand is referred to as) sthe 
delta-volume method. This method calculates the primary 
dose, and adds fimst—-scattier dose, corrected TOM 
heterogeneities. It further adds second-scatter dose, but 
this component is treated geometrically as first-scatter. 
Pinavlve ‘as vesivdualotermstor multiple (52) cearters doses 
added.) Itvis, very similar> to the volume integration of 
ditierential seatter-air “ratios, Dibte vd dst Pe tS denente beri 
OMLyr treats second=scatter as tirst-—scatter,, whereas vthe 
dSARo method treats all multiple scatter as first-scatter. 

With present day computers, calculations with the 
delta-volume method are long (e.g. a 3D calculation on a 
VAK=1/7 508 computers tforea l6xl6xl6 cm medium takess 1-2 


hours "per™ beam (141), but with development of very large 
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seale ‘integration *(VUSL) “chips” and/or array °®*erocessors; 
rowtine '“clintcal* -wse may be “possible inthe near feature: 
Furthermore this method has the potential to account for 
electron “transportstfior ‘hieh “energy X-rays) and ©is “seit 
being actively improved. Another method with similar 
capabilities, but Sowhich uses ta convolution principles ts 
being developed *in this laboratory (€71). 

The latter three methods are complicated and must be 
performed on computers. The “inapet": inftermatiou dis: fhe 
detaited Vanatomical (information “available tirom tC) *tsean¢e. 
Bebhorerithe adventivaf «CT ssceatning,~ this metatled anatompecal 
information was not available; internal body contours 
were estimated and this gave rise to errors in dose 
cakeulation whteh were greater than terrors ~tdue® $to." tthe 
simpitstic” “tnhomogeneity ‘corrections. per btse “G06 )- 0 tke 
present, however, detailed anatomical information is 
available, and the accuracy of dose calculations is 
limited by the quality of the inhomogeneity correction. 
Therefore, the complex methods must be validated 
experimentally and it must be determined whether their use 
truly improves the accuracy of dose calculation. 

The method which can, in principle, overcome all the 
deficiencies of the above methods uses Monte Carlo 
simulations (88). This method can "follow" individual 
photons as they are scattered throughout the patient, and 
Pan “ailso consider the energy 9 transpore op a charged 


particles from sites of interaction. In order to obtain 


143 


iad Ty ~ 
nae J 


_ocbese pats SOTae SRUGhE “eqagia” grer¥) nottadee 
Lanse) rehome ah SRG ae ESR 20H baste: ov 


thee of bee 9 OSeUAR Opens Cgth » 707 ‘preqanes2 © 2 
gis bedsen Smtynnk \ chaverqul” Giakiose genet 
simiveleg ee) tetows6a’ a" ereg (esi Jud -aei31i¢¢eunm 
(2%) $sweawodat Atus 21° baqultsvst ented 


s 


hesewilgecs Ge ebodiee seta? rea3es aay | : 


pavetad, “tuqgnt” eft \apsuqnes. oe poxrelzeg 
G ' eae 


* eldeitess woltansaiss Teéatnosaae beiiatek- 
. 

2oteh sede (pehenace TS Ic seevbe 403 s1atet 
* tas7o 

jefdelsava ton, esv aobsueretat © - 


Aa ieansgsie; 


e é ce 


.. «* @efy 486g @107 Sne “Sd001900 —BiRy 
- OG agwer 


* ¢942s70 made  7wraetg’ S19” Aahiy natin 
a mS we 


jebee os vee guePaoexses \yatenegomatal cisetiqete 
- pe eT 


eusogn? tena@eiesn teidesek tsteten epee 5 
is om 


pietovlen  tteh Se qhaaus7® si? bas ,@idat 


,gjpeeeqonatal edt —26 vattawp e072 <¢ 


id 


beds ° +909 gbovtagd. | watgeus, ods. pares 


ese ban e@ped >. botinver sb ee aed - ta ytiae 


ae 
. 


_etiseluaiae. enet te Thaemaae od3 
» ~Ge 


sit ile ewwsveve insgapetsg gi (aes 
site? @386F.~ AeSR | 
caulbiw iio) ee 
hen ,Agedseg “at hae 


reasonable results, however, millions of photons must be 
considered! ein ™ sordere Com\obtain. ~aceurate @dosemeresuits 
CL7piS4) saeWithepresent sday dcomputerss.uvlt i potiPleetakes 
hours @towdays.of computation to obtain ~acdosemdistri bution 
in a simple inhomogeneous mediun; this method is not ‘yet 
practiecad eons routine elinical use. For example, 
calculations for a LOxlO0ecm- —fieid of 6: MV K=raystawer —00 
Hours | tO “obtain  ~ardose accuracy of 0252 0nma VAM iy o0 
computer (private “communication —1.R. Mackie). The 
implementation of the complex methods for HBI is still in 
the future. We are presently restricted to the simpler 
inhomogeneity Corr ecitd on methods, with possible 
modifications, and realize that one method may yield good 
results» in sa“given ‘situation, but? poorer “resules ein 
another. This is the rationale behind the experiments 


described in the following sections. 
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3.3 Homogeneous Phantoms 

As noted earlier, the) “calculation of Sose nee 
heterogeneous medium is preceded by a knowledge of dose in 
a homogeneous water medium. The dose, Ds» £or yt ive 
situation shown in Figure 3.8.1s: 


Dose (D,) = ICF x Dose (D,) 


(A) Equivalence of the homogeneous phantoms to water 

im order to caleulate and measure correction) tactonrs 
in a heterogeneous phantom, various size cork slabs are 
inGvroduced. "into: 4 homogeneous mediun. For these 
experiments, homogeneous phantoms composed of slabs are 
most convenient. Thus, polystyrene and prestwood were 
used as the homogeneous material. Their equivalence to 
water and muscle has already been discussed in section 
een) met Cy weeds Leste) ce Nevertheless, tissue-maximum ratios 
(TMR's) were measured in polystyrene and prestwood for 
Eve ldwesizes, 4x4 cm (tow o0x60" cm2 and are civenv in slabple 
SoD ee IMRas for the @larger tields were” not. fpreviioucily 


avadtilablestor 6° MY X=rays. For £leld sizes > larcer= than 


Joo ecm 2s, mthere - Gs ay minimalm “changer tne uk with 
imcereasdue. field» “si.ze7, because Scattering ~from |e long 
range” is negligible. 


The ilMR valueg splotteds “abe di filerentedentic in 
polystyrene and prestwood are within 1% of those in water 
wp to depths, of 20 cm. Beyond this depth, they) are jwithin 
a OMNES ore These are shown in Figure 3.9 for a standard 
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Distance) of 100 cm. Similar agreement is obtained for 


TMR Ys at the other field sizes. 


(B) DMR*s at extended distance 

Several experiments were performed at larger 
distances from the source (SAD= 200 ecm) to achieve the 
larger field size. THe IMR “seat “this ‘textended distance 
aceyewithin» IZ of IMR's meastired at the normal) SAD ot 100 


cm (see Figure 3.10). 


GOD sbocatron wor cpaae at extended distance 

For larger distances from the source, the position of 
the depth of maximum dose dae in a phantom may be 
Cotre Gent etnemecdt wan SAD. or LOOM Ee me oizel. This ss ides Co 
the different scattering “conditions woiilche acount, suLe 
secondary electrons. Therefore the build-up to maximum 
dose was investigated for field sizes 10x10 cm? and 40x40 
ecm in order to determine whether there is a shift in cae 
at. asthe 9 larger SAD. the build-up curve was found toe be 
reasolapLly wupitormyto. Within 0.59%, from ta depthmot. iow cm 


tom 2.0) cm for’ both small and large fields — Theretoresmany. 


; and 


Value "between l.5 cm to 2.0 cm ds> acceptable as doe 


shifts within this range are inconsequential. 

Many of the dose calculation methods do not account 
forselecctronic disequilibrium in the mediums The eange for 
electronic disequilibrium will be longer in a medium “of 
lower density and the geometric location of d is 


max 


expected to be deeper (143). In order to determine this 
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range For cork, the dose build-up was measured for field 
sizes 1LOxlO cm? and 5x5 cm2. For these measurements, 
there savas no “build-up cap on the “air Lonizatvon prove. 
d for both fieldsewas found, to be 6.5 “em Vasyeshown ein 


max 


Ed cure v3)... 


COMMER Lecter of backsicatrer i neumater tial 

The effect of material behind the measurement’ point 
see of bene 1snored in / dose Mcalcutlatrons, SO;e that) deal 
backscattering is implicitly @ssumeda  tihes diftterence fin 
backscatter due to either cork or polystyrene behind the 
probe was also investigated. Readings were obtained with 
10.8 cm of either cork or polystyrene behind the probe for 
several field sizes. The difference in backscatter was 
Pound toppe pl.3% for the 20x20) em? fteld and Mess) thane. 
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3.4 Simple Heterogeneous Phantoms 

Radiation doses were measured in simple heterogeneous 
phantoms irradiated with 6 MV X-rays and compared with 
those calculated using several algorithms implemented on 
radiotherapy planning computers. ThesRAri©roniT eee i2 
computer calculates inhomogeneity corrections using the 
effective SSD method. The AECL TP-11 computer uses the 
modified Batho method GL09)) to Cale urate tissue 
inhomogeneity corrections. Furthermore,edose calculations 
were also performed using the equivalent TAR method as 
implemented at the Ontario Cancer Institute. All these 
methods were described in section 3.2. Dose measurements 
were ‘made with the Capintec PR-O06C air ionization chamber 
described “ine-esection) ~ 3.17 The probe was used with 
build-up = cap= in polystyrene, whivle Che build-up cap was 
removed for measurements in cork. The geometry of 
trradtatl one as Shown. eli Kigure so 12. The probe is 
positioned at a fixed source-to-probe distance (SPD) ealong 
the central axis. The distance to the entrance side of 
theicorksasection, a; is increased, thus letitting the entire 
Cork section. rise “throuch the poimtron mtenesteater. 
Measurements were taken for different cork thicknesses, t, 
geometric depths, ds and field sizese Wap as Distedaean 
Table 3.6. These geometries were selected in order to 
sample the » range (of lung ‘geometries encountered =) in 
¢linical radiotherapy. For the larger fleld size), an (Sad 


of 9200 “emiwas used. The inhomogeneity correction factor 
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is plotted versus parameter “a” as shown in Figure = 3.12 
both inside and beyond the cork inhomogeneity. Results 


are shown in this €ieure for fleld sizes® 5x5 “em*>) 10x10 


em?, and the large field of 35x35 cm’ and cork thickness 
t=7.070'. 2) (cme Measured values are compared to ICF's 
Cabeuwlated Susinge these ratic Sof. TAR ts (effective SSD), 
generalized Batho, and equivalent TAR methods. From the 
experimental values one notes an initial sharp increase in 


the correction factor which then levels off. ING Tela se LIL 


size increases, the, impoptance of = Scattered “radiatvYon 


Table 3.6: Summary of experimental set-ups 


Cork jthickness Geometric depth Field size 
Shed regi 13.4) cen oe 5, Cm 
LOx LO) em 

20x20 “cm- 


35x 3:5. (cme 


5.0 em 15 “em 5x. em- 
10x10 cm? 
20x20 cm 


35x35 “em 


Teoeem 5S em Secs) fehuk- 
LOs10) em 
20x20 cm 
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nebative to the primary, increases; «but the -inhomogeneity 
afbbects (the  primary,-radiatfon .~more than )the scattered. 
For this reason, experimentally a greater LCF is observed 
Por therssmaller fields. For, field sizes .lareer than ,l0«l0 
em2, the Batho correction factors agree more closely, with 
the. sexperimental. values, presumably because the Batho 
cornection sindivecthy (by .use of *TAR!sawhich, sinc bude Mee 
scatter) accounts in some way for the change in scattering 
conditions produced by the inhomogeneity. However, it 
underestimates suthe qaecorrect Lon, «factor Aboot hieinsidemand 
beyonudweork «aby 24..5 The ratio of ~ TAR ‘s.*overprediicus etre 
eonnectionsfactor invcorkeaby as much asyes2c near thagcork 
polystyrene interface for the 10x10 ecm? field. Beyond the 
eonk polystyxetteinterface this method predicts a constant 
ICF which overpredicts the experimental values. 

The equivalent TAR method performs better than the 
simple ratdo of «TAR's method, .but sis inferior ta the —Batho 
method in this geometry both within and beyond cork. Lt 
generally overpredicts the correction factor especially 
near the cork polystyrene interface where the 
overnoredicotLon | isis .27.67 5 fore gag 10RLO, vem? Bite dd. athis 
performance of the equivalent TAR method at 6 MV is not 
Gomes bs bewk) with is Eeported superior performance at 
CobaLt-60. EFEor Cobalt—-60,..acreement with the jexpenimental 
daiitaeees anwitiins2s ebothsinse ide. and beyond lune gGbom, LLOo. 
Lhe smeater discrepancies, for,6 (MV. Xorays, abe «sattnibuted 


to. improper calcqulation of .the contribution jof gscattened 
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Fadtation tosdose.s2Therweighting« factors, haver ~tittie: 
effect Onesthe Correction «factors: eebuteinaceuracies 10 
Ssteatter-alt=ratios greatlylaffectethe corrections | factors. 
Bnesedithtculty »sidecer int *obtainingsuscorrects pezetocsred 
tissue=airy Cor phantomje ratioss Practicallye theee Ace 
found «by* (extrapolation: “based “on small Eteldesizesdatac 
This methodology (54) becomes difficult to apply for high 
energy X-rays because electron disequilibrium develops for 
smaller fields. 

These experiments: at 6 MV support the suggestion by 
Henkelman and Wong (43) that pixel-by-pixel methods which 
tzmeate Sscatter §tontribationss -[rom/eseatterine elements 
independently, may be inferior to simpler methods which 
use a bulk correction for tissue inhomogeneity. 

Eor -thessmaditer field sizet off 5x5acm (>the =rattot tat 
TAR's (effective SSD) method as well as the equivalent TAR 
MethodRoverpredact correct tones tactorseebDOtny Wa tothe aod 
beyond cork. The Batho method again gives the best 
agreement, but exhibits some izreguilarities® The 


Correctdon»s factor is" now overpredicted in lunge, and the 


curverpota calculated correetinny yiactours crosses the 
experimental curve Bahortly= behind Jthe, dntertacespanad 
underpredicts the correction: factors sthereatter. The 


ecorrection- “factors ins lingithus) aresoverpredicted by ail 
correction methods to a greater extent than expected. For 
the smaller fields incident on low density materials there 


isgesoss?,of, lateralgeelectronit «transport CL.268 more 
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electrons leave the central region than enter from the 
sides) C143). Therefore a decrease in dose, and a 
decrease “in | the correction “factor occurs. (‘Thismettece 
becomes more pronounced with higher energy X-rays and is 
not handled by the correction methods discussed here. The 
above experiments were repeated for 15 MY X-rays and are 
included in a paper to be submitted for publication; the 
eptects of electron ‘disequilibrium “are davesticated 10 
detail at this energy. 

When the depth dis less than the depth of maximum 
Dudiid Sup in tissue, the econds E107 of ~electronic 
eqn Dri umeiOt thew bah On me th OdmdSmn Ivara tied. Lepod eeeies 
used with IMR's measured in the build-wp region, the 
Correction “factors ‘so ‘calculated scive  vwery anomalous 
results. Therefore care must be taken when correction 
factors are calculated in tthe Jnilttal paild-up | restonwm.07 
at the intertace “of two media. Bor Or MV eX—i7a creo iis 
region extends over 1.5 cm near the interface. For higher 
energy beams, there is a greater distance in which the 
Bathory conrectioneshould snot besapplied. 

The 5 cm cork-10 “em polystyrene combination. (shows 
Similar trends as) sshown in FP teure 3.12, except stiateen 
this case, the equivalent TAR method sli ghey 
underpredicts correction factors beyond cork. 

Por the 3 em cork=-i5 cm polystyrene “combination |) the 
measured and calculated correction factors are shown in 


Figure 3:43 Thess Batho. scorrection “eives the best 
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agreement both within andr beyoud cork. A ‘sivtvent 
underprediction occunstin cork ifor field sizes 10x10 cm 
and beyond. At slight overprediection oceurs for thaisx> 
em? field size. The equivalent TAR method overpredicts 
aeses, ginvelung  byierowghly, .2.52,. crosses oven unearnthe 
interface and underpredicts the doses beyond cork by 2.5% 
hoGwagthepesmaller™s of teldwesizes. This) (hanudidne) fom the 
connection fLactorsbeyondacork by the equivalent TARamethod 
PSUEnoOoIesekbirconsistenrt *gwilth= resubussobtained gtormoturear 
cork thicknesses. For example, it overpredicts the effect 
Hor eaerthickermccork section by soverestimating ucorrec tion 
fiactoust beyond (Fisuren3.l2a)5 but underprediects thevertect 
fox a oesmaller thickness. of “cork -by underestimating 
correction sfactors beyond GFigure 3.13) The agreement ost 
thes ereation.of (TAR Lsecaleulatedycornection (factousewittiea the 
measured values is not better than with the other 
combinations of cork-polystyrene. 

Similar cork-polystyrene combinations were used _ to 
Measure “correction factors at an SAD of 200 (cm and for 
ELela sizes, ups to-o0x60> cme. Theresiss not much change sin 
the measured and calculated correction factors from tield 
size 35x35 cm* to 60x60 cm?. This is expected since, the 
TMR | curve versus field size changes little beyond a 35x35 
emcatield. In saddttton) ay “0. cm \cork—3> ems) prestwood 
combination was also _ used. IMS Vexopat groin ston FeeWeh Sevalce Gacole 
thisscombination and a large Gteld size of 50x50 ~cem> are 


Shown deride re solo ofeoection 3.3. 
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The Batho correction performed very well giving 
values within 22 both within and beyond .conk. This is 
contrary to what would be expected based on similar 
experiments performed at Cobalt-60 (125,139) where the 
Batho correction underpredicted dose in lung by as much as 
LO, sinielarre= fields: The effective energy of the 6 MV 
X-ray beam ( 2 MeV) is not greatly different from the 
monoenergetic photon emissions of Lily and) 1.33) Mev £2 
Cobalt-60 beam. Therefore the improvement for in-lung 
dose calculations at large fields cannot be due to the 
difiierence in. beam “quality, .. In vorder “tovy expladnaments 
unexpected good performance of the Batho “correction, 
experiments were performed with Cobalt-60 radiation. This 
Usmpranalysed inrdetatl inwsecelon. 3.5, awhdtich conciicmom orm 


paper to be published in Medical Physics (May/June 1984). 
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3.5 Improved Lung Dose Calculations Using TMR'S 


(vs= TAR®S) “ine tne Satho Correction 


INTRODUCTION 
Dose calculations are usually performed by first 
considering the patient as being composed uniformly of 


water and theng applying a correction, to “account. for u0n 


Unit=—density #structures, such as lung. With computed 
tomographic (CT) scanning, the in-vivo geometry and 
density Ot the lung can be accurately determined, 


resulting in improved accuracy of calculated doses CLO)". 
The avaitability of “such JVanatomical, “iIniormation . nas 
prompted a re-evaluation of the accuracy of “ caleulatdons 
within heterogeneous tissue (21). 

The power law tissue-air ratio (TAR) ” method was 
developed adinittially for calculating dosevat points lying 


bo the “shadow of wide lung “trradtated @witth” Cobalc=60 
photons onl G2e) ee Tt haswsinece been refined to consider 


dosee points: within» lung, Cl09), Sand to account sore bune sot 


Linke s width 6C69>).. Application of the method to hichex 


photon energies has been attempted (113). A theoretical 
analysis of the Batho expression has recently been 
Penportced  G/0,139)" but emus E validations have been 


experimental, as summarized in Table 3.7. 


Recording to Sontag (LOO) ss thes Bathos —conrect boner or 


the simple geometry shown in Figure 3.14 is: 


sf) 
2 ji-P2 1A 
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for’ points beyond” ling, -and 


CF(TAR) = [ TAR(dy,r) 1°37 1B 


For points within lune. 

where: 

TAR(d. ,r) is) the tissue-air ratio at depth - and One 
Beeld wsige reat depen he ; 


fe) P, are the ‘electron densities (e/em™>) “of thea two 


53 
media, relative to water. 

These expressions assume that 

(a) the lateral extent of the inhomogeneity is larger than 
the incident radiation bean, 

Cb) "elé@ctronitc equilibrium exists: atthe edose: =polntcumot 
interest. the equation only “corrects for ~changes) in 
photon fluences and ignores the finite range of electrons 
set’ in motion. Ihis creates ditfticulties in ‘extending, the 
Batho method to higher photon energies such as 10 MV 
ee ray se C43.) 

Co)ssthes tissue ds) water—like in, its Jenersy absorption sand 
properties, 

(d) changes in backscattering from materials beyond the 


point "Of Anterest are negiigibie. 


Despite these shortcomings, the method ~continues® “to?” be 
used clinically’ with Cor “without )P -adjustments= for 
situations in which the above assumptions are untenable. 


Ageceen ine Tablewou , the Batho, Conrectilone factor nis 
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acecurabe sto within@.tingide and beyond lung for Cabale—o0 
irradiation with fleldusizges up to l5x)5 ucmi6gadoweverpeac 
the field size is increased, the accuracy degrades and the 
Batho “correction “underpredicts sthe “dose within lunge thy *eas 
mich Sas ElOZ fore titelds asedw im *halbti-pody irradiation 
Gl 5.)es The *proposed texplanattion is that scattered 
tadtationwiplays *aetgreater roleeat large fields eand Stne 
Bathe ~GAR@mnethod® fails o> account Georreculy |) foreyrnis 
component of dose. 

In ‘our “study with 6° MV X-rays, IMR” (tissue-maximum 
Eatios ) were substituted im the batho correction sand we 
found agreement to within 2% both inside and beyond lung 
for? [eld esizes “up “to 50x50 cm’. In orden to investiga ve 
whether the “ise “of “IMR "s “itstedd of " TAR.’ s “accounted ™ for 
the improvement in accuracy at the larger fields, we 
performed experiments with Cobalt-60 irradiation and used 
bocnaelaAR sw and) IMR eS@iwhich are both easily measured at 
this energy. Experiments confirmed a sient &icant 
improvement by 54 in accuracy when TMR's were used for 
inv=-Lungicorrections "for "all-field «sizes up “to “S0x50> Femmes: 
Tables Gof measured TMR data-for “Cobalt-60 are reported in 
section Vi. 

The theoretical analysis "of the “tmproved “*performance 
Ghnithe 'IMR=Batho )dcorrecticn the epresented Pein sanaly tic 


carcuiatious of primary and singly-scattered radiation sin 
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IIT. EXPERIMENTAL METHODS 

Doses were measured 1g) heterogeneous phantoms 
Coneistine — ot, islabse ot = Cork to ps samilate melume a dds on 
polystyrene or prestwood to simulate unit density tissue. 
The electron densities of these materials was determined 
Dy computed tomography tor ben0.2on 0.325 41.01.) Bande l.02 
respectively, relative to water. The phantoms were 
Peradiated “wi tn 20) chy X-rays* and ~ Cobalt—o0 eer aee 
Dose measurements were made with an air ionization chamber 
WLtne ane nacti ver “volume 61Of )) 70). 658 ml. connected to a 
calibrated electrometer. . Measurements in polystyrene 
and prestwood were obtained with a standard “build-up cap” 
GOL Oeeecmaeot 8 Lucite). while in cork, thesbudilda—-up scape was 
purposely “removed in order to observe the “effects, “or 
electronic disequilibrium. Under these circumstances the 
dfamteter ot “the air. cCavaty was 6 mm,. wath sa) plastic mwa 
thickness of only, 0%) mm. 

The geometry of the irradiation is shown in Figure 
Shey Beye iie=sprobe is post tioned wate alti xed ssourcesto-prope 
CSED oe dustance along thercenttal axis. For small field 
Sizes sup to 20x20 em- the SPD was 100) em for Oo MY (and) Gov 
em for Cobalt—60 radiations: the SPD was increased to 200 
CHL Ot ETO MY sand) oo... cm omom CObalt-C0Mmradd atiomaset o 


achieve larger fields. The independence of TMR values on 


a) Siemens Mevatron VI linear accelerator 

Deinerarronme —cUe and.) ‘—Theratron 780" <Atonic Biereye vom 
Canada) 

©) Gapintec PRO6C chamber and Model 192A electrometer 
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Fig. 3.15 Geometry used to obtain the experimental data. 


Dose at P 4s measured as parameter “4 is 


WeveiL@el saeieem 1 ji@ ID Ci 


SPD was validated to within 1%. The major variable in our 
experiment is the distance to the entrance side of the 
conk  'slab; as The “results “are most “easily “interpreted “by 
considering that the cork slab effectively rises “throwgh 
the ‘measurement point PP) thus allowing validation “of “both 
equations lA andwib «fer isodentric irradiations. 

For all field *sizges; “the phattem was larger® than. one 
P2ecd eeinw worden ® to esatissy “aesunptwon Ca) OLe thet nacnG 
equation. » Assumption” (b) is -gsenerally “valid, except «=f6r 
transition “pointes near the tissue interfaces. Assumption 
Ge) eis valid tco®within 22 Foro sope tissues “and ling sat 
megavoltage energies (54). The energy absorption 
properties of polystyrene and prestwood, the materials 
used in the experiments, were found to be similar to those 
ob water’ (54). Assumption Cd) is approximately valid= for 
our experiments since the reduction “in “backscatter 
produced when cork was substituted for polystyrene behind 


the probe was measured and found to be less than 24. 


IIIT. EXPERIMENTAL RESULTS 

Measurements were made at a constant depth of lo een 
alone “the” central axis below the entrance “surface of tne 
phantoms" "A cork thickness of LO *am was used”) to**sinulate 
tingevroek typical “thickness. *°Réesults "were obtained erored 
snavimGss>- cm) yermedium),(LOxl0e cm’) andes*large Piieaid 
(50x50 em)))4 The data expressed as “inhomogeneity 


Cortecthicneractcors “G2 )etor G&G. MV “sxA=-Tays are shown in 
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Piagure 3.176 A,B,C. The inhomogeneity correction taetor 6 
PpLObted as aufunction or la, the distance between the 


measurement point and the “top” cork-polystyrene interface 
(see Figure 3.15). The discrepancy between measured and 
calculated doses within and beyond lung is approximately 
Z6 OL smal las iwell (as. lerge stields. | This, tist'contrarye co 
expectations based on the sesults ) ote Van se Dyk aor 
itmadtarLonsewa theelamee Cobalt=60mt1e ldsiu(is Jp). 

simivar reswilts stor’ Cobalt=60 are shown in Figure ge1/ 
AB Cr The) calculated Batho correction factor using tthe 
usual TAR values and the proposed TMR values are also 
shown for comparison. For points beyond cork, there is no 
diftietencesbetween, the two correction, factoms, calculated 
using either TAR*s or <IMR's, as expected from their 
invatiant. ratio in equation LA. The maximum deviations of 
them “calculated correction ftactor f£Lrom measured, values™ is 
34. However; for points inside cork there is a marked 
difference between the TMR-Batho correction and _ the 
TAR-Batho correction. The deviation from measured values 
PSmeeCOnSuStenbtiy sworse “Lorethe TAR-Batho Correction, sand 
Canmoc masa tiaras 69. “for “large “fivedide. This poor 
DerLornance | Of. thes sTAR-Batho correction at large tields 
haswalready beens documented (125), but the simple remedy 
Off isine =~ the= “IMR=Batho correction has not been reported 
previously. Nevertheless, ft has» been noted by K-vuen 
(privates communication): (that for Cobalt-00 inadtlationu. van 
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would give better Sor een cart with experimental values. The 
feason “given was that TAR (0. 5cem) and MSF Cmissine scatter 
EACEOT) VCOLKeEGts ore ichangess wing scatrer behind the 
measurement point. 

They reasons for the better results observed with the 
TMR-Batho method will now be examined. Note that the 
original ~Batho correction was NOT imtended “tor fuse wicn 
TMR values which are suggested for higher energy X-rays 
(46,58). However, many clinics have (bldndiy) “adopred the 
Batho correction using TMR values for both manual and 
computer-aided radiotherapy planning. 

For (the, reader “s -conventence,, (tables of Cobalkt—60n IME: 
values for field sizes 5x5 cm* to 50x50 cm’, measured on a 
Pheratron o0fF and Theratron /S0+ are provided... “IMR  valwes 
asipameasured. on Doth units <form small cand. Varce @tdieddeiscivesc 
were identical to within 12. Comparinge our measured’ IMR"s 
witir those “calculated” from. published” TAR @s=) “wee tind 
similar agreement for the smaller field sizes (54) and the 
larger field size = ClS0)™= sor depths: * ups *to 208 cm. 
Differences are attributed to variations in the Cobalt 
source and collimator (Ciba a and the experimental 
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where De is the’ dose’ to a’ small isolated mass of tissue 
invainr. “We varenthustusing “a consistent set of) TAR and iin 


datasin our calculation of Batho correction factors. 


IV. THEORETICAL CALCULATIONS 

fn order to explain the unexpected performance” of the 
TMR-Batho method, we analysed how the Batho’ method 
approximates the dose in an inhomogeneous’7 medium. For 
comparison, we also calculated the dose due to primary and 
singly-scattered photons analy tically. ls Saye For 
Cobalt-60 raddation, the dose due to primary, “and 
singly-scattered radiatvon accounts for 93% of "the toatl 
dose” at ai depth of LO cm in, watem and) a, larce ficlapotelo 
CieeadiwGe lia Os). consideration of, only the primary.) and 
first-scatter components thus yields a reasonable estimate 
Of the structure of the solution. Consideration of higher 
order scattering is possible to second order s@l40py but 
beyond thas fan analytic solution, 1s)" intractable ms Glormsa 
finite heterogeneous’ medium). Phe vapplicatton ofr Monte 
Garlo techniques is" possible, although time-consuming 
(153). sHowever. this? approach doespnoteyield a structured 
solution which can be compared term-by-term with the Batho 
equation. Thewmelattere ts) pimport ant eet ten etc mem pied care 
linttattons of ‘the Batho tcapproach are to besidentitited and 


understood more clearly. 
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The total dose at P is decomposed into 
Die Dig tabue + Dot pte.) aD 
where D; = component of dose due to photons which have 
Scattered’ si times prior cto vinteractine at vrei nethe 
following derivation, we only solvesfor DD, and) D7 - 


Assuming monoenergetic photons of energy E incident 


fe) 
on a homogeneous” mediun, and electronic equilibrium at 


point P (see Figure 3.18), the» dose due” toy primary 


PadtacLons is atrnaw a liliy.s 


Dy OH eet d ao E 


0 p 0 
where d = depth in phantom (cm) 
@, = incident photon number fluence (photons/cm? ) 
(Cu, /P) = mass energy absorption coefficient (cm? /g) 
ilo =lLinear attenuation coefficient Geneon for the medium 


at photon energy E (MeV). 
fe) 
The dose at P, due to single scattering at elemental 


volume (QO also assuming electronic equilibrium, is: 


dA 5 
& d- ) fo} xs 
D, = toe vo(d-r cos§ ) do “'p n, dV e Har | ee (2) 
dr? p dA 
p 
where: = = differential Compton eros¢e "section per 
unit solid angle (cm /e) 
Ey = once scattered photon energy (MeV) 
n = electron density (e/cm?) of medium 
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(Cu, ,/P) = mass energy absorption coefficient for 
once-scattered photons (cm?/g) 
Hi = linear attenuation coefficient Cenk 
for firstsscattered photons 
ae, =element of area at’ Ps 
ASGMmMine, Sa eparallel= Scireular Tadtations Jfiietd and a 


spherically) coordinate systemy the total firse seatteredaoce 


at P from all surrounding elements rearranges as: 


r 
ayn fh Ti) 
D, = 2 mn, Soe Hod ee) sino 29 el \ Ho COS - yn) dr dé 
0 ae Ome a] do 


where: Le = maximum madivayil distance from 12 to the 


eeometric edge of beam, or, integrating over 1; 


r_( wocos =p) 
Sere m 
Ds = 2 ™N, dye Hod Ey aa sind do | ee i. 
0 doy (nijeose: =. ia) 


dé 


Consider now the homogeneous phantom as two separate 
layers of water equivalent material as “shown in Figure 
Sy RS in order to calculate separate Scatter 
Comtri pielonsam Lhe= first-seant er doses Da Peerenuine poddts & 


can be divided into two patus; ) Bh) the 


lw 
wer P ; EE A 

Eitet—cscatter dose at) P oricinating in Layer 1 De ie is 

the tirst-scatter dose “at Po “originating “in » Layer il. 


Dede restriction of the Limits of) integration toswiativin 


iavere tlesand) LL yields: 
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Picesselo eGCeometry used for calculation) Ofemi ns Gos cater 


dose Co; pointi-r. from two  separatemilay erse sand sl. 
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where Eee maximum distance to the edge of region I 
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where ro «Coma ximum distance from =P to the edge of rection 


TL. 
For compatiscon with the dose calculated by the Bathe 

method, we also require an expression for the scatter 

contribution ftom the second Layer when the first) Vayer is 

removed: 
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Suppose now the two layers are replaced by two 
difterent layers “of electron densities Pp») and © P3)> 
Felativestoe (water sas adin  \Figure®> (3414. Then, in the 
expressions (4) and (C519) POmmetLns t= scathie tad ose 1a 
distances in the layers must be scaled by the relative 
electron density appropriate for .each layer CS:0 5 
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tnlorder to vexpress dD, (d)) in terms ‘of the dose in water, 
Deelah ok we use a Taylor “expansion of the’ exponential 
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In order to simplify the expression for the anaiytic dose, 
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B. The Batho Equation 
For tthe, situation. depicted in) Fisunme 3-14) thes seacho 


correction factor using TAR’ s is given by: 
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where Dig = dose to a small mass of tissue in air, 
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where De is now the reference maximum dose, measured in 
phantom. The only difference between CF, (TAR) and 
CF, (TMR ) isiune the use of difterent. calibrations whicuieare 
based on dD, (no phantom scatter) or Die (phantom scatter), 
respectively. 

We can now proceed to calculate ns For 


simplitications again. we dividesthe dose wn) water si ntowlts 


primary component, Doe and .<its £Lirst— scatter “compomesty, 
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where Dads) and Dae are given by analytic 


expressions (4) “and (3). The dose calculated using (Chie 


Bathoscorrection tactor, equation (12)iJis “siven by: 
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Dp (TMR) (TMR) D (d,) 


B W wed 


=H ye <I BRACE py SA eet aetna 


eoeeee 


wea. (15) 


iain D, is approximated by a first-order binomial expansion, 


then 


u - = ss 
D,(TMR) = on(— "JE, J a7 Holp2(d, See ele P3 
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Komethe hatho correction «factor to ytleld a‘@dose iwhteh. ts 
analytically “correét “for “primary. stand tfirst scatters we 


require that equations 4(10)icand (Cl?) be equal: 


S + [et 
res etolo2(d, d.) p3do] D,(d,) : 5 ent} D3 y (dy) i sod, ( 1-3) 
u ul Diietenhers 
en en 
%5 (£2) Ey eo(-£2) Ep 
lac en try Sey aed pbe: es (2M) & yrosp! (d,) + sdltd 1}; 
eo. 03) >n/ hm Nee Oa owe ol eae) (18) 


For the “ease when the measuring point P Lies in water-like 
tissue “beyond lung (ie. P, =e lee Pp, * 1) 


D, (4, ) = o (1-02 )uo(d,-da) te (d 


lw 1) 


and condition (18) is exactly satisfied as shown by Wong 
and Henkelman (139). Noite! sthate un) Vthiss Pease. the 
reference dose, De or oh cancels and does not appear in 
mela tason who») « This suggests that agreement beyond lung 
between measured and calculated values would not depend on 
the, choice ‘of TAR “sor IMR s> as shown experimentadiy an 
igure) 3.17 . However, for polimts P within lung 
Gare Cos P5 = 1 and P, x 1 ) eat Onl ee oe) Cannot pe 
simplified. Win — @oyreal(ayie to verify the equality Dadi 
Daca. Deedes) and the dose at the reference point, Ss 
Ea must be evaluated. 


A computer program was written in FORTRAN 7/7 to 
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I ii TE 
calculate Di Cd, )/ Do ; Dy ,64,)/ 9 , +and Dy 64,)/ Oo 
analytically for some of the experimental conditions used. 
Results “are* ™ shown in® Table 3.9. Equation (18). can be 
rearranged so that we can calculate BE ae 


-1 


S_/P_ = ———— 
mm (1-93) 60000 


Fen (dredp) spade 
ee i CDV e a eNO GE) ey 


(19) 


“ud (1-p3)r4 
(1-ps) e7Hedn't-P2)f20) £ + cody to, oF (ds )/0 + os0!1(d,)/00)] 


Stee ean also be obtained trom a Monte Carlo calculation. 
Bor sa Cobalt—-60 beam and a field size of T0xl0 cms Sie 
rs, OF O59 (Chak. Mackie, private communication) .. wide he lsmet a 
very good agreement with the values as calculated by 
equation §C19)).  Wone (140) cives Sr = 0.043" for an 
tupinite. otield = size. Theretore, using, CheCIMR® we expect 
good agreement with measured values in lung. Lt CPE ra? 
had been used, Dd. would appear in equation (16) inetead or 
OF and see. would be used instead of Se dace im equations 
Gia sols). and).@lo) .. ther rieht hand sidesot equatd onm@la) 
would remain unaltered but would now equal Beds Ens 
expression was calculated by computer and its numerical 
Value is) given in the last column of Table 39. 
Lyprcalay. the layers for? Celectron, Fequil librium used an 
Cobalt—O60O is) only. 0. cmeson that dhe ise nNesili gible sand 


does not. give the correct numerical result. The correct 
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Table 3.29:  Analytically calculated doses at :P 
(Figure,3.19)tfor different «thicknesses of the 


Layers il .andsIl; 


ICAL an i 
Dasa si. (Org Dee Cae 20 Datla i eo 
lwool lw sl 4 bl ; calculated S 
(keV cm?/g) (keV cm? /g) (keV cm?/g) or Se. 
Fileld#diameter.= 5.6 %cm 
d, =l5cm 
d, =10cm 2.7600 One hU9 0.4213 OF OSy 2 
d, =l5cm 
d,=7.5cm 2080 Oso 2a On vow 0.0326 
d, =l5cm 
d, =5cm Dele Oa 0.9088 Le aS On ersysttt 
d, =l5cm 
d, =3cm 6015 1.42.95 126495 0.03338 
doen 
d, =lem Oe E70 PRN Ms ee) Dearie 050340 
Fieldvdiameter = Ll.2 cm 
d, =l15cm 
d,=10cm 4.0045 O27 642 eee 2305 0.0348 
d, =l5cm 
d,=7.5cm Be oO.) len Sewer) Pe9r3g 0.0367 
d, =l15cm 
d, =5cm DEO Oily), 25037.0 DOO OF 03:87 
d, =l5cm 
do =3cm 1.9040 2.8648 Seo a9 0.0394 
d, =l5cm 


d, =lem 0.9089 3.8599 4.0967 0.0372 
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Valuegtolusemin equation (19) is S/F and therefore 


TMR's should be used in. the Batho corrections. De and 
CF (TAR) would consistently underpredict doses within 
lung. Erom),) Table, 369 "one “cansobserve-thateasathevrield 
size issincreased the right+hand “side. <of. equations (19) 
increases in value. Soa also increases forpilarecer 
Fields. but S/o. remains negligible. Therefore, the 
deviation sof CF (TAR) from measured values should get 
lanperpas ihe field sizes is” itncreased: this is also 
observed experimentally. 

Thewiactethat CF, (TAR) is lower than CF ,(TMR) can 


also be deduced from equation 1B. 


CF ,(TMR) [TMR(d,,r)] : 


[TAR(d,,r)/TAR(d ,r)]°37} 
Zz m 


WAR Caye 1 2h nee ARG tian: 


i] 


CE, (TAR)[TAR(d_ ,r)] “78 

Bopelunc wel ps 2b and [TAR(d_,r)] >1. Therefore CF .(TMR) 
2 CF, (TAR) and since TAR(d_, 1) increases with increasing 
field size; (the inequality increases witheincrneasing sitelad 
size. Aftpartial physical explanation why CF, (TMR) should 
give better agreement is that; 

(a) the TMR value includes no inherent backscatter 

(b) the TAR value includes inherent backscatter. 

In tung of density 0.3 ¢/cem-, backscatter 1s, reduced and 


(a) is more appropriate. 
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V. CONCLUSIONS 

In this paper we have investigated the Batho 
power-law method for “calculating “dose within ) tune. 
Experiments indicate good agreement between measured and 
calculated dose “for “points “in water-like tissue beyond 
lung. If the measurement point lies within lung however, 
there are large discrepancies between calculated and 
measured doses only if TAR's are used. When TMR's~ are 
substituted, the agreement improves considerably for 
Cobale-60 irradiation. Agreement tsePalso eood for =o) = My 
X-rays where TMR's are normally used. 

The reasons for these different results using TMR or 
TAR's were studied by comparison with analytic 
Calculations of primary and =tirst—-sceatter | doce. Ere Seire 
reference dose calibration is performed in phantom, there 
is closer agreement between the Batho- and ava bytie 
expression. 

An improvement in accuracy for dose calculations for 
points within lung by as much as 5% can be achieved simply 
Dye supsthteutangm= (MRV Ss) @in® the? “Bathor “correct tons, factor 
instead of the customary TAR's. Ehis=1s  pacticularivy 


important for lange tieild trradiation, which is “now @used 


More routinely (32): 
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3.6 Complex Heterogeneous Phantoms 

fhe human thorax is by no means similar to the simple 
geometries discussed in the previous sections. The lunge 
are finite win width, separated by a mediastinal ssece1on .0t 
water-like tissue and may be smaller than the cadtation 
Pela The simpler inhomogeneity correction methods 
bestéd  =fand (discussed “in Vsection 13°4) do not tavemints 
account the lateral extent of the inhomogeneity (except 
for the equivalent TAR method). The effect of changing 
the lateral size of the inhomogeneity is now investigated. 
Since we are primarily interested in large field 
dosimetry, a field size of 40x40 cm* was chosen for these 
experiments. The experimental set-up is shown in Figure 
S225. A cork-thackness of 10 em was used) and this) “cork 
was “shrunken” laterally so that its cross section ranged 
from 30x30 cm? to 2x2 cm?. For the simple slab geometries 
discussed amine BSect lous 64, the Batho correctionwmecnod 
performed veryewell for small as well) as tance Pieldss and 
various thicknesses Of inhomogeneities. Therefore 
measured inhomogeneity correction factors were compared 
with those calculated with the Batho method. For a point 
5Semebelow the cork secttonye the correction factors “were 
plotted 9 for decreasing Bilateral size Sot (sthe vcorkm( (sce 
Picures 3.20). For cork areas of less than 320x300 cm-7 the 
measured ICF increases gradually and deviates by 74 from 
the constant correction factor predicted by “the Batho 


methods torea SOnS cm- area “of the cork inhomogeneity. for 
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cork’ sizes ‘smaller than 5x5 em2, thes correct ton, practon 
decreases, as is expected since it should return to 1.00 
when the inhomogeneity width approaches zero. 

Another type of phantom which was investigated 
consisted of two separated cork volumes on each side of a 
unit-density “mediastinum= along the ‘central axis The 
cork thickness was kept at LOD em, “whiylese tied 
eross-sectional areas were 20x20 cm2.” In this experiment, 
the cross-sectional area of mediastinum was varied from 
3220 cm2 to 20x20 cm2. This set-up 4s Shown in’ ~Fieure 
S21 ye Although, Chis) “geometry is) Still” rectilinearnw et 
resembles more closely the configuration of the human 
thorax. ine this case, thes Bathory methody "predicts. a 
constant correction factor of 1:0) since the, imhomogenecaty 
does not, lve directly abovemthe, pointe-of, calculation along 
the jcentval axis, sExperimentally, and tor both) sizecmuce 
mediastinum, thes corrections factor. (varies? rome Jee co 
OL99 (rieune 3.21) for all positions sof thew cork volume. 
This lower factor is expected because there is a decrease 
imescattered tadiationalatetally » ine) this ‘ceometre. LG 
there is only one volume of cork on one side of the 
mediastinum, correction factors of 0.99 were measured for 
distances of the cork between 2 to 4 cm from the central 
axis. 

in order to’ calewlate doses correctly an the above 
EWO Sioa tuat-Lons, only inhomogeneity correction methods 


Whtehmaccount) tor the lateral extent of the anhomogeneity 
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have the potential to be accurate. However, the simpler 
methods could and have been modified for this purpose. 
For the situation shown in Figure 3.23 a scaled IMR could 
be determined for an equivalent water depth and width. 


Then a correction factor can be calculated as follows: 


TMR(d" Wi) 
ICF = 
TMR (d,W,) 
d 
where d' and Wa are the water-equivalent depth and field 


size. scaled according to density €50))u and ad; Wa are the 
geometric depth and field size. The field size is scaled 
according toy density ini "both lateral directious. A 
simivarvapproachewas usedsinitiadly in athe ftorcmulattonmo. 
the equivalent TAR method (108). This factor can be used 
instde) cor ki. But. Lor calculation points beyoud conkeeweh.s 
method would not. accountator the distance sof the corksirom 
themcatculatrilonepotnt. 

Another possibility has been suggested by Lulu and 
Byarngard, ¢69),e and involves dividing Che® total yohantom 
intos constituents parts: ase shown in) @ieure S225 The 
imhomoveneity corrections tactorels calculated) as susualy tor 


phantom 4 which now has an inhomogeneity of dimensions 


equal) to’ the’ field. The dose at yi desis 


Thea copnectilonetactor for situation» Li is then: 
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Ke Do > D. + De CICF) 
TClS= = Def Doe 2 o 3 
eee $$ ———— = 1 -(D./D.)(1 - ICF) 
D Saez. 
Z 
TMR(d,W.)D(d_,W ) 
=l1- CL = ICE) 


TMR(d,W, )D(d_ ,W) 
iC¥ can also be replaced by any other ‘correction factor 
which assumes’ the lateral extent of the inhomogeneity is 
Bheater than! the, beamasize (GoGo, In our case we use the 
Batho correction factor for ICP {7 as it had given et hes pect 
Cesules inesimd lar sttuatitonus. 

tne Lulu and “Bijarngard correction | tactor ICF 2 as 
well as the scaled TMR factor ICFo¢ were tested on the 
phantom shown in Figure 3°23. For cork) sizes 10x10 em- geo 
35859" "'cm>) “in? a 40x40 “cm field, the correction factor 
using ICF 2 was within 4% of measured values, 
underpredicting ithe “correction (“factor within Scork) but 
Ooverpredieting le beyond.) as seen vin, Fioure  U325) pA oc. 


ve 


The ICF is plotted as a function of parameter "a" as shown 


ce eour ens or LOT meticCee, Bosam ple Neon kc sizes. The 
Correct onesactor ,caleulated@usins ICF oo is) within. 4 OL 
the measured values in cork. However, beyond scork.)” it 


predicts al’ constant ICF as does the simple ratio of TAR 's. 
This is expected since this correction method cannot sense 
ties, proximity "of the “cork to the point of calcudiation: 
Therefore (this method fs appropriate fom calculation) of 
doses within cork (or lung), but not beyond. 

Since the detailed equivalent TAR method accounts for 
lateral dimensions of the inhomogeneity smaller than the 


fieldmsize, the above situation 1s oné in whieh it) should 
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perroru very well. fudeed," “within ~lungs = scorrectrom 
factors are within “22% "of" "the “measured “values "(Figure 
Diet ie Beyond "lung, “and"for cork "sizes of = !O0x10 em sand 
smaller, discrepancies are larger. For comparison, 
correction factors calculated using the generalized Batho 
method are also shown in Figure 3.23. 

When the cork is "only 5x5 cm*\falloneethe central axis 
ia a GOx40 cre field > Figure 3923 9As irresuLarities occur 
and both ICF and ICF,, underpredict the experimental 
BCReby "as *much as -o2 "in thevitrst few centimeters of m@cork. 
in thts scace, most. of the “cork, thickness less behind. mtite 
probe. Mose* tof "the “material above’ -the probereisceot 
untiv=density which scatters radiatiion@and “sets Vellectrons 
in motion more than cork does. However, some cork lies 
directly above the probe and transmits more of these 
scattered photons “and” electrons” set" in “motion "than 
unit=density material. “This *cork (volume™ “is” tnsuffictent 
tO=- Achieve “electronic SVequilibrium, but elecetrous seuein 
motton in the surroundins wnhit-density material can = *reach 
thes®aprobe® -from’ Tthreevsides.*9 Therefore the siltiteatilouess 
suchethate ean Séexcess* "of electrons ~reaches  ? 7the seprope, 
thereby increasing the dose. The calculation methods 
witeh edo not “account for dose *dwe*to electron transport, 
do not sgive the correct “LG. 

For the mediastinum ™ phantom shown in Figure "3%21,; “the 


Correction stactor Mis =~plotted versus patameter ~a~s= for 


the two mediastinum sizes, the correttion factor” ” is ©0:97 
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equivalent TAR method “senses” this situation and predicts 
the, correct correction factor. ‘The Batho methodyas well 
as ICF . Wield gael ac HO hmOt 6 lO. ltiwili be sng erLrouenDs 
less, than 34 for mediastinum) sizes .5x20 cm ,ton20220 ema. 
These are the typical mediastinum sizes encountered in 
humany hemi-body eixnradiation.: The ICF o deviates by as 
much as -6%4 from the measured values. The effect of the 
cork) =<on) thesedosesatethe centtalpaxiswiseexacceratedpand 
thisebactor is not tusetulsiorsudose MealculLations | einemt ie 
mediastinum. 

inw conchusion) «then; form, vale these rectangular 
geometries, Ene ssimpleveortectionstactorsmoane Der uscamic 
calculate doses= beyond, within, and) “lateral | too lune 
inhomogeneities with an accuracy of 3% along the central 
axis. The only exception is for a small inhomogeneity 
along. the. central axis ina large fields “near thevtinst 
cork-prestwood interface, discrepancies of 8% were 
detected. This can be attributed to electronic 
disequilibrium above and to the sides of the probe which 
is not handled by any of these methods. Geometries such 
as inhomogeneities smaller than field size along the 
eentral axis and water-equivalent material along the 
central axis with inhomogeneities to the side will be less 
of a problem as the energy of the beam is increased (114). 


At higher energies the contribution to dose by scattered 


radiation is reduced, and is more forward peaked (i.e. 
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behaving more like the primary beam). Therefore, there 
will “be less%photon’ scatter from*points distant from the 
central axis. However, for higher energy beams, the range 
Of Selectronserseteein emotion As teereaten, aandweneccrore 
electronic (“disequilibrium may= develop; especiallye! tor 
small fields. Inhomogeneity correction methods which 
handlesscatteredsphoton radiation corréectly= twilleenote abe 
aitected Sby a change in the *dtirectt#onvof tthe scatter, put 
most of these methods do not handle electronic 
disequilibrium. Therefore there will be problems with 
higher energy beams at interfaces of two media where 


etlectvonvtransport isedisturbed. 


SUMMARY 

In the above experiments, various inhomogeneity 
correction methods were tested in a variety of situations 
similar to those encountered in radiotherapy of the 


thorax. None of the methods performed very well in all 


situations. It is therefore recommended that before any 
Ones method 4:6 “chosen” to™ ‘correct “doses » fore lung ss aite 
performance be verified in the situations likely to be 
encountered iIclinically. Moreover, if a method performs 
Well, atone energy Of radiation, it isi notstoOe be assumed 
that it performs well at another. Finally, the accuracy 
Of these calculation methods should "be “checked in jaa 
humanoid phantom. This] will be ther subject vor the next 
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3.7 Humanoid Phantom 


50 far, only central-axis doses in simple rectilinear 


phantoms have been investigated. The accuracy of the 
calculation methods for “off-axis” points, which are 
relevant for hemi-body irradiations “was “therefore 


investigated. The Alderson RANDO phantom (62) was used to 
simulate the human thorax. Tt is» constructed of "materniats 
simulating the anatomical structures and densities as well 
as the typical external shape of the human. 

If TMR's measured along the central axis are to be 
used for) off-axis ‘calculations, corrections, jus pe 
applied= for variations in photon fluence across they bean. 
For \Cobalt-60;5) this does "not “create a iwseriouss provwem 
because the primary photon fluence is very uniform 
(Corse eimih ce ise pe However, for X-ray beams from linear 
accelerators, the fluence "profile depends on. the X-ray 
target sm ithesinattenine filter, “and -cold@imnation.. hor stne no 
MV accelerator used in our experiments, the profile was 
measured at several depths in a water phantom irradiated 
with a 40x40 cm? field. Results are shown in Figure $8.24: 
tites path | lengths® of) primary rays will Sg bes -ditterent 
off-axis due to the divergence of the beam. However, for 
a 40240 “em* field at van SSD of 200 cm, and a depth of 15 
Cipin asphantom, there will be wa change dnuy pathtleneth or 
only 10.07 cm wp to the outer edge of (the beam.) Thus the 
divergence of the beam produces a negligible» ertect fox 


the: large fields at a distance of 200 cm: Thus, the major 
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diftferénees™ in’*tiuence* at “al @depth= are Sduceeto the 
characteristics of the accelerator and beam-shaping 
devices. 

In order to simulate the hemi-body treatment a field 
size of 40x40 en? at’ an SSD™ of 200) cmP incident, on ene 
posterior™ side of the? phantom was’ used. Cin’ practacerer 4 
pair of parallel opposed fields is used.) The 
posterior-anterior “distante™ "in. the: ~thorax= =) feetonm was 
bétween 21°55 to 23 cm. The Lateral dimenstons ranged trom 
BONecm at the shoulders toe 3255 "cm at Sthe waist. For 
calculation purposes, the field dimensions are determined 
by*the®sitzerof the field “or® the “sizer =of Pthe® phanren 
whichever is smaller (26). However, since TMR's change 
little” for field sizes rangine from “30x30 ‘cm tow 40540 
cm?, iit, is valid™ to use™ TMR”’ s¥ measured! for al 40x40) cn 
field) to caleulate doses in?’ the RANDO thorax. The 40x40 
em* “radiation field encompassed RANDO’ slices 12-26% with 
the central axis alone slice "16." ihe dnner? dimensions eas 
well as the densities in the thorax were determined from a 
CT’ Scan Of the phantom.” Slices’ 14 to 23” weres” scanned} 
with a CT?" siitvce’ ethickness® fand= indexing orm cine 
average relative lung density of 0.37 was determined by 
tracing the entire” lung. “The posterior-antertor diameter 
ote luney ine the*®mid=thoraxy slites* was in®the* range OFfFl35to 
16°°ems “while” © thePlateral size was typically 8” tor lO ecm. 
Similar values were obtained from human thorax scans “and 


are’ also*published by Van Dyk etal. ClZ 6.) 
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Several capsules filled with thermoluminescent (TLD) 
powder, were inserted in various places into the Rando 
slices number 15; WW and 18. The precision “and ‘accuracy 
of this dosimeter is 3Z and has been discussed in section 
3-153.2 The phantom’@was “irradiated and doses were 
determined. The positions of the TLD capsules are shown 
ip the Cl iseans of Pieure 3.25). 

TO verity the jcalibration, of the TLD) “dosimeters. a tne 
dose at qa in the geometrical prestwood phantom was 
measured both with TLD “and “with “the | tonization “probe. 
Brom. Chis. comparicou,s (as table of LMRins £oOt a> 0 <O meme 
field andy the known: profile lacrosse, the; beam, deses ame ai 
depths in a water phantom can be obtained (3). Doses were 
calculated for" each of ‘the’ “positions of the TLD 
dosimeters. Lf the sdosimeter was in Lune, “a correcelon 
factor was calculated using the modified Batho method 
GLOSo)n the BathoOsmeCoOLrrectwon  <aAcecOrding a tC OsmaaW Umea nid 
Bojarngard, and the scaled@iIMR methods descuilbed in section 
ery oe Two computer calculation methods were also used to 
obtain complete dosetdistributions in — the Pentire. RANDO 
slice. From these, correction factors were calculated at 
the positions of the TLD“dosimeters. The calculated, (and 
measured correction factors for several positions in and 
Outwoty lume are given invTable 3.10. (One of) thes <compater 
calculations was performed on the computer developed at 
the Cross Cancer Institute for radiotherapy planning (8). 


ANCT images@ot the Rando slice was entered. The lungs» were 
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Posterior 

Rando Slice |8 
lO Posterior 

Rando Slice |7 
Posterior 

Rando Slice |5 


Vite d.25 Positions of sthe TUD dosimeters are indicated 
im the Cli scans of Rando slices (a)51S ) Ubjal7. 


and (cy) 15. 
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outlined by joystick on the video screen and assigned Avigl 
average density of 0.5/) ¢/ cm 25 jihe doseudisine dpi coOlssuu 
lung is calculated using the generalized Batho Cowie ct vom 
method. such Tasdistribution 1s shown. inv Ficeurewos 20 0G 
Rando siice 16. The correction factors thie obtained ms are 
simitar to those “<calliculated manually sine ethessacno 
correction, as expected. 

The second computer calculation was performed at’ the 
Ontario Cancer Institute, using the detailed equivalent 
TAR method based on pixel—-by—pixel intormation available 
from a CT scan. All slices in the radiation field were 
scanned and used for a three-dimensional dose calculation. 
Such a dose distribution superimposed on the CT slice is 
shown, in Figure S.27 for Rando slice 138.) Formecompardsommea 
two-dimensional, calculation... wsing ondly ther central isiiiice 


and assuming symmetry in the third dimension was also done 


Gi 2a Gorrece1.on factors Cavcwlated using the 
two-dimensional) ‘calculation are | within, “17 of those 
calculated using a three-dimensional Calcule ta one 


rnenerore ime hemi—body “irradiations <a” twosdimensitonadl 
calculation “would “appear to be sufficient provided the 
CaLculationesia ce. Ls nots near (thee lungs apex emote tie 
diaphragm. 

The agreement between the detailed equivalent TAR 
values and those measured is better than +42. (The 
uncertainty of the TLD measurement, however, is itseln or 


this order). The simple scaled TMR method also gives 
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Fig. 3.26 The dose distribution in a Cl cican of Rando 
slice 18 is shown. The inhomogeneity correction 


was performed using the generalized Batho 


me tirodaal09))e 


fa 


Fie. 3.2) The dose distribution, in, a Cl sean, of Rando 
slice 18 is shown. The inhomogeneity correction 
was performed using the equivalent TAR method (110) 


as implemented at the Ontario Cancer Institute. 
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values within 42 0f the measured correction’ factors: This 
indicates that for the hemi-body treatment, doses within 
lung are obtained with reasonable accuracy using a 
simplified TMR method provided the correct internal 
dimensions cand densities; as obtained onvas GTI° stanze are 
entereds,) Ther detailed pixel=by-pixel correctionimethod’ is 
thereiore notwessenttalvfor accurate dose “caleulation tn 
hemi-body treatments. tne) Lula and) Biarngard) ebatho 
correction performs as well “as the Sesimple wscaled Sir. 
However, the generalized Batho method gives correction 
factors 2 to 3% lower than the other correction “methods. 
Ehasp is) ‘expected! since, 1 doesinot account. for the, finite 
extent of the lungs. Therefore if the generalized Batho 
CORLeci Onpeel Ges wil sed one can expect doses in lung in the 
hemi-body irradiations to be underpredicted by roughly 6). 


More accurate doses are obtained using any of the other 


methods. However, a pixel-by—pixel correction "method ) is 


not Sessential. The simple scaled TMR yields doses with 
equal, “accuracy. This simple method could easily be 
implemented for dose calculations on a treatment planning 
Computer. Then. “lsodosedispributions throughout the 
thorax can be obtained. The method can also be used 
manually to calculate doses, at specified) locations’ in the 
thorax. However, the accuracy of this method, as well as 
the other methods, depends on the accurate anatomical 
dimensions and densities obtained from CT scans. The 


influence of anatomical dimensions and densities on lung 
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dose calculations has been analysed in detail and 
discussed in several papers by Van Dyk et al. 
(123,124,125) sand jhas@ibeen. found sato sabe, important. In 
summary, the use of a simple lung correction calculation 
in conjunction with correct specific anatomical dimensions 
and average density has been found to be _ sufficiently 
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4. FUTURE WORK AND CONCLUSION 


This thesis deals with radiation-induced damage to 
lung injglargeslitield  pitrradtatcion:. The incidence of 
radiation pneumonitis asa function of. doses toygclungy vin 
radiotherapy patients. has been documented by Fryer (30) 
and Van Dyk (128). (Recently; fractionated!) radiotherapy, 
Which®¢isianore connon tim oma lieiteld einnadiations, £isealsa 
being used forthemi-body imradiations 9 (120). Chemotherapy 
may also be given to supplement the radiotherapy. Under 
thesevconditions; thesresponse ofelungwas swell) sase Stunord 
will be» different. Therefore dose response data for 
radiation pneumonitis need to be determined for these 
dtfiterent conditions. 

ConttLcosteroids vand santibilothes | have: been found 
useful potioeftreat- certain sgeaspects. of the acute radiation 
pneumonitis syndrome, but its treatment remains an area of 
active , research. When clinical fsymptoms- of radiation 
pneumonitis develop, the, patient experiences severe 
respiratory complications and the syndrome progresses 
rapidly and gis,potentially fatal. Entenvention “is suthen 


too, litthe gatoo) yiate” 2) fOnesof thertobjactivesisot 


often 
this thesis was to investigate a non-invasive test based 
on computed tomography (i-e. Cl scans) which couldas be 
used agean early indicator tof tiradtiation® damages tot tluugs 
Lorvseordenws ston? this itestatcsbe susefuly sit eshowld indicate 
accurately the changes occuring in lung, preferably before 
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clinical symptoms occur. Lt should also lbe ssupertor ongat 
least as good as other non-invasive tests in terms of 
sensitivity and specificity. The mouse lung was first 
chosen to monitor radiation response using CT 
densitometry. The advantage of using mice is that a large 
number of animals can easily be studied. When animals of 
the same strain (which are often inbred) and age are used, 
differences in radiation response due to age or _ genetics 
are eliminated and good statistics are obtained. 
Moreover, the radiation changes in mouse lung are well 
documented in the literature. The severity of the changes 
as a function of dose as well as the development of these 
changes with time are known. Before embarking upon these 
experiments, however, the CT densitometry was tested with 
large and small objects and was found to be reliable for 
thesstudy of mouse tlingese(sectionsm2 al 56 82.126 58) 2.82) 
Indeed, CT number changes observed in mouse _ lung 
corresponded to the time course of microscopic changes 
observed post-mortem. More severe changes were detected 
for those mice irradiated to the higher doses. Thus, dose 
response data can’ accurately “be “obtained “using “Cr 
densitometry (section 2.2.2). This test could now be used 
to investigate other aspects of the radiation response of 
lung. 

Most of the irradiated mice did not survive until the 
late ‘phase of ‘radiation damage. Data for this phase is 


therefore limited. Also it is not known whether mice with 
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less severe symptoms during the early phase are those 
capable of recovering. CT densitometry experiments with 
mice® arew ’now beitie®carried cutPby Dr. 9 ees) “Mitier@ot the 
Radiobiology” Croup*®tomSstudy =theY@later “phases™ of * lung 
damage. These mice are irradiated to lower doses and are 
expected to  “surviver® thel acute) ophase) Fofe damage. Dr. 
Miller eis alsotstudying theseffect (of radiations protector 
dfugs: Rats with®a lune tvolume approximately = “ten*® times 
larger than that©of ‘mice might also tbe “considered” (section 
DLS SO) 

Once CT densitometry was established as” %an “accurate 
indicator of radiation response in lung, the second phase 
of our experiments was to compare the usefulness’ of this 
fest tosother tn=viivo’ tests. For these experiments, dogs, 
with a larger lung volume, were more useful (section 2.3). 
CT densitometry was compared to conventional X-radiography 
and two nuclear medicine tests CSO CELGM eS e24 5500s The 
objective was to determine which of these tests was the 
earliest "indicator of radiation damage to lung. in} ethats 
respect, conventional eX—-radtogpraphy twas imihentoneto tall 
tests closely followed by the macroaggregated albumin 
Periuvstonettest. CT changes were observed two weeks 
earlier, and changes in radioaerosol clearance from lung 
was the earliest indiLleator sof? vraditattron =re'sponse 
However, ethe wltimate *soallsof these tests *fs to apply “them 
to radiotherapy patients. Often in these patients 


radiation eifects ©"are  isuperimposed on” "recurrent tumor 
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growth. This eiwiligebecrdtiticuloem togpdisctinguien berom 
radiation changes if a radioaerosol clearance test is 
being used. However if ssa §epre=in radvatton (W-scan ie 
available ss "the Location of the tumor can be noted. Then, 
on subsequent scans, a solid tumor can be distinguished 
from radiation changes and changes due to radiation alone 
can be quantified. However, *"for” ac *ditiuse® ttumerc; 
differentiation of radiattonwchanee "Strom Stumor. is emore 
diiaicit lt Nevertheless the merit of the CT scan as an 
accurate Tandseanly indicatortof® radiation ehange © imeelume 
isis not ito ibe underestimated ; Another advantage of the CT 
sciamilss thateletcrs ae qulickers procedure: EGLo2emins) faa 
compared to 40 minutes or longer for the radioaerosol 
scan. In the above experiments the dogs were sacrificed 
asi “soion® asm all ‘tests. showed “abnormalinmy:. Thus, the 
eonrelation ofthe tearly effects!’ withtlatern eirecvs during 
pneumonitis “could not )be determined: This could be the 
basis of a future study where there should be longer-term 
follow-up to the late phases. At present, however, the 
tests detecting radiation damage in the “prepneumonitis” 
phase  Vare rof* "lintted Sivyalue * since no therapeutic vagent 
exists *that will modify "the radiation® ‘reaction. On the 
other hand, there ® fis “no® incentive “to investigate such 
therapeutic tagentsi tac long tas sit (is "not tknown §iwhere sand 
when to intervene. The pathogenesis | Sof “iradiation 
pneumonitisigac wwell agittherdpeutic agents sto \emodify ithics 
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studies, an early test for radiation response would be a 
valuables stool < This test should be non-invasive so that 
it ‘can -sultimately! be lemployed= sin. thes follow-up of 
radiotherapy patients. 

The CT densitometry. test, <as wells <as “the snuclear 
medicine. testy \still» need to, «ber investigated  undersa 
Vanlety wots conditions. Clinical ,trialsvare also necessary 
in -patients*® with wand jwithouts tumor, sin ~ lune.) lhe nev 
imaging modality NMR (Nuclear Magnetic Resonance) may 
offer promise in the detection of edema in lung. Besides 
Lessimacing capability pjeNMRahasethe, potenttal® toyaedetece 
inaviyo (chemistryn duceeto Biteqgability eto identity etene 
structure and properties of complex molecules (138). NMR 
cane thus enot yonlys aldtinediagcnosissebut maywdetectwthe 
molecular changes’ leading to radiation pneumonitis and may 
thus *®@point, to, ethe) undertyings cause and process of gthe 
syndrome. This methodeappliedstogmedicine isestitivingits 
imtanecyoeand its clinicaljvalue vandgpotentialsremain toupe 
explored. 

Finally, «in. thiseywork, the} «usefulness, sof) the CT 
densitometry test was investigated) in,clinical situations. 
CT scans of radiotherapy patients irradiated with a high 
Vocal yudose to Lung, or a lower dose to the entire lung 
(hemi-body irradiation) were analysed to obtain density 
information CoOL diagnostic purposes (section 2.4). 
Averages densities, obtaineds forse our patients' lungs 


excluding the diseased regions agree well with the 
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publishedivaluese(129)25 In thespatients (recefivinggavehten 
local dose to lung CT number changes indicated radiation 
changes and/or tumor volume changes. Ine thesSpatients 
receiving hemi-body irradiation no density increases were 
observed, <assexpected “for tdosevidelivered COCEGCyS <insas 
fractions) > Thiseestudyesis ongoing @with a new group sor 
patients irradiated to higher dose levels. 

Although "padiation pneumonitis will mot be, *ratal if 
only “sa small, Epartew cof Sethe telune |) Ys Sinradiated those 
patients who have a smaller volume of lung irradiated 
receiver "much higher Hocalized( doses. Jln vthesetcaces,, Sthe 
area of irradiation is well defined; CT scan can be used 
toy lobtainy pdose °eresponses datas atl) such sites. ein sien 
cases, an accurate dose distribution should be computed 
usinoe Gl eeseans: FOLLY 252 30 This data can also be 
Obtained for fractionated radiation. Li accurate” Fdese 
response datas ™are® availapletaionewall® these @eclinical 
situations, then dose to lung can be increased to achieve 
amegreater “therapeuticw “gailne) “This presupposes of scourse 
that dose canbe accurately “calculated candy delivered? to 
lung. Chapter 30°0f this thesis therefore ‘deals with line 
dose calculations. 

The existence of4a Lowedensi ty tmatertalb=<such @as lung 
in unit-density tissue alters the dose to points within 
andubeyond lung*inta complicated ways" elhe primary vasewell 
as thevescattered (Gradiathon®vis antected. Lung dose 


ecnliculatronsalgorithns@ane often! *crudepfatvaccounting® ony 
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for primary changes in dose caused by the inhomogeneity. 
However, inaccuracies in delineation of inner anatomical 
contours and densities can introduce greater errors. With 
Chiscanning ‘the *possibility ‘exists <to toptain™® the @icorrect 
internal geometry and densities of tissue structures. The 
accuracyof lune dose calculations 1s Inowtlimited® by wetne 
accuracy ™*ot thetcaleutlation “algorithms and much 9résearen 
Yeebeingedevotedatonimprove, Ehese &@2:eo9 7 si VO Sian) 

ny tthis ppthests, svappl ication a of wethe lung dose 
calculation methods to) p elas e:9 (ike td Bair padirat omnis, © ais 
discussed. “The simple tcalculatiion algorithms easiewe Mineas wa 
More =tconplex i“ wpixel-—by-pixel” “dose scomrectiiom method iwas 
tested. Depending upon the type of radiation and _ the 
complexity of the phantom, these methods may perform well 
im one vsituation bute not tinvanether. Electron atransport 
becomes a problem with higher energy beams. None of the 
simple inhomogeneity correction methods account Lor 
changes» inateliecttrom transport ath sthe tintertace of stwo 
media. Monte *4Garlo ) Gand ® nay >eracine Cor » «<convol ution) 
complex computer caleculat fons can overcome this 
limi tat Bot, 

Methods such as the simple ratio of TAR's (effective 
SSD )> €¢20) tewhich@Peonlyetaccount*® fonvechanves sain. primary 
radiation are mot sufittietently “accurate (seven " invesimple 
situations (section 3. aay) « Methods Pewh ron vindireculy 
Aceountnton (change tim Iscatterteuch "as the aba tho s*method 
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pixel-by+pixel.correction methods (equivalent TAR) any the 
simple geometries, giving agreement with measured values 
within»and beyond»: lungwto 22%.. This» method; however’, does 
not “account "for a’ situation» where: the lateral extenteort 
thes inhomogenei ty)is? “smaller “thane, the Pradiacionee fieda 
size, and does not perform well in these more complex 
Sadtuations~ For such cases the *olxel=-by=pixels 
equa vatents PTAR@method (107 7110) senses? this sttuation and 
gives better agreement with measured values. Lf thes Batho 
méthode is modifiedvas proposed by) Ludiuw ands Busrneard) (oD 
tor account) for the) lateral extent” Jot the inhomogeneity, 
agreement is within 44 for complex situations (section 
3.6). A simple scaled TMR method developed during this 
work, gives agreement to within 3%. 

Por, ther: 6yMV .X=rayse used) in, our experiments: the 
region of electron disequilibrium extends over 1.5 cm in 
ula tedensaty: material.) Thusy + there wile bess inaecutractes 
fife dose calculation, “near® the siutertaces fof) twotmediar 
However, this’ region is not are@el enoughs Comwarrant weer of 
Momte! /Carlo® or ray) tracineymethods.*s Atty present, use, OF 
these methods cannot be contemplated sey hemi-body 
irradiations due to the time involved (section 3.2). 

Finally, se dose! Pcalculationsy were performed “for 2 
humanoid phantom (RANDO). This phantom was irradiated 
usine the same “set-up “as for the hemi-body patients 
(section 3.7). Calculated doses were compared to measured 


doses. in this» complex configuration, dose= computations 
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performed using the equivalent TAR method, the Lulu and 
Bjarngard version of the Batho method and the scaled TMR 
all gave good agreement with measured values within and 
GUESOL Lune (+4). eb orahkl these mecalculatrons. anatomic 
information was obtained from CT scans. With this correct 
information even the simpler inhomogeneity correction 
methods yield § dose vicalculatitons = wrth 842 accuracy. im 
hemi=pody trraditatiomn sr. 

The work) reported «ins schise ethesis, has validated 
computed tomography as a Valin ble lool) scorn woot 
radiotherapy planning and diagnostic follow-up. Coe tire 
protocol “for the = hemi=body™ irradiations “at, the Grogs 
Gancer Institute the total. dose to lung was kept at 6 Gy. 
This dose was prescribed to mid-lung assuming unit-density 
Pore Milne. such ae calculation, where tie density. of “lung 
is not considered will underpredict the dose by more than 
DOV rate some: Locat Tons: inmate me WesCo toe C ClO nm onidC GO hammane 
mid-lung in’ RANDO Slice 16 1s 1.073 calculated using the 
equivalent TAR method. The dose at this position would 
therefore be 6.44 Gy as compared to 6 Gy prescribed to 
ties potut, py =the 7CCl "hemi—body (irradiation. protocol, 
Subsequent to the present work, dose at any point in’ lung 


can be calculated with 4% accuracy. 
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PROGRAM SINGLEFOR 
Elien  Et-Khatib 


June#L7/33 


This program calculates the dose due to first scatter 


from two separate slabs of a phantom. 


real*4 pi,rsq,eleens hn lm. xl ez tieid TAND COS) 6 END 
Parameter(PI=3.1416) 

Parameter(elecm3=3.343E23) 

Parameter(rsq=/.941E-26) 

type*,'Please enter primary linear attenuation coefficient 
1 Cem) se 

accept*,xmuo 

suml=0.0 

sum2=0.0 

type*,'Please enter primary energy (Mev):' 

accept*,hno 

alo=hno/0.511 

type *,' Please enter the field diameter (cm):' 

acecepy =, fieid 

type*,'Please enter distance from surface to probe xl (cm):' 
Accept ;xs 

type*,'Please enter thickness of lower layer x2 (cm):' 
accepE*;x2 

type*,'Please enter electron density relative to 


1 water of second layer' 


er seretuste' ae rgorg aiaT 


~~ we 
eepebdy w 16 dtule eea0hqea” ow2,m6T? 
| 7 


tw)», eel, Senne, wey, te s*fena : 
(olbi. Ealtisevennte® - 
(eva At. fofurel 2) Ie) sas7e4: 
(Os-a50e pes eetomarat 

(* teentd Svawisg 16545 n06elT!,*sgys 7 
“((l-aa) Hy 
suse, *9498a8 _ 
bseiade | 

0, 0-aue 


54 


©) Gene Opamtd§ teree: sess kd! canes | 
odiy*ognao8 : 
£480 aadwe ne 
or), Vedenbt Gleb? adt Seated erate rie 


ie on 


enh 
— - - rons & 
ios) (4 4089 91 eoaleee worh aawenth a iF pimees ia’ s ao 


Perwt 4 


, 


‘kee? Qu fawqed 
; | 


238 


accept*,ro3 

do .a#050 9179 50 50801 
ec=COSD Ca) 
al=alo/(1lt+alo*(l-ct)) 
hn=al*0.511 


hn is the energy of the once scattered photon 


DitrerentialmProebabilkity 
cem=l-ct 
ep=lt+ct*ct 
tem=cp*1.0/((lt+alo*cm)*(1l+alo*cm) ) 
tem=tem*(lt+alo*alo*cm*cm/((ltalo*cm)*cp)) 


dsig=(rsq/2.0)*tem 


integral “(rio bapi lt oy 
temp=(€2* (l+al)/ (142441 )-ALOG(14+2*%al)/al)* 
1. Gli-ral)/Calal) 
temp=temp+ALOG(1+2%al)/(2%*al)- 
l Gis <a yy Gla 1) (Gls 2 als) 
linear attenuation coefficient of the once scattered photon 


xmul=temp*rsq*2.0*PI*elecm3 


imceeralvenersy stransler coetitiicient 
temtr=2.*(1+al)*(l+al)/(Cal*al*(1+2*al) )—(1+3*al)/ 

1 GU+2 tal) / Gl+AZ al) 
temtr=temtrt+(ltal)*(1+2*al-2*al*al)/ 

1 Galea Cl+2 41) * (lao xa) 


temtr=temtr—-4*al*al/(3*(14+2¥*al)*(1+2%al)*(1+2¥*al)) 
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temtr=temtr—((ltal)/(al*al*al j-l/(@*al +i 


1 (2*al*al*al))*ALOC(1L+2 <a) 


Mass energy absorption coefiltcient for scattered photon 


xmulabs=temtr*rsq*2.0*PIl*elecm3 


evaluation of integral 
y=xmuo*C0OSD(a)-xmul 
TECATAND CEt eld) (2. 0* x2) )..Givea) sthen 
maximum angle for first slab is smaller than 180 degrees 
determine whether this angle has been reached 
TECATAND GEieldy/ (2 70*x))).GEaa) then 
lm=(xl-x2)/COSD(a) 
else 
Im=((field/2.0 )—x2*TAND (a) )/SIND(A) 
end if 
x=xmul*ro3*x2/COSD (a) 
if(y.-eq.0.0)then 
r=0.0 
else 
r=(EXP(1m¥*y)-1)/y 


end if 


Evaluation of swl 
Su=hn*xmulabs*SIND(a)*dsig*EXP(-x)*r*0.01*(PI/180) 
suml=suml+Su 
end if 


LECATAND Cf teld/ (2. 0*x2)))) CEs a) then 
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rm=x2/COSD(a) 
else 
rm=field/2.0/SIND(a) 
end if 
if (y.eq.0.0) then 
st=0.0 
else 
st=(EXP(rm*y)-1)/y 


end if 


Evaluation of Sw2 
ad=hn*xmulabs*SIND(a)*dsig*st*0.01*(PI/180) 
jot Carb S (Cal 00s) merce ern Orel erat tem 

type*, shn=", no 
type*,'dsig=',dsig 
type*,'xmul=',xmul 
type*,'xmulabs=',xmulabs 


! 


type*,'1lm=',1m 


’ ' 


type*,  rm= 


types. .50—"_ ou 
type*, ad=",ad 
endif 


sum2=sum2+ad 


end do 


swlflu=suml*2*PI*elecm3*EXP(-xmuo*(x1l-x2)) 
type*, This, 1s first scatter dose per 
rf unit flux from first. layer’ 


type*, sswitlu=*,swltflu 


erliv(yPavit2a)ega . 
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sw2flu=sum2*2.0*PIl*elecm3*EXP(-xmuo*xl ) 

Eype*; "This is the first seatter dose 

iL per unit flux from 2nd Javer 

type*, 'sw2flu=',sw2flu 

type*,'Please enter primary energy absorption ‘coefficients’ 
accept*, xmuoabs 

sp=sw2flu/(xmuoabs*hno*¥EXP(-xmuo*x2 ) ) 

types, lhnis: 1s) the fivsteacatterm toa pramarty ratlo sor she 
‘ second layer' 

types)  Ssp=* ,058 p 

Sitiop 


end 
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